clean

DiffractionGratings/grating3D.geo

@@ -22,30 +22,30 @@ e = 5*nm;
 E = 10*(period_x+period_x/2);
 
 Macro SetPBCs
-  masterX() = Surface In BoundingBox{.5*(-period_x-dys)-e,-dyc/2-e, PML_bot_hh-e,(-period_x+dys)/2+e, dyc/2+e, PML_top_hh+PML_top+e};
-  slaveX()  = Surface In BoundingBox{.5*( period_x-dys)-e,-dyc/2-e, PML_bot_hh-e,( period_x+dys)/2+e, dyc/2+e, PML_top_hh+PML_top+e};
-  masterY() = Surface In BoundingBox{.5*(-period_x-dys)-e,-dyc/2-e, PML_bot_hh-e,( period_x+dys)/2+e,-dyc/2+e, PML_top_hh+PML_top+e};
-  slaveY()  = Surface In BoundingBox{.5*(-period_x-dys)-e, dyc/2-e, PML_bot_hh-e,( period_x+dys)/2+e, dyc/2+e, PML_top_hh+PML_top+e};
-  // For k In {1:#masterX()-1}
-  //   Printf("masterX surf %g",masterX(k));
+  BlochXm() = Surface In BoundingBox{.5*(-period_x-dys)-e,-dyc/2-e, PML_bot_hh-e,(-period_x+dys)/2+e, dyc/2+e, PML_top_hh+PML_top+e};
+  BlochXp() = Surface In BoundingBox{.5*( period_x-dys)-e,-dyc/2-e, PML_bot_hh-e,( period_x+dys)/2+e, dyc/2+e, PML_top_hh+PML_top+e};
+  BlochYm() = Surface In BoundingBox{.5*(-period_x-dys)-e,-dyc/2-e, PML_bot_hh-e,( period_x+dys)/2+e,-dyc/2+e, PML_top_hh+PML_top+e};
+  BlochYp() = Surface In BoundingBox{.5*(-period_x-dys)-e, dyc/2-e, PML_bot_hh-e,( period_x+dys)/2+e, dyc/2+e, PML_top_hh+PML_top+e};
+  // For k In {1:#BlochXm()-1}
+  //   Printf("BlochXm surf %g",BlochXm(k));
   // EndFor
-  // For k In {1:#masterY()-1}
-  //   Printf("masterY surf %g",masterY(k));
+  // For k In {1:#BlochYm()-1}
+  //   Printf("BlochYm surf %g",BlochYm(k));
   // EndFor
-  // For k In {1:#slaveX()-1}
-  //   Printf("slaveX surf %g",slaveX(k));
+  // For k In {1:#BlochXp()-1}
+  //   Printf("BlochXp surf %g",BlochXp(k));
   // EndFor
-  // For k In {1:#slaveY()-1}
-  //   Printf("slaveY surf %g",slaveY(k));
+  // For k In {1:#BlochYp()-1}
+  //   Printf("BlochYp surf %g",BlochYp(k));
   // EndFor
   If (tag_geom==6) // bi-sin : BoundingBox does not catch BSpline Surfaces?
-    masterX()+={17,12};
-    slaveX()+={19,14};
-    masterY()+={15,20};
-    slaveY()+={13,18};
+    BlochXm()+={17,12};
+    BlochXp()+={19,14};
+    BlochYm()+={15,20};
+    BlochYp()+={13,18};
   EndIf
-  Periodic Surface{slaveX()} = {masterX()} Translate{period_x,        0,        0};
-  Periodic Surface{slaveY()} = {masterY()} Translate{    dys, dyc,        0};
+  Periodic Surface{BlochXp()} = {BlochXm()} Translate{period_x,        0,        0};
+  Periodic Surface{BlochYp()} = {BlochYm()} Translate{    dys, dyc,        0};
 Return
 
 If (tag_geom==6)
@@ -230,10 +230,10 @@ Else
   Physical Volume("PMLTOP"        ,8) = {8};
   Physical Volume("SCAT"          ,9) = {9};
 EndIf
-Physical Surface("BXM" ,101 ) = masterX();
-Physical Surface("BXP" ,102 ) = slaveX();
-Physical Surface("BYM" ,201 ) = masterY();
-Physical Surface("BYP" ,202 ) = slaveY();
+Physical Surface("BXM" ,101 ) = BlochXm();
+Physical Surface("BXP" ,102 ) = BlochXp();
+Physical Surface("BYM" ,201 ) = BlochYm();
+Physical Surface("BYP" ,202 ) = BlochYp();
 Physical Surface("STOP",301 ) = Surface In BoundingBox{-period_x/2-E,-period_y/2-E, PML_top_hh-e,period_x/2+E, period_y/2+E, PML_top_hh+e};
 Physical Surface("SBOT",302 ) = Surface In BoundingBox{-period_x/2-E,-period_y/2-E, hh_L_6-e,period_x/2+E    , period_y/2+E, hh_L_6+e};
 Physical Surface("SPMLTOP",401 ) = Surface In BoundingBox{-period_x/2-E,-period_y/2-E, PML_top_hh+PML_top-e,period_x/2+E, period_y/2+E, PML_top_hh+PML_top+e};

DiffractionGratings/grating3D.pro

@@ -346,21 +346,10 @@ PostProcessing {
   { Name postpro_helmholtz_vector; NameOfFormulation helmholtz_vector; NameOfSystem M;
     Quantity {
       { Name u      ; Value { Local { [ {u}       ]; In Omega; Jacobian JVol; } } }      
-      { Name uper   ; Value { Local { [       {u} *Exp[-I[]*(skx1[]*X[]+sky1[]*Y[])]  ]; In Omega; Jacobian JVol; } } }
-      { Name uperx  ; Value { Local { [ CompX[{u} *Exp[-I[]*(skx1[]*X[]+sky1[]*Y[])]] ]; In Omega; Jacobian JVol; } } }
-      { Name upery  ; Value { Local { [ CompY[{u} *Exp[-I[]*(skx1[]*X[]+sky1[]*Y[])]] ]; In Omega; Jacobian JVol; } } }
-      { Name uperz  ; Value { Local { [ CompZ[{u} *Exp[-I[]*(skx1[]*X[]+sky1[]*Y[])]] ]; In Omega; Jacobian JVol; } } }
-      { Name E1perx ; Value { Local { [ CompX[E1[]*Exp[-I[]*(skx1[]*X[]+sky1[]*Y[])]] ]; In Omega; Jacobian JVol; } } }
-      { Name E1pery ; Value { Local { [ CompY[E1[]*Exp[-I[]*(skx1[]*X[]+sky1[]*Y[])]] ]; In Omega; Jacobian JVol; } } }
-      
-      { Name Etotpery ; Value { Local { [ CompY[({u}+E1[])*Exp[-I[]*(k1x[]*X[]+k1y[]*Y[])]] ]; In Omega; Jacobian JVol; } } }
-      
       { Name Etot   ; Value { Local { [ {u}+E1[]  ]; In Omega; Jacobian JVol; } } }
       { Name Htot   ; Value { Local { [ H1[]-I[]/(mur[]*mu0*om0)*{Curl u}]; In Omega; Jacobian JVol; } } }
-      { Name Htotx  ; Value { Local { [CompX[ H1[]-I[]/(mur[]*mu0*om0)*{Curl u}]]; In Omega; Jacobian JVol; } } }
       { Name Edif   ; Value { Local { [ {u}+E1d[] ]; In Omega; Jacobian JVol; } } }
       { Name E1     ; Value { Local { [     E1[]  ]; In Omega; Jacobian JVol; } } }
-      { Name H1x     ; Value { Local { [CompX[H1[]]  ]; In Omega; Jacobian JVol; } } }
       { Name Em     ; Value { Local { [     E1d[] ]; In Omega; Jacobian JVol; } } }
       { Name source ; Value { Local { [  source[] ]; In Omega; Jacobian JVol; } } }
       { Name epsr_xx; Value { Local { [  CompXX[epsr[]] ]; In Omega; Jacobian JVol; } } }
@@ -393,13 +382,6 @@ PostProcessing {
           { Name numbering_ij~{i}~{j}   ; Value { Term{Type Global; [Vector[i-Nmax,j-Nmax,0]] ; In SurfIntBot ; } } }
         EndFor
       EndFor
-
-      // For i In {0:Nb_ordre-1}
-      //   For j In {0:Nb_ordre-1}
-      //     { Name alpha~{i}~{j} ; Value { Term{Type Global; [alpha~{i}[]-alpha~{i}[]] ; In SurfIntBot ; } } }
-      //     { Name beta~{j}~{j}  ; Value { Term{Type Global; [beta~{j}[]] ; In SurfIntBot ; } } }
-      //   EndFor
-      // EndFor
     }
   }
 }
@@ -485,7 +467,6 @@ PostOperation {
                                   {0.5*( period_x-dys), -dyc/2, hh_L_1+thick_L_1/2}
                                   {0.5*(-period_x+dys),  dyc/2, hh_L_1+thick_L_1/2} }
                                   {npts_checkpoyX-1,npts_checkpoyY-1} , File StrCat[myDir,"Poy_inc_gd.pos"], Format Table];
-      Print[ Abs_scat2[Scat]  , OnGlobal, File > StrCat[myDir,"temp-Q_scat2.txt"], Format Table ];
 
       For k In {2:6}
         Print[ Abs_L~{k}[L~{k}], OnGlobal, File > StrCat[myDir,Sprintf["temp-Q_L_%g.txt",k]], Format Table ];