vortex.h

Go to the documentation of this file.

#pragma once
#include "objectshape.h"
namespace GOAT
{
    namespace raytracing
    {
        class VortexPlate :
            public ObjectShape
        {
        public:
            VortexPlate();
            VortexPlate(ObjectShape& os);
            VortexPlate(VortexPlate& c);
            VortexPlate(
                const maths::Vector<double>& P,
                double r,
                double h,
                double dh,
                int m, 
                std::complex<double>  n,
                double r0 = 1.0,
                const maths::Matrix<std::complex<double> > alpha = maths::CUNITY,
                const maths::Vector<double>& Ex = maths::ex,
                const maths::Vector<double>& Ey = maths::ey,
                const maths::Vector<double>& Ez = maths::ez
                 );
            double radius() { return r; } 
            double height() { return h; }
            void setRadius(double r); 
            void setHeight(double h);
            void setr0(double r0); 
            void setm(int m); 
            void binWrite(std::ofstream& os);  
            void binRead(std::ifstream& is);
            maths::Vector<double> norm(const maths::Vector<double>& p);
            bool isInside(const maths::Vector<double>& p);
            double volume();  
            void initQuad(); 
            maths::Vector<double> calcCoM() { return P; } 
            maths::Matrix<double> computeInertia();   
            void setPos(maths::Vector<double> r) { P = r; }; 
            void setPos(double x, double y, double z) { P = maths::Vector<double>(x, y, z); };
            bool next(const maths::Vector<double>& Ps, const maths::Vector<double>& K,
                maths::Vector<double>& pout);
            VortexPlate& operator =  (VortexPlate& f);
            VortexPlate& operator = (VortexPlate f);
            void scale(double sf); 
            int order() {return m;} 
            double vortexHeight() {return dh; } 
        protected :
            double h = 1.0; 
            double r = 1.0; 
            int m = 1; 
            double dh = 1.0; 
 
        };
    }
}