more gnu parallel

DiffractionGratings/grating3D_parallel_Mmatrix.sh

@@ -2,13 +2,12 @@
 export OPENBLAS_NUM_THREADS=1
 export OMP_NUM_THREADS=1
 
-# export NPROC=18
-# export nb_phi=3
-# export nb_lam=3
-export NPROC=100
-export nb_phi=50
-export nb_lam=100
-export myDir="res_Matrix_nb_phi"$nb_phi
+export NPROC=96
+export lambda_min=400
+export lambda_max=800
+export FLAG_TOTAL=0
+export GRATING_CASE="halfellipsoid"
+export myDir="res_Matrix_nb_lam"$nb_lam"_nb_phi"$nb_phi"_total"$FLAG_TOTAL
 export myFile="run_phideg_values_"$nb_phi".out"
 
 rm -r $myDir
@@ -19,12 +18,12 @@ gmsh grating3D.geo -3 -o grating3D.msh
 myfunc() {
     local mysubDir=$myDir/run_lam$1_phi$2_psi$3
     mkdir $mysubDir
-    cp grating3D.msh grating3D.pro grating3D_data_halfellipsoid.geo grating3D_data.geo grating3D_materials.pro $mysubDir
+    cp grating3D.msh grating3D.pro grating3D_data_$GRATING_CASE.geo grating3D_data.geo grating3D_materials.pro $mysubDir
     local lam=$(echo "scale=10;400+400/($nb_lam-1)*$1" | bc )
     local phi=$(echo "scale=10;360/($nb_phi)*$2" | bc )
     local psi=$(echo "scale=10;90*$3" | bc )
     cd $mysubDir
-    getdp grating3D.pro -pre helmholtz_vector -msh grating3D.msh -cal -pos postop_helmholtz_vector -petsc_prealloc 200 -setstring test_case halfellipsoid -setnumber lambda0 $lam -setnumber thetadeg 50 -setnumber phideg $phi -setnumber psideg $psi
+    getdp grating3D.pro -pre helmholtz_vector -msh grating3D.msh -cal -pos postop_helmholtz_vector -petsc_prealloc 200 -setstring test_case $GRATING_CASE -setnumber lambda0 $lam -setnumber thetadeg 50 -setnumber phideg $phi -setnumber psideg $psi -setnumber FLAG_TOTAL $FLAG_TOTAL
     if [ $3 -eq 0 ]; then cp res3D/rs.txt ../r_pin_sout_lam$1_phi$2.out ; fi
     if [ $3 -eq 0 ]; then cp res3D/rp.txt ../r_pin_pout_lam$1_phi$2.out ; fi
     if [ $3 -eq 1 ]; then cp res3D/rs.txt ../r_sin_sout_lam$1_phi$2.out ; fi