ellipsoid.h

Go to the documentation of this file.

#ifndef ELLIPSOID_H
#define ELLIPSOID_H
 
#include <complex>
#include "vector.h"
#include "matrix.h"
#include "objectshape.h"
#include <iostream>
#include <fstream>
 
namespace GOAT
{
    namespace raytracing
    {
        class Ellipsoid : public ObjectShape
        {
        public:
            // Ellipsoid &operator = (const Ellipsoid& e);
            ~Ellipsoid() {};
            Ellipsoid(); 
            Ellipsoid(const ObjectShape&);
            Ellipsoid(const Ellipsoid& E);
            Ellipsoid(
                const maths::Vector<double>& P,
                const maths::Vector<double>& r,
                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
            );
            void scale(double sf); 
            bool next(const maths::Vector<double>& Ps, const maths::Vector<double>& K,
                maths::Vector<double>& pout);
            Ellipsoid& operator =  (Ellipsoid& f);
            Ellipsoid& operator = (Ellipsoid f);
            maths::Vector<double> norm(const maths::Vector<double>& p);
            bool isInside(const maths::Vector<double>& p);
            double volume();  
            void initQuad(); 
            friend std::ostream& operator << (std::ostream& os, Ellipsoid E);
            void binWrite(std::ofstream& os); 
            void binRead(std::ifstream& os); 
            double a() { return r[0]; } 
            double b() { return r[1]; }
            double c() { return r[2]; }
            maths::Vector<double> getr() { return r; }  
            void setPos(maths::Vector<double> r) { P = r; } 
            void setPos(double x, double y, double z) { setPos(maths::Vector<double>(x, y, z)); } 
            void setr(maths::Vector<double>& r); 
            void setr(double a, double b, double c); 
            void seta(double a, bool VConst = false);
            void setb(double b, bool VConst = false);
            void setc(double c, bool VConst = false);
            void setr0(double r0);
            maths::Vector<double> calcCoM() { return P; } 
            maths::Matrix<double> computeInertia();   
 
            // protected :
            maths::Vector<double> r;   
            maths::Vector<double> r_2;      
            maths::Vector<double> P2;      
        };
    }
}
#endif