aim at gmsh/getdp directly

DiffractionGratings/grating3D_parallel_Mmatrix.sh

@@ -6,21 +6,21 @@ export OMP_NUM_THREADS=1
 export NPROC=96
 
 ### set gmsh/getdp path if necessary
-# export GMSHPATH='/path/to/onelabdir'
-# export GETDPPATH='/path/to/onelabdir'
-export GMSHPATH=''
-export GETDPPATH=''
+# export mygmsh='/path/to/onelabdir/gmsh'
+# export mygetdp='/path/to/onelabdir/getdp'
+export mygmsh='gmsh'
+export mygetdp='getdp'
 
 ### choose configuration
-### theta loop study for the "U" case
-### config 1
+## theta loop study for the "U" case
+## config 1
 # export flag_angle_study="theta"
 # export loop_angle_max=50
 # export fixed_angle=90
 # export GRATING_CASE="U"
 
-### phi loop study for the "U" case
-### config 2
+## phi loop study for the "U" case
+## config 2
 export flag_angle_study="phi"
 export loop_angle_max=360
 export fixed_angle=50
@@ -35,7 +35,7 @@ export myDir="res_Matrix_nb_lam"$nb_lam"_nb_"$flag_angle_study$nb_angle"_total"$
 
 rm -r $myDir
 mkdir $myDir
-$GMSHPATH/gmsh grating3D.geo -setstring test_case $GRATING_CASE -3 -o grating3D.msh
+$mygmsh grating3D.geo -setstring test_case $GRATING_CASE -3 -o grating3D.msh
 
 myfunc() {
     local mysubDir=$myDir/run_lam$1_$flag_angle_study$2_psi$3
@@ -51,7 +51,7 @@ myfunc() {
     fi
     local psi=$(echo "scale=10;90*$3" | bc )
     cd $mysubDir
-    $GETDPPATH/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 $theta -setnumber phideg $phi -setnumber psideg $psi -setnumber FLAG_TOTAL $FLAG_TOTAL
+    $mygetdp 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 $theta -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_$flag_angle_study$2.out ; fi
     if [ $3 -eq 0 ]; then cp res3D/rp.txt ../r_pin_pout_lam$1_$flag_angle_study$2.out ; fi
     if [ $3 -eq 1 ]; then cp res3D/rs.txt ../r_sin_sout_lam$1_$flag_angle_study$2.out ; fi