/dbfs1/fh22/linux/src/FMMLIB/Grating1D.h

Go to the documentation of this file.

#pragma once

#ifndef GRATING_H
#define GRATING_H


#include "GratingLayer1D.h"
#include <vector>

#include "tbci/f_bandmatrix.h"
#include "tbci/cscmatrix.h"

using std::complex;
using std::string;
using std::vector;
using TBCI::F_Matrix;
using TBCI::F_BandMatrix;
using TBCI::Vector;
using TBCI::CSCMatrix;


class Grating1D{

private:
        double theta;
        double lambda;
        double period;
        double gratingSize; 
        double padding;
        Vector< complex<double> > r;
        Vector< complex<double> > t;
        Vector< complex<double> > refl;
        Vector< complex<double> > tran, transS;
        Vector <double> rEFF;
        Vector <double> tEFF;
        complex <double> n0;
        complex <double> n3;
        complex <double> nPad;
        //CSCMatrix< complex<double> > lhsMatBand, precond;
        F_BandMatrix< complex<double> > precond;
        F_Matrix< complex<double> >invA;
        F_BandMatrix< complex<double> > lhsMatBand;
        Vector < complex<double> > rhsVec;
        Vector < complex<double> > solVec;
        Vector < complex<double> > incid;
        Vector < complex<double> > gamma;  
        complex <double> I;
        complex <double> za, zb;
        int nLayers;
        int nOrds
        int minOrd
        int maxOrd;
        int nUnknowns;
        int nPermSegs, nPadSegs, gratingSegs;
        string polMode;
        string initMode;
        F_Matrix<int> material, MCoeffs;
        int nMat;
        int matPad;
        F_Matrix < complex<double> > indices;
        Vector< complex<double> > refIndMat;
        F_Matrix <double> transpt;
        Vector<double> slabThick;
        F_Matrix<double> fieldAmplitude;
        int dimX;
        int dimY;
        // for calculation of external field of embedded structure
        double embedding, embPeriod, propAfter;
        int nOrdsEmb;
        Vector < complex<double> > TransEmb, tEmb, gammaEmb;
        double plotSize;
        // for the scattered matrix approach
        F_Matrix< complex<double> > Q1, Q2, G, tau, R, T, R_tilde, T_tilde, F;
        Vector < complex<double> > expTerm;
        double incidPowFactor;
public:
        //public variables
        vector<GratingLayer> layer;
        
        //public functions

        Grating(int truncation, int layers, int nSegments);
        ~Grating();
        void setParameters(double,double,double,complex<double>,complex<double>,complex<double>,double,Vector< complex<double> >,string);
        void initLayers(string,string,string);
        void initLayers(F_Matrix<int>,Vector<double>,string);
        void initLayers(F_Matrix< complex<double> >,Vector<double>);
        void solve(string);
        void calcIntField(int,int);
        void calcExtField(int,int,double,double,double);
        void calcEmbField(int,int,double,double,double);
        complex<double> getIntField(double,double);
        complex<double> getExtField(double,double);
        complex<double> getEmbExtField(double,double);
        void writeField(string);
        void calcTransEmb(double embSize, int embOrd);
        double getEfficiencyR(int);
        double getEfficiencyT(int);
        complex<double> getTrans(int);
        complex<double> getRef(int);
        
        void calcIncid(Vector < complex <double> >,int);
        
        
        void changeRefInd(int,int,complex<double>);
        void changeMaterial(int slab, int segment, int mat);
        void addFrame(complex<double>,double);
        void writeParameters(string);




        void setTheta(double);
        double getTheta();
        void setLambda(double); 
        double getLambda();
        void setPeriod(double); 
        double getPeriod(void);
        double getGratingSize(void);
        void setnPermSegs(int); 
        int getnPermSegs(void);
        void setPolMode(string);
        string getPolMode();
        void setN0(complex<double>);
        complex<double> getN0();
        void setN3(complex<double>);
        complex<double> getN3();
        void setNPad(complex<double>);
        complex<double> getNPad();
        void setIncid(Vector< complex<double> >);
        complex<double>getIncid(int);
        void setMatPad(int);
        void setPadding(double);        
        double getPadding(void);
        int getnOrds();
        int getnLayers(void);
        double getEmbSize();
        complex<double> getLhsMat(int,int);
        complex<double>getRhsVec(int);
        complex<double> getA(int,int);
        complex<double> getAdash(int,int);
        complex<double>getGamma(int);   
        
        


        void setEquations(void);
        void setLhsMatTE(void);
        void setLhsMatTM(void);
        void  solveEquations(void);
        void setSMatrix(void);
        void findEfficiency(string);
        void findEfficiencyT(void);
        void findEfficiencyS(void);
        void calcSTR(void);
        void importMaterials(string);
        void calcGamma(void);
        void setLhsMatBand(int,int,complex<double>);
        // find the number of lines of a file
        // 1: filename
        int nLines(string);
        // return 1 if both parameters are equal, zero else
        double delta(int,int);

        inline F_Matrix < complex<double> > unit(int i)
        {
                F_Matrix< complex<double> > m(i);
                m.setunit();
                return m;
        }

};

#endif

---