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

Go to the documentation of this file.

#pragma once

#ifndef GRATINGLAYER_H
#define GRATINGLAYER_H

#include <iostream>
#include <math.h>
#include <complex>
#include <string>

#include "tbci/std_cplx.h"
#include "tbci/vector.h"
#include "tbci/matrix.h"
#include "tbci/f_matrix.h"
#include "tbci/lapack_stdcplx.h"

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






class GratingLayer1D{
private:                
        complex <double>  I;
        complex<double> n0;
        complex<double> gamma0;
        // thickness of the layer, period of structure, wavelength
        double thickness, period, lambda, theta;
        // Indicates to which order the fourier series of the field and the permitivity is truncated
        int nOrds;
        // the minimum/maximum Order of the fourier series = (nOrds-1)/2
        int minOrd, maxOrd;
        // number of Segments per layer
        int nPermSegs;
        // refractive index of each segments
        Vector<double> transPt;
        Vector < complex <double> >  refInd;
        F_Matrix< complex<double> > epsMat, invEpsMat, epsMat_inv, invEpsMat_inv, B;
        // the gamma values of the field
        Vector < complex <double>  > gamma;
        // the permitivity coefficients of the fourier series
        Vector < complex <double> > permCoeff;
        Vector < complex <double> > permCoeffInv;
        // The eigensystem in row-major order
        F_Matrix < complex <double> > eigMat;
        // The eigenvectors: m_th column contains the m_th eigenvector
        // Eigenvalues
        Vector < complex<double> > eigVal;
        // alpha[i] = sqrt(EigVal[i])
        Vector <complex <double> > alpha;
        // polarization mode of the incident wave (TE/TM)
        string polMode;
        int minPermCoeff;
        int maxPermCoeff;
        int nPermCoeffs;
                
public:
        F_Matrix < complex<double> > eigVecE;
        F_Matrix < complex<double> > eigVecH;
        F_Matrix < complex<double> > invEigVecE;
        F_Matrix < complex<double> > invEigVecH;
        
        

        GratingLayer(int,int);

        ~GratingLayer();

        void setParameters(double,double,complex<double>,double,double,Vector<double>,Vector< complex<double> >,string);

        void solve(string);

        void writeIndex(int, string);




        // set and get data members
        void setPolMode(string polarization);   
        string getPolMode(void);
        void setThickness(double);              
        double getThickness(void);
        // incident wavelength in vacuum
        void setLambda(double);                 double getLambda(void);
        void setN0(complex<double>);            
        complex<double> getN0(void);
        void setPeriod(double);         
        double getPeriod(void);
        void setnPermSegs(int);         
        int getnPermSegs(void);
        void setTransPt(Vector <double>);       
        double getTransPt(int);
        void setAngle(double);          
        double getAngle(void);
        void setRefInd(Vector< complex<double> >);
        complex<double> getRefInd(int);

        int getnOrds(void);
        complex<double> getPermCoeff(int);
        complex<double> getAlpha(int);
        complex<double> getEigVal(int);
        complex<double> getEigMat(int,int);
        complex<double> getGamma(int);
        // 1: vector
        // 2: element
        complex<double> getEigVecE(int vec, int comp);
        complex<double> getEigVecH(int vec, int comp);
        complex<double> getInvEigVecE(int vec, int comp);
        complex<double> getInvEigVecH(int vec, int comp);



        void CalculatePermitivity(void);
        void createEigMat(void);
        void solveEigMat(void);
        void calcSMatrix(void);
        inline F_Matrix < complex<double> > unit(int i)
        {
                F_Matrix< complex<double> > m(i);
                m.setunit();
                return m;
        }

};

#endif

---