scaling tests

DiffractionGratings/grating3D.geo

@@ -136,7 +136,7 @@ If (tag_geom==2)
 EndIf
 
 If (tag_geom==3)
-  Torus(9) = {0,0,hh_L_3+ry-1,rx,ry};
+  Torus(9) = {0,0,hh_L_3+ry-1*nm,rx,ry};
   Dilate { { 0,0,hh_L_3 }, { 1, 1, rz/ry } } { Volume{9}; }
   BooleanDifference{ Volume{9}; Delete; }{ Volume{4}; }
 EndIf

DiffractionGratings/grating3D_data_2Dlamellar.geo

-nm  = 10000;
+nm  = 1000;
 pp1 = "1Incident Plane Wave";
 pp2 = "2Layers Thicknesses";
 pp3 = "3Scatterer Properties";
@@ -45,12 +45,6 @@ DefineConstant[
     eps_im_L_5    = {2.9524       , Name StrCat[pp4,"/layer 5: imag part of relative permittivity"]},
     eps_re_L_6    = {-44.9757 , Name StrCat[pp4,"/layer 6: real part of relative permittivity"]},
     eps_im_L_6    = {2.9524   , Name StrCat[pp4,"/layer 6: imag part of relative permittivity"]},
-    // eps_re_L_4    = {1  , Name StrCat[pp4,"/layer 4: real part of relative permittivity"]},
-    // eps_im_L_4    = {0  , Name StrCat[pp4,"/layer 4: imag part of relative permittivity"]},
-    // eps_re_L_5    = {1  , Name StrCat[pp4,"/layer 5: real part of relative permittivity"]},
-    // eps_im_L_5    = {0  , Name StrCat[pp4,"/layer 5: imag part of relative permittivity"]},
-    // eps_re_L_6    = {-44.9757 , Name StrCat[pp4,"/layer 6: real part of relative permittivity"]},
-    // eps_im_L_6    = {2.9524   , Name StrCat[pp4,"/layer 6: imag part of relative permittivity"]},
      
     og            = {0          , Name StrCat[pp5,"/0geometrical order [-]"]  , Choices {0="1",1="2"} },
     oi            = {1          , Name StrCat[pp5,"/0interpolation order [-]"], Choices {0="1",1="2"} },
@@ -68,11 +62,11 @@ DefineConstant[
     
     InterpSampling     = { 30   , Name StrCat[pp6,"/0Interpolation grid step [nm]"]},
     Flag_interp_cubic  = { 0    , Name StrCat[pp6,"/1Interpolate on cubic grid?"], Choices {0,1} },
-    FlagOutEtotCuts    = { 0    , Name StrCat[pp6,"/2Output Total Electric Field cuts?"] , Choices {0,1} },
+    FlagOutEtotCuts    = { 1    , Name StrCat[pp6,"/2Output Total Electric Field cuts?"] , Choices {0,1} },
     FlagOutHtotCuts    = { 1    , Name StrCat[pp6,"/3Output Total Magnetic Field cuts?"] , Choices {0,1} },
     FlagOutEscaCuts    = { 0    , Name StrCat[pp6,"/4Output Scattered Electric Field cuts?"] , Choices {0,1} },
     FlagOutPoyCut      = { 0    , Name StrCat[pp6,"/5Output Poynting cuts?"] , Choices {0,1} },
-    FlagOutEtotFull    = { 1    , Name StrCat[pp6,"/6Total Electric Field Full Output?"] , Choices {0,1} },
+    FlagOutEtotFull    = { 0    , Name StrCat[pp6,"/6Total Electric Field Full Output?"] , Choices {0,1} },
     FlagOutEscaFull    = { 0    , Name StrCat[pp6,"/7Scattered Electric Field Full Output?"] , Choices {0,1} },
     FlagOutPoyFull     = { 0    , Name StrCat[pp6,"/8Poynting Full Output?"] , Choices {0,1} }
 ];

DiffractionGratings/grating3D_data_bisin.geo

-nm  = 10000;
+nm  = 1000;
 pp1 = "1Incident Plane Wave";
 pp2 = "2Layers Thicknesses";
 pp3 = "3Scatterer Properties";
@@ -48,9 +48,9 @@ DefineConstant[
     
     og            = {0          , Name StrCat[pp5,"/0geometrical order [-]"]  , Choices {0="1",1="2"} },
     oi            = {1          , Name StrCat[pp5,"/0interpolation order [-]"], Choices {0="1",1="2"} },
-    paramaille    = {7     , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
+    paramaille    = {8     , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
     lc_scat       = {lambda0/(2*paramaille)    , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
-    PML_top       = {lambda0, Name StrCat[pp5,"/4PML top thickness [nm]"]},
+    PML_top       = {lambda0*1.2, Name StrCat[pp5,"/4PML top thickness [nm]"]},
     PML_bot       = {lambda0*1.2, Name StrCat[pp5,"/5PML bot thickness [nm]"]},
     Nmax          = {2      , Name StrCat[pp5,"/6Number of non specular order to output [-]"]},
     refine_mesh_L_1= {1      , Name StrCat[pp5,"/7refine layers/1refine mesh layer 1 [-]"]},

DiffractionGratings/grating3D_data_halfellipsoid.geo

-nm  = 10000;
+nm  = 1000;
 pp1 = "1Incident Plane Wave";
 pp2 = "2Layers Thicknesses";
 pp3 = "3Scatterer Properties";

DiffractionGratings/grating3D_data_hole.geo

-nm  = 10000;
+nm  = 1000;
 pp1 = "1Incident Plane Wave";
 pp2 = "2Layers Thicknesses";
 pp3 = "3Scatterer Properties";

DiffractionGratings/grating3D_data_pyramid.geo

-nm  = 10000;
+nm  = 1000;
 pp1 = "1Incident Plane Wave";
 pp2 = "2Layers Thicknesses";
 pp3 = "3Scatterer Properties";

DiffractionGratings/grating3D_data_solarcell.geo

-nm  = 10000;
+nm  = 1000;
 pp1 = "1Incident Plane Wave";
 pp2 = "2Layers Thicknesses";
 pp3 = "3Scatterer Properties";
@@ -50,7 +50,7 @@ DefineConstant[
     
     og            = {0          , Name StrCat[pp5,"/0geometrical order [-]"]  , Choices {0="1",1="2"} },
     oi            = {1          , Name StrCat[pp5,"/0interpolation order [-]"], Choices {0="1",1="2"} },
-    paramaille    = {7       , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
+    paramaille    = {8       , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
     lc_scat       = {lambda0/(4*paramaille)    , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
     PML_top       = {lambda0, Name StrCat[pp5,"/4PML top thickness [nm]"]},
     PML_bot       = {lambda0, Name StrCat[pp5,"/5PML bot thickness [nm]"]},
@@ -90,23 +90,6 @@ PML_top       = nm * PML_top;
 PML_bot       = nm * PML_bot;
 InterpSampling= nm * InterpSampling;
 
-lambda0       = nm * lambda0;
-period_x      = nm * period_x;
-period_y      = nm * period_y;
-thick_L_1     = nm * thick_L_1;
-thick_L_2     = nm * thick_L_2;
-thick_L_3     = nm * thick_L_3;
-thick_L_4     = nm * thick_L_4;
-thick_L_5     = nm * thick_L_5;
-thick_L_6     = nm * thick_L_6;
-rx            = nm * rx;
-ry            = nm * ry;
-rz            = nm * rz;
-lc_scat       = nm * lc_scat;
-PML_top       = nm * PML_top;
-PML_bot       = nm * PML_bot;
-InterpSampling= nm * InterpSampling;
-
 lambda_m = lambda0;
 
 og+=1;

DiffractionGratings/grating3D_data_torus.geo

-nm  = 10000;
+nm  = 1000;
 pp1 = "1Incident Plane Wave";
 pp2 = "2Layers Thicknesses";
 pp3 = "3Scatterer Properties";