box.h

Go to the documentation of this file.

#pragma once
#include "objectshape.h"
#include <iostream>
 
 
 
 
namespace GOAT
{
    namespace raytracing
    {
        class Box :
            public ObjectShape
        {
 
        private:
            bool isOctree = false; 
 
        public:
            Box();
            Box(const ObjectShape& F);
            Box(const Box& B); 
            Box(const maths::Vector<double>& P,
                const maths::Vector<double>& d,
                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);
            Box(const maths::Vector<double> bound0, const maths::Vector<double> bound1); 
            ~Box(); 
            void binWrite(std::ofstream& os);  
            void binRead(std::ifstream& is);   
            void scale(double sf); 
            double distance(const maths::Vector<double>& p, const maths::Vector<double>& k); 
            bool next(const maths::Vector<double>& p, const maths::Vector<double>& k, maths::Vector<double>& pout); 
            maths::Vector<double> norm(const maths::Vector<double>& p); 
            bool isInside(const maths::Vector<double>& Ps); 
            double volume() { return abs2(d); } 
            double getSize() { return abs(bounds[1] - bounds[0]); }  
            void initQuad();  
            void setr0(double r0); 
            void calcDiag() 
            {
                diag[0] = bounds[1] - bounds[0];
                diag[1] = maths::Vector<double>(-diag[0][0], diag[0][1], diag[0][2]);
                diag[2] = maths::Vector<double>(-diag[0][0], -diag[0][1], diag[0][2]);
            }
 
            void setPos(maths::Vector<double> r) 
            {
                P = r;
                bounds[0] = P - d / 2.0;
                bounds[1] = P + d / 2.0;
                pul = bounds[0];
                por = bounds[1];
            }
            void setPos(double x, double y, double z) { setPos(maths::Vector<double>(x, y, z)); } 
            maths::Vector<double> getD() const { return d; } 
            void setD(maths::Vector<double> D) 
            {
                d = D;
                bounds[0] = P - d / 2.0;
                bounds[1] = P + d / 2.0;
            }
 
            
 
            maths::Vector<double> bounds[2];  
            maths::Vector<double> d; 
            maths::Vector<double> diag[3]; 
            maths::Vector<double> calcCoM() { return maths::dzero; }
            void setOctree(bool isOctree) { this->isOctree = isOctree; } 
        };
 
        std::ostream& operator<< (std::ostream& os, Box B); 
    }
}