asphericLens.h

Go to the documentation of this file.

#pragma once
#include "objectshape.h"
namespace GOAT
{
    namespace raytracing
    {
        typedef struct
        {
            double k; 
            double R; 
            std::vector<double> A;
            double shift=0; 
            bool isPlano = false;
        } asphericLensSide;
 

        typedef struct
        {
            asphericLensSide left,right;
            double offset;
            double radius;
        } asphericLensParms;
        
   

 
        class asphericLens :
            public ObjectShape
        {
         public:
             asphericLens(const maths::Vector<double>& P,
                 std::complex<double>  n, asphericLensParms lensParms, 
                 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 binWrite();
            void binRead();
            void scale(double sf);
            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>& p);
            double volume();
            void binWrite(std::ofstream& os) { } 
            void binRead(std::ifstream& os) { } 
            void initQuad() { } 
            void setr0(double r0) { } 
            void setPos(maths::Vector<double> r) { } 
            void setPos(double x, double y, double z) { } 
            maths::Vector<double> calcCoM() {
                return maths::Vector<double> (0, 0, 0);
            } 
 
         
            enum Side {left,right,lateral};
            asphericLensParms lensParms;
            double zleft(maths::Vector<double> P, double lambda);
            double zright(maths::Vector<double> P, double lambda);
            double z(double r, asphericLensSide side);
 
            double dzdrleft(maths::Vector<double> P, double lambda);
            double dzdrright(maths::Vector<double> P, double lambda);
         
/*
            double dzright(double lambda);
            double dzleft(double lambda);
            */
            double drdlambda;
            maths::Vector<double> P, k;
           maths::Vector<double> currentnorm;
        };
    }
}