diff --git a/DiffractionGratings/grating3D_parallel_Mmatrix.sh b/DiffractionGratings/grating3D_parallel_Mmatrix.sh index 9fc960ccb5b8864fb6d6611cf23916e2619ab62e..73dac5e36cacfab085e06829bda9966884d027eb 100644 --- a/DiffractionGratings/grating3D_parallel_Mmatrix.sh +++ b/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