gnu parallel

DiffractionGratings/grating3D.pro

@@ -436,11 +436,6 @@ PostProcessing {
 PostOperation {
   { Name postop_helmholtz_vector; NameOfPostProcessing postpro_helmholtz_vector ;
     Operation {
-      Print [ uz , OnElementsOf Gama, File StrCat[myDir,"uz.pos"]];
-      Print [ u  , OnPlane { {0.5*(-period_x-dys), -dyc/2,hh_L_6}
-                             {0.5*( period_x-dys), -dyc/2,hh_L_6}
-                             {0.5*(-period_x+dys),  dyc/2,hh_L_6} }
-                              {npts_checkpoyX-1,npts_checkpoyY-1} , File StrCat[myDir,"uz_check.pos"]];
       If (FlagOutEscaFull==1)
         If (Flag_interp_cubic==1)
           Print [ u , OnBox { {-period_x/2,-period_y/2,hh_L_6-PML_bot} {period_x/2,-period_y/2,hh_L_6-PML_bot} {-period_x/2,period_y/2,hh_L_6-PML_bot} {-period_x/2,-period_y/2,hh_L_1+thick_L_1+PML_top} } {npts_interpX,npts_interpY,npts_interpZSca} , File StrCat[myDir,"u_grid.pos"], Name "u_grid"];

DiffractionGratings/grating3D_data_halfellipsoid.geo

@@ -65,11 +65,11 @@ DefineConstant[
     FlagLinkFacets = {0      , Name StrCat[pp5,"/8FlagLinkFacets? [-]"], Choices {0,1}, Visible 0},
         
     InterpSampling     = { 10   , Name StrCat[pp6,"/0Interpolation grid step [nm]"]},
-    Flag_interp_cubic  = { 1    , Name StrCat[pp6,"/1Interpolate on cubic grid?"], Choices {0,1} },
+    Flag_interp_cubic  = { 0    , Name StrCat[pp6,"/1Interpolate on cubic grid?"], Choices {0,1} },
     FlagOutEtotCuts    = { 1    , Name StrCat[pp6,"/2Output Total Electric Field cuts?"] , Choices {0,1} },
     FlagOutHtotCuts    = { 0    , Name StrCat[pp6,"/3Output Total Magnetic Field cuts?"] , Choices {0,1} },
-    FlagOutEscaCuts    = { 1    , Name StrCat[pp6,"/4Output Scattered Electric Field cuts?"] , Choices {0,1} },
-    FlagOutPoyCut      = { 1    , Name StrCat[pp6,"/5Output Poynting 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    = { 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_parallel_Mmatrix.sh

@@ -2,30 +2,41 @@
 export OPENBLAS_NUM_THREADS=1
 export OMP_NUM_THREADS=1
 
-export NPROC=64
-export nb_phi=64
+# 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 myFile="run_phideg_values_"$nb_phi".out"
 
 rm -r $myDir
 rm $myFile
 mkdir $myDir
+gmsh grating3D.geo -3 -o grating3D.msh
 
 myfunc() {
-    local mysubDir=$myDir/run$1_psi$2
+    local mysubDir=$myDir/run_lam$1_phi$2_psi$3
     mkdir $mysubDir
-    cp grating3D.geo grating3D.pro grating3D_data_halfellipsoid.geo grating3D_data.geo grating3D_materials.pro $mysubDir
-    local phi=$(echo "scale=10;360/($nb_phi)*$1" | bc )
-    if [ $2 -eq 0 ]; then echo $1"  "$phi >> $myDir/$myFile ; fi
+    cp grating3D.msh grating3D.pro grating3D_data_halfellipsoid.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
-    echo gmsh grating3D.pro -setstring test_case halfellipsoid -setnumber thetadeg 50 -setnumber phideg $phi -setnumber psideg $2 -setstring Solver.SocketName id$1$2 -
-    if [ $2 -eq  0 ]; then cp res3D/rs.txt ../r_pin_sout_phi$1.out ; fi
-    if [ $2 -eq  0 ]; then cp res3D/rp.txt ../r_pin_pout_phi$1.out ; fi
-    if [ $2 -eq 90 ]; then cp res3D/rs.txt ../r_sin_sout_phi$1.out ; fi
-    if [ $2 -eq 90 ]; then cp res3D/rp.txt ../r_sin_pout_phi$1.out ; fi
+    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
+    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
+    if [ $3 -eq 1 ]; then cp res3D/rp.txt ../r_sin_pout_lam$1_phi$2.out ; fi
+    cp res3D/eff_t1.txt ../eff_t1_lam$1_phi$2_psi$3.out
+    cp res3D/eff_r1.txt ../eff_r1_lam$1_phi$2_psi$3.out
+    cp res3D/eff_t2.txt ../eff_t2_lam$1_phi$2_psi$3.out
+    cp res3D/eff_r2.txt ../eff_r2_lam$1_phi$2_psi$3.out
+    cp res3D/temp-Q_scat.txt ../Q_scat_lam$1_phi$2_psi$3.out
     cd ../..
+    rm -r $mysubDir
 }
 
 export -f myfunc
-# parallel --lb -j $NPROC myfunc ::: $(seq 1 $nb_phi) ::: 0 90
-parallel -j $NPROC myfunc ::: $(seq 1 $nb_phi) ::: 0 90
+parallel -j $NPROC myfunc ::: $(seq 0 $(($nb_lam-1))) ::: $(seq 0 $(($nb_phi-1))) ::: 0 1