oops

DiffractionGratings/grating3D.geo

@@ -291,9 +291,9 @@ For k In {0:8}
 EndFor
 
 
-// If (tag_geom==3) // Split U weird otherwise
-//   Mesh.Algorithm = 6;
-// EndIf
+If (tag_geom==3) // Split U weird otherwise
+  Mesh.Algorithm = 6;
+EndIf
 
 
 // Mesh.SurfaceEdges = 0;

DiffractionGratings/grating3D_data.geo

 
 DefineConstant[
-  skewed_lattice = {"skewed_lattice",
+  test_case = {"skewed_lattice",
     Name "1Geometry",
     Choices {"half_ellipsoid", "hole", "pyramid", "U", "bi_sinusoidal", "retrieve_2D_lamellar", "nanowire_solarcell", "convergence", "skewed_lattice"},
     GmshOption "Reset", Autocheck 0

DiffractionGratings/grating3D_parallel_Mmatrix.sh

@@ -2,18 +2,19 @@
 export OPENBLAS_NUM_THREADS=1
 export OMP_NUM_THREADS=1
 export NPROC=96
-
-export nb_lam=100
-export nb_phi=59
+export nb_lam=1
+export nb_phi=1
+# export nb_lam=100
+# export nb_phi=59
 export lambda_min=350
 export lambda_max=800
 export FLAG_TOTAL=0
-export GRATING_CASE="U"
+export GRATING_CASE="half_ellipsoid"
 export myDir="res_Matrix_nb_lam"$nb_lam"_nb_phi"$nb_phi"_total"$FLAG_TOTAL
 
 rm -r $myDir
 mkdir $myDir
-gmsh grating3D.geo -3 -o grating3D.msh
+gmsh grating3D.geo -setstring test_case $GRATING_CASE -3 -o grating3D.msh
 
 myfunc() {
     local mysubDir=$myDir/run_lam$1_phi$2_psi$3