add solarcell test case

DiffractionGratings/grating3D.pro

@@ -72,7 +72,7 @@ Function{
         epsr2[] = Complex[eps_re_L~{i} , eps_im_L~{i}];
       EndIf
     Else
-      For j In {2:nb_available_materials}
+      For j In {1:nb_available_materials}
         If(flag_mat~{i}==j)
           epsr[L~{i}]  = Complex[epsr_re_interp_mat~{j}[lambda0/nm*1e-9] , epsr_im_interp_mat~{j}[lambda0/nm*1e-9]] * TensorDiag[1,1,1];
           If (i==1)
@@ -151,7 +151,7 @@ Function{
   AmplEis[] =      spol[];
   AmplErs[] = rs[]*spol[];
   AmplEts[] = ts[]*spol[];
-  AmplHis[] = Sqrt[ep0*epsr1[]/mu0]      *spol[];
+  AmplHis[] = Sqrt[ep0*epsr1[]/mu0]     *spol[];
   AmplHrs[] = Sqrt[ep0*epsr1[]/mu0]*rp[]*spol[];
   AmplHts[] = Sqrt[ep0*epsr1[]/mu0]*tp[]*spol[];
 
@@ -337,9 +337,9 @@ PostProcessing {
       { Name Poy_tot; Value { Local { [ 0.5*Re[Cross[{u}+E1[] , Conj[ H1[]-I[]/(mur[]*mu0*om0)*{Curl u}]]] ]; In Omega; Jacobian JVol; } } }
       
       For k In {2:6}
-        { Name Abs_L~{k} ; Value { Integral { [ ep0*om0 * 0.5*Fabs[eps_im_L~{k}]*(SquNorm[{u}+E1[]]) / (Pinc*period_x*period_y) ] ; In L~{k} ; Integration I1 ; Jacobian JVol ; } } }
+        { Name Abs_L~{k} ; Value { Integral { [ ep0*om0 * 0.5*Im[CompXX[epsr[]]]*(SquNorm[{u}+E1[]]) / (Pinc*period_x*period_y) ] ; In L~{k} ; Integration I1 ; Jacobian JVol ; } } }
       EndFor
-      { Name Abs_scat ; Value { Integral { [ ep0*om0 * 0.5*Fabs[eps_im_Scat]*(SquNorm[{u}+E1[]]) / (Pinc*period_x*period_y) ] ; In Scat ; Integration I1 ; Jacobian JVol ; } } }
+      { Name Abs_scat ; Value { Integral { [ ep0*om0 * 0.5*Im[CompXX[epsr[]]]*(SquNorm[{u}+E1[]]) / (Pinc*period_x*period_y) ] ; In Scat ; Integration I1 ; Jacobian JVol ; } } }
 
       For i In {0:Nb_ordre-1}
         For j In {0:Nb_ordre-1}
@@ -410,6 +410,7 @@ PostOperation {
         Print [ Poy_tot , OnPlane { {0,-period_y/2,hh_L_6} {0,period_y/2,hh_L_6} {0,-period_y/2,hh_L_1+thick_L_1} } {npts_interpY,npts_interpZTot} , File StrCat[myDir,"Poy_tot_cut_X=0.pos"], Name "Poy_tot_cut_X=0"];
       EndIf
       // // // // For DEBUG
+      // Print [ epsr_xx , OnElementsOf Omega, File StrCat[myDir,"epsr_xx.pos"]];
       // Print [ H1x , OnElementsOf Omega, File StrCat[myDir,"H1x.pos"]];
       // Print [ Etot , OnElementsOf Omega, File StrCat[myDir,"Etot.pos"]];
       // Print [ Htotx , OnElementsOf Omega, File StrCat[myDir,"Htotx.pos"]];