5BIN_PATH="/Users/likask/MyBuild/um_debug/solid_shell_prism_element"
6SRC_PATH="/Users/likask/MyBuild/um_debug/users_modules/solid_shell_prism_element"
8MPIRUN="/opt/petsc/arch-darwin-c-debug/bin/mpirun -np $NB_PROC"
15# echo LOG: Mesh 9 - Convergence in plane - locking
17# for ODISP in `seq 1 8`;
19# $MPIRUN $BIN_PATH/solid_shell \
20# -my_file $SRC_PATH/meshes/pinched_cylinder_9.h5m \
21# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
22# -my_order $ODISP -my_order_thickness 0 -my_order_3d 0 \
23# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
25# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 0 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
28# echo LOG: Mesh 9 - Convergence in plane - thinckes 2 - plane locking
30# for ODISP in `seq 1 8`;
32# $MPIRUN $BIN_PATH/solid_shell \
33# -my_file $SRC_PATH/meshes/pinched_cylinder_9.h5m \
34# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
35# -my_order $ODISP -my_order_thickness 2 -my_order_3d 0 \
36# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
38# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 2 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
41# echo LOG: Mesh 9 - Convergence in plane - thinckes 2 - plane 3
43# for ODISP in `seq 1 8`;
45# $MPIRUN $BIN_PATH/solid_shell \
46# -my_file $SRC_PATH/meshes/pinched_cylinder_9.h5m \
47# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
48# -my_order $ODISP -my_order_thickness 3 -my_order_3d 2 \
49# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
51# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 2 2 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
54# echo LOG: Mesh 6 - Convergence in plane - locking
56for ODISP in `seq 5 5`;
58 FILE=$SRC_PATH/meshes/pinched_cylinder_6.h5m
59 ./solid_shell_partition_mesh -my_file $FILE -my_out_file tmp.h5m -my_nb_parts $NB_PROC
60 $MPIRUN $BIN_PATH/solid_shell \
61 -my_file tmp.h5m -my_is_partitioned 1 \
62 -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
63 -my_order $ODISP -my_order_thickness 3 -my_order_3d 2 \
64 -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
66 echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 0 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
69# echo LOG: Mesh 6 - Convergence in plane - thinckes 2 - plane locking
71# for ODISP in `seq 1 8`;
73# $MPIRUN $BIN_PATH/solid_shell \
74# -my_file $SRC_PATH/meshes/pinched_cylinder_6.h5m \
75# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
76# -my_order $ODISP -my_order_thickness 2 -my_order_3d 0 \
77# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
79# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 2 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
82# echo LOG: Mesh 6 - Convergence in plane - thinckes 3 - plane 2
84# for ODISP in `seq 1 8`;
86# $MPIRUN $BIN_PATH/solid_shell \
87# -my_file $SRC_PATH/meshes/pinched_cylinder_6.h5m \
88# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
89# -my_order $ODISP -my_order_thickness 2 -my_order_3d 3 \
90# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
92# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 3 2 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
97# echo LOG: Mesh 4 - Convergence in plane - locking
99# for ODISP in `seq 1 8`;
101# $MPIRUN $BIN_PATH/solid_shell \
102# -my_file $SRC_PATH/meshes/pinched_cylinder_4.h5m \
103# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
104# -my_order $ODISP -my_order_thickness 0 -my_order_3d 0 \
105# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
107# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 0 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
110# echo LOG: Mesh 4 - Convergence in plane - thinckes 2 - plane locking
112# for ODISP in `seq 1 8`;
114# $MPIRUN $BIN_PATH/solid_shell \
115# -my_file $SRC_PATH/meshes/pinched_cylinder_4.h5m \
116# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
117# -my_order $ODISP -my_order_thickness 2 -my_order_3d 0 \
118# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
120# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 2 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
123# echo LOG: Mesh 4 - Convergence in plane - thinckes 3 - plane 2
125# for ODISP in `seq 1 8`;
127# $MPIRUN $BIN_PATH/solid_shell \
128# -my_file $SRC_PATH/meshes/pinched_cylinder_4.h5m \
129# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
130# -my_order $ODISP -my_order_thickness 3 -my_order_3d 2 \
131# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
133# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 3 2 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
138# echo LOG: Mesh 2 - Convergence in plane - locking
140# for ODISP in `seq 1 6`;
142# $MPIRUN $BIN_PATH/solid_shell \
143# -my_file $SRC_PATH/meshes/pinched_cylinder_2.h5m \
144# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
145# -my_order $ODISP -my_order_thickness 0 -my_order_3d 0 \
146# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
148# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 0 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
151# echo LOG: Mesh 2 - Convergence in plane - thinckes 2 - plane locking
153# for ODISP in `seq 1 6`;
155# $MPIRUN $BIN_PATH/solid_shell \
156# -my_file $SRC_PATH/meshes/pinched_cylinder_2.h5m \
157# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
158# -my_order $ODISP -my_order_thickness 2 -my_order_3d 0 \
159# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
161# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 2 0 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
164# echo LOG: Mesh 2 - Convergence in plane - thinckes 3 - plane 2
166# for ODISP in `seq 1 6`;
168# $MPIRUN $BIN_PATH/solid_shell \
169# -my_file $SRC_PATH/meshes/pinched_cylinder_2.h5m \
170# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
171# -my_order $ODISP -my_order_thickness 3 -my_order_3d 2 \
172# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3 \
174# echo LOG: `grep 'Problem SHELL_PROBLEM Nb' tmp_log | awk '{print $5}' | uniq` $ODISP 3 2 `grep "DISPLACEMENT \[ 1 \]" tmp_log | awk '{print $5}'`
179#/opt/local/bin/valgrind --dsymutil=yes --track-origins=yes --trace-children=yes
181# $MPIRUN $BIN_PATH/solid_shell \
182# -my_file $SRC_PATH/meshes/pinched_cylinder_4.h5m \
183# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
184# -my_order 5 -my_order_thickness 2 -my_order_3d 0 \
185# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3e6 \
188# $MPIRUN $BIN_PATH/solid_shell \
189# -my_file $SRC_PATH/meshes/simple_plate.h5m \
190# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
191# -my_order 2 -my_order_thickness 0 -my_order_3d 0 \
192# -my_shell_thickness 3 -my_poisson_ratio 0.3 -my_young_modulus 3e6 \
198# $MPIRUN $BIN_PATH/solid_shell_nonlinear \
199# -my_file $SRC_PATH/meshes/pinched_cylinder_4_long_tension.h5m \
200# -snes_monitor -snes_max_it 24 -snes_atol 1e-6 -snes_rtol 1e-5 -snes_stol 0 \
201# -snes_linesearch_type l2 -snes_linesearch_monitor -snes_linesearch_minlambda 0.01 \
202# -snes_converged_reason \
203# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
204# -ksp_converged_reason \
205# -my_step_size 1 -my_arc_beta 0 \
206# -my_numerical_damping 0 \
207# -my_order 6 -my_order_thickness 2 -my_order_3d 0 \
208# -my_shell_thickness 0.094 -my_poisson_ratio 0.3125 -my_young_modulus 10.5 \
209# -ts_monitor -ts_dt 0.001 -ts_final_time 100000 \
210# -ts_max_snes_failures -1 \
214# $MPIRUN $BIN_PATH/solid_shell_nonlinear \
215# -my_file $SRC_PATH/meshes/pinched_cylinder_6_long_compression.h5m \
216# -snes_monitor -snes_max_it 24 -snes_atol 1e-6 -snes_rtol 1e-5 -snes_stol 0 \
217# -snes_linesearch_type l2 -snes_linesearch_monitor -snes_linesearch_minlambda 0.01 \
218# -snes_converged_reason \
219# -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -ksp_atol 1e-12 -ksp_rtol 1e-12 \
220# -ksp_converged_reason \
221# -my_step_size 1 -my_arc_beta 0 \
222# -my_numerical_damping 0 \
223# -my_order 5 -my_order_thickness 2 -my_order_3d 0 \
224# -my_shell_thickness 1 -my_poisson_ratio 0.3 -my_young_modulus 30 \
225# -ts_monitor -ts_dt 0.001 -ts_final_time 100000 \
226# -ts_max_snes_failures -1 \