v0.10.0
test_bash_codes.sh
Go to the documentation of this file.
1 # # #
2 # # #
3 # # # ODISP="4"
4 # # #
5 # # # ~/build_moFEM1/usr/tools/mofem_part -my_file impinging_cylinder_$ODISP.cub -my_nparts 8
6 # # #
7 # # # echo "Start numerical calculation field split is $FIELD_SPLIT !!! with $NB_PROC processes, mesh h $ODISP and order p 2 ..."
8 # # # mpirun -np 8 ./fe_field_split \
9 # # # -my_file out.h5m \
10 # # # -my_is_partitioned true \
11 # # # -wave_number 10 \
12 # # # -wave_direction 1,0,0 \
13 # # # -save_postproc_mesh false \
14 # # # -my_order 4 \
15 # # # -my_max_post_proc_ref_level 0 \
16 # # # -amplitude_of_incident_wave 1 \
17 # # # -duffy true \
18 # # # -lobatto true \
19 # # # -fieldsplit_1_ksp_type fgmres \
20 # # # -fieldsplit_1_pc_type lu \
21 # # # -fieldsplit_1_pc_factor_mat_solver_package superlu_dist \
22 # # # -fieldsplit_1_ksp_max_it 30 \
23 # # # -fieldsplit_0_ksp_type fgmres \
24 # # # -fieldsplit_0_pc_type lu \
25 # # # -fieldsplit_0_pc_factor_mat_solver_package superlu_dist \
26 # # # -fieldsplit_0_ksp_max_it 30 \
27 # # # -pc_fieldsplit_type schur \
28 # # # -ksp_type fgmres \
29 # # # -ksp_atol 1e-8 \
30 # # # -ksp_rtol 1e-8 \
31 # # # -ksp_max_it 800 \
32 # # # -ksp_monitor \
33 # # # -add_incident_wave false 2>&1 |
34 # # # # echo $(grep --line-buffered -i "Problem ACOUSTIC_PROBLEM Nb. rows\|max_value_order_p\|min_value_order_p\|avg_value_order_p\|nz_used\|info.memory\|Time =\|Memory usage =" | sort | uniq) | tee -a -i hard_cylinder_field_split_${FIELD_SPLIT}_h_${ODISP}_k${WAVE_NUMBER}_lobatto_${LOBATTO}.txt
35 # # # # echo $(grep --line-buffered -i "Problem ACOUSTIC_PROBLEM Nb. rows\|max_value_order_p\|min_value_order_p\|avg_value_order_p\|nz_used\|info.memory\|Time =\|Memory usage =" | sort | uniq) | awk '{print $2,$4,$6,$8,$10,$15,$25,$32,$37}' | tee -a -i hard_cylinder_field_split_${FIELD_SPLIT}_h_${ODISP}_k${WAVE_NUMBER}_lobatto_${LOBATTO}.txt
36 # # # # echo $(grep --line-buffered -i "Problem ACOUSTIC_PROBLEM Nb. rows\|nz_used\|Time =\|Memory usage =" | sort | uniq) | awk '{print $2,$15,$22,$23,$27,$28}' | tee -a -i hard_cylinder_field_split_${FIELD_SPLIT}_h_${ODISP}_k${WAVE_NUMBER}_lobatto_${LOBATTO}.txt
37 # # # echo $(grep --line-buffered -i "Problem ACOUSTIC_PROBLEM Nb. rows\|max_value_order_p\|min_value_order_p\|avg_value_order_p\|nz_used\|Time =\|Memory usage =" | sort | uniq) | awk '{print $2,$5,$9,$16,$21}' | tee -a -i hard_cylinder_field_split_true_h_${ODISP}_k${WAVE_NUMBER}_lobatto_${LOBATTO}.txt
38 # # # # echo $(grep --line-buffered -i "Problem ACOUSTIC_PROBLEM Nb. rows\|max_value_order_p\|min_value_order_p\|avg_value_order_p\|nz_used\|Time =\|Memory usage =" | sort | uniq) | awk '{print $2,$7,$17,$24,$29}' | tee -a -i hard_cylinder_field_split_true_h_${ODISP}_k${WAVE_NUMBER}_lobatto_${LOBATTO}.txt
39 # # # # echo " "
40 # #
41 # #
42 # #
43 # #
44 # #
45 # #
46 # #
47 # # # run 10MHZ droplet problem, and impinging sphere problem in time domain
48 # #
49 # #
50 # # #
51 # # #
52 # # # # remove all the previous time dependant results.
53 # # rm -rf pressure_real_time_step_*.h5m
54 # # #
55 # #
56 # # # rm -rf droplet_pressure
57 # # #
58 # # #
59 # ~/build_moFEM1/usr/tools/mofem_part -my_file ./mix_droplet.cub -my_nparts 8
60 # # #
61 # # # # 4 MHz
62 # # # # mpirun -np 8 ./fe_approximation \
63 # # # # -my_file out.h5m \
64 # # # # -my_is_partitioned true \
65 # # # # -duffy true \
66 # # # # -adaptivit false \
67 # # # # -error_level 1 \
68 # # # # -my_order 6 \
69 # # # # -save_postproc_mesh true \
70 # # # # -wave_number 6283.1 \
71 # # # # -my_max_post_proc_ref_level 2 \
72 # # # # -add_incident_wave true \
73 # # # # -wave_direction 1,0,0 \
74 # # # # -wave_oscilation_direction -0.3756,-0.9268,0.000000001 \
75 # # # # -monochromatic_wave false \
76 # # # # -signal_length 0.001 \
77 # # # # -radiation_field true \
78 # # # # -frequency 4000000 \
79 # # # # -signal_duration 0.00000025 \
80 # # # # -velocity 1500 \
81 # # # # # -complex_wave_number 79.8387 \
82 # # # # -complex_wave_number 2768 \
83 # # # # -transmission_coefficient 0.1491 \
84 # # # # -rayleigh_wave false \
85 # # # # -density 998 \
86 # # # # -reynolds_stress true \
87 # # # # -radiation_force true \
88 # # # # -material_coefficient1 -1348.8 \
89 # # # # -ksp_atol 1e-8 \
90 # # # # -ksp_rtol 1e-8 \
91 # # # # -ksp_max_it 2000 \
92 # # # # -ksp_monitor \
93 # # # # -space_data acoustic_data.txt \
94 # # # # -nb_of_time_step 14 \
95 # # # # -material_coefficient2 \
96 # # # # -ksp_converged_reason \
97 # # # # -ksp_monitor_true_residual \
98 # # # # -pc_type lu \
99 # # # # -ksp_type fgmres \
100 # # # # -pc_factor_mat_solver_package superlu_dist \
101 # # # # -droplet_radius 0.00210 \
102 # # # # -attenuation 0.067077277922418
103 # # #
104 # # #
105 # mpirun -np 8 ./fe_approximation \
106 # -my_file out.h5m \
107 # -my_is_partitioned true \
108 # -ksp_final_residual \
109 # -ksp_monitor \
110 # -ksp_converged_reason \
111 # -pc_type lu \
112 # -ksp_type fgmres \
113 # -pc_factor_mat_solver_package superlu_dist \
114 # -duffy true \
115 # -adaptivit false \
116 # -error_level 1 \
117 # -my_order 6 \
118 # -save_postproc_mesh true \
119 # -wave_number 15707.96326794897 \
120 # -my_max_post_proc_ref_level 2 \
121 # -add_incident_wave true \
122 # -wave_direction 1,0,0 \
123 # -wave_oscilation_direction -0.3756,-0.9268,3.2051e-09 \
124 # -monochromatic_wave false \
125 # -signal_length 0.0004 \
126 # -signal_duration 0.0000001 \
127 # -frequency 10000000 \
128 # -velocity 1500 \
129 # -complex_wave_number 2768 \
130 # -transmission_coefficient 0.1491 \
131 # -rayleigh_wave false \
132 # -density 998 \
133 # -droplet_radius 0.00210 \
134 # -attenuation 0.265146221848709 \
135 # -reynolds_stress true \
136 # -radiation_force true \
137 # -material_coefficient1 -1348.8 \
138 # -ksp_atol 1e-8 \
139 # -ksp_rtol 1e-8 \
140 # -ksp_max_it 2000 \
141 # -space_data acoustic_data.txt \
142 # -nb_of_time_step 13 \
143 # -material_coefficient2 | tee log_droplet
144 # # #
145 # # #
146 # # # ./do_vtk1.sh pressure_real_time_step_*.h5m
147 # # # mbconvert reynolds_stress.h5m reynolds_stress.vtk
148 # # # mkdir droplet_pressure
149 # # # mv pressure_real_time_step_*.vtk droplet_pressure
150 # # # mv reynolds_stress.vtk droplet_pressure
151 # #
152 # #
153 # # # -complex_wave_number 199.5968 \
154 # #
155 # #
156 # # rm -rf sphere_pressure
157 # #
158 # # ~/build_moFEM1/usr/tools/mofem_part -my_file ./impinging_sphere.cub -my_nparts 8
159 # #
160 # #
161 # #
162 # #
163 # #
164 # # mpirun -np 8 ./fe_field_split \
165 # # -my_file out.h5m \
166 # # -my_is_partitioned true \
167 # # -duffy true \
168 # # -adaptivit false \
169 # # -error_level 1 \
170 # # -my_order 6 \
171 # # -save_postproc_mesh true \
172 # # -wave_number 5 \
173 # # -my_max_post_proc_ref_level 2 \
174 # # -add_incident_wave true \
175 # # -wave_direction 1,0,0 \
176 # # -monochromatic_wave false \
177 # # -space_data acoustic_data.txt \
178 # # -nb_of_time_step 12 \
179 # # -signal_length 1.2566 \
180 # # -signal_duration 0.0008377580409572781 \
181 # # -amplitude_of_incident_wave 1 \
182 # # -fieldsplit_1_ksp_type fgmres \
183 # # -fieldsplit_1_pc_type lu \
184 # # -fieldsplit_1_pc_factor_mat_solver_package mumps \
185 # # -fieldsplit_1_ksp_max_it 45 \
186 # # -fieldsplit_0_ksp_type fgmres \
187 # # -fieldsplit_0_pc_type lu \
188 # # -fieldsplit_0_pc_factor_mat_solver_package mumps \
189 # # -fieldsplit_0_ksp_max_it 45 \
190 # # -pc_fieldsplit_type schur \
191 # # -ksp_type fgmres \
192 # # -ksp_atol 1e-8 \
193 # # -ksp_rtol 1e-8 \
194 # # -ksp_max_it 800 \
195 # # -ksp_monitor \
196 # # -density 998 \
197 # # -reynolds_stress true | tee log_sphere
198 # #
199 # # # -wave_number 10 \
200 # # # -signal_length 0.6283 \
201 # # # -signal_duration 0.0004233827493261457 \
202 # #
203 # #
204 # # ./do_vtk1.sh pressure_real_time_step_*.h5m
205 # # mbconvert reynolds_stress.h5m reynolds_stress.vtk
206 # # mkdir sphere_pressure
207 # # mv reynolds_stress.vtk sphere_pressure
208 # # mv pressure_real_time_step_*.vtk sphere_pressure
209 # #
210 # #
211 # #
212 # #
213 # #
214 # #
215 # #
216 # #
217 # #
218 # #
219 # # # END
220 #
221 # # mbpart -t -p PARTKWAY 8 mix_droplet.cub mix_droplet_8.h5m
222 # #
223 # # ############### START UNIFORM ADAPTIVITY calculation ####################
224 # # ADAPTIVITY="false"
225 # #
226 # echo "Start uniform P enrichment numerical calculation with $NB_PROC processes, P Non-uniform ? ${ADAPTIVITY}. mesh mix_droplet.cub and order 7 ($COUNT) and wavenumber $WAVE_NUMBER, error level $ERRORLV"
227 # mpirun -np 8 ./fe_approximation \
228 # -my_file mix_droplet_8.h5m \
229 # -my_is_partitioned true \
230 # -wave_number 5 \
231 # -wave_direction 1,0,0 \
232 # -wave_oscilation_direction -0.3756,-0.9268,3.2051e-09 \
233 # -material_coefficient1 -1348.8 \
234 # -save_postproc_mesh false \
235 # -ksp_type fgmres \
236 # -pc_type lu \
237 # -pc_factor_mat_solver_package mumps \
238 # -ksp_monitor \
239 # -my_order 3 \
240 # -my_max_post_proc_ref_level 0 \
241 # -amplitude_of_incident_wave 1 \
242 # -duffy true \
243 # -lobatto true \
244 # -adaptivity false \
245 # -error_level 2 \
246 # -error_type 1 \
247 # -add_incident_wave false 2>&1 |
248 # echo $(grep --line-buffered -i "Problem ACOUSTIC_PROBLEM Nb. rows\|max_value_order_p\|min_value_order_p\|avg_value_order_p\|nz_used\|Time =\|Memory usage =" | sort | uniq)
249 # # | awk '{print $2,$9,$16,$21}'
250 # # echo $(grep --line-buffered -i "Problem ACOUSTIC_PROBLEM Nb. rows\|max_value_order_p\|min_value_order_p\|avg_value_order_p\|nz_used\|Time =\|Memory usage =" | sort | uniq) | awk '{print $2,$7,$17,$24,$29}' | tee -a -i p_refinement_adaptivity_${ADAPTIVITY}_hard_droplet_k_lobatto_${LOBATTO}_errorlv_${ERRORLV}_type_${ERROR_TYPE}.txt
251 #
252 # echo " "
253 # # echo " begin error calculation "
254 # echo " check if data is correct "
255 # #echo $(grep --line-buffered -i "Problem ACOUSTIC_PROBLEM Nb. rows\|nz_used\|Time =\|Memory usage =" | sort | uniq) | awk '{print $2,$7,$17,$24,$25,$29,$30}' | tee -a -i hard_droplet_k_${WAVE_NUMBER}_lobatto_${LOBATTO}_errorlv_${ERRORLV}_type_${ERROR_TYPE}.txt
256 # #echo $(grep --line-buffered -i "Problem ACOUSTIC_PROBLEM Nb. rows\|nz_used\|Time =\|Memory usage =" | sort | uniq) | tee -a -i hard_droplet_k_${WAVE_NUMBER}_lobatto_${LOBATTO}_errorlv_${ERRORLV}_type_${ERROR_TYPE}.txt
257 # #BEGIN error calculation awk '{print $6,$12}' grep first 10 lines |head -10|
258 # #grep --line-buffered -i "Problem ACOUSTIC_PROBLEM Nb. rows\|Time =\|Memory usage =" | tee -a -i hard_droplet_k_${WAVE_NUMBER}_lobatto_${LOBATTO}_errorlv_${ERRORLV}_type_${ERROR_TYPE}.txt
259 #
260 # # # mpirun -np 2 ./error_norm \
261 # # # -my_file ./fe_solution.h5m \
262 # # # -my_is_partitioned false \
263 # # # -norm_type $ERROR_TYPE \
264 # # # -relative_error false \
265 # # # -ksp_type fgmres \
266 # # # -pc_type lu \
267 # # # -pc_factor_mat_solver_package superlu_dist \
268 # # # -ksp_monitor \
269 # # # -my_order 1 \
270 # # # -my_max_post_proc_ref_level 0 \
271 # # # -save_postproc_mesh false 2>&1 | echo $(grep --line-buffered -i " realtive error ") | awk '{print $6,$12}' | tee -a -i p_refinement_adaptivity_${ADAPTIVITY}_error_hard_droplet_k_lobatto_${LOBATTO}_errorlv_${ERRORLV}_type_${ERROR_TYPE}.txt
272 # # #
273 #
274 #
275 # # mbpart -t -p PARTKWAY 8 dense_droplet.cub dense_droplet.h5m
276 # #
277 # #
278 # #
279 # #
280 # # echo "Start H refinement numerical calculation with $NB_PROC processes, order = 2, H Refinement mesh dense_droplet.h5m and order 2 and wavenumber $WAVE_NUMBER, error level $ERRORLV"
281 # # mpirun -np 8 ./fe_approximation \
282 # # -my_file dense_droplet.h5m \
283 # # -my_is_partitioned true \
284 # # -wave_number 5 \
285 # # -wave_direction 1,0,0 \
286 # # -wave_oscilation_direction -0.3756,-0.9268,3.2051e-09 \
287 # # -material_coefficient1 -1348.8 \
288 # # -save_postproc_mesh false \
289 # # -ksp_type fgmres \
290 # # -pc_type lu \
291 # # -pc_factor_mat_solver_package mumps \
292 # # -ksp_monitor \
293 # # -my_order 2 \
294 # # -my_max_post_proc_ref_level 0 \
295 # # -amplitude_of_incident_wave 1 \
296 # # -duffy true \
297 # # -lobatto true \
298 # # -adaptivity false \
299 # # -add_incident_wave false 2>&1 |
300 # # echo $(grep --line-buffered -i "Problem ACOUSTIC_PROBLEM Nb. rows\|max_value_order_p\|min_value_order_p\|avg_value_order_p\|nz_used\|Time =\|Memory usage =" | sort | uniq)
301 # # # echo $(grep --line-buffered -i "Problem ACOUSTIC_PROBLEM Nb. rows\|max_value_order_p\|min_value_order_p\|avg_value_order_p\|nz_used\|Time =\|Memory usage =" | sort | uniq) | awk '{print $2,$7,$17,$24,$29}' | tee -a -i h_refinement_hard_droplet_k_lobatto_${LOBATTO}_errorlv_${ERRORLV}_type_${ERROR_TYPE}.txt
302 # # # echo " "
303 # # # echo " begin error calculation "
304 # # echo " check if data is correct "
305 #
306 #
307 #
308 #
309 # mpirun -np 8 ./fe_approximation -my_file mix_droplet.cub -my_is_partitioned false -wave_number 6283.185307179586 -wave_direction 1,0,0 -wave_oscilation_direction -0.3756,-0.9268,3.2051e-09 -material_coefficient1 -1348.8 -save_postproc_mesh true -ksp_type fgmres -pc_type lu -pc_factor_mat_solver_package mumps -ksp_monitor -my_order 8 -my_max_post_proc_ref_level 2 -amplitude_of_incident_wave 1 -duffy true -lobatto true -adaptivity true -error_level 2 -error_type 2 -add_incident_wave false -complex_wave_number 79.8387 -frequency 4000000 -velocity 1500
310 
311 rm -rf fe_solution_mesh_post_proc_adaptivity_10mhz.vtk
312 
313 mpirun -np 8 ./fe_approximation \
314 -my_file mix_droplet.cub \
315 -my_is_partitioned false \
316 -wave_number 15707.96326794896 \
317 -wave_direction 1,0,0 \
318 -wave_oscilation_direction -0.3756,-0.9268,3.2051e-09 \
319 -material_coefficient1 -1348.8 \
320 -save_postproc_mesh true \
321 -ksp_type fgmres \
322 -pc_type lu \
323 -pc_factor_mat_solver_package mumps \
324 -ksp_monitor \
325 -my_order 1 \
326 -my_max_post_proc_ref_level 3 \
327 -amplitude_of_incident_wave 1 \
328 -duffy true \
329 -lobatto true \
330 -adaptivity true \
331 -error_level 2 \
332 -error_type 2 \
333 -add_incident_wave false \
334 -frequency 10000000 \
335 -velocity 1500 \
336 -complex_wave_number 199.5968
337 
338 mbconvert fe_solution_mesh_post_proc.h5m fe_solution_mesh_post_proc.vtk
339 mv fe_solution_mesh_post_proc.vtk fe_solution_mesh_post_proc_adaptivity_10mhz.vtk
340 
341 
342 
343 
344 rm -rf fe_solution_mesh_post_proc_p_10mhz.vtk
345 ~/build_moFEM1/usr/tools/mofem_part -my_file ./mix_droplet.cub -my_nparts 8
346 
347 mpirun -np 8 ./fe_approximation \
348 -my_file out.h5m \
349 -my_is_partitioned true \
350 -wave_number 15707.96326794896 \
351 -wave_direction 1,0,0 \
352 -wave_oscilation_direction -0.3756,-0.9268,3.2051e-09 \
353 -material_coefficient1 -1348.8 \
354 -save_postproc_mesh true \
355 -ksp_type fgmres \
356 -pc_type lu \
357 -pc_factor_mat_solver_package mumps \
358 -ksp_monitor \
359 -my_order 8 \
360 -my_max_post_proc_ref_level 3 \
361 -amplitude_of_incident_wave 1 \
362 -duffy true \
363 -lobatto true \
364 -adaptivity false \
365 -error_level 2 \
366 -error_type 2 \
367 -add_incident_wave true \
368 -frequency 10000000 \
369 -velocity 1500 \
370 -complex_wave_number 199.5968
371 
372 mbconvert fe_solution_mesh_post_proc.h5m fe_solution_mesh_post_proc.vtk
373 mv fe_solution_mesh_post_proc.vtk fe_solution_mesh_post_proc_p_10mhz.vtk
374 
375 
376 rm -rf fe_solution_mesh_post_proc_h_10mhz.vtk
377 ~/build_moFEM1/usr/tools/mofem_part -my_file ./dense_droplet.cub -my_nparts 8
378 
379 mpirun -np 8 ./fe_approximation \
380 -my_file out.h5m \
381 -my_is_partitioned true \
382 -wave_number 15707.96326794896 \
383 -wave_direction 1,0,0 \
384 -wave_oscilation_direction -0.3756,-0.9268,3.2051e-09 \
385 -material_coefficient1 -1348.8 \
386 -save_postproc_mesh true \
387 -ksp_type fgmres \
388 -pc_type lu \
389 -pc_factor_mat_solver_package mumps \
390 -ksp_monitor \
391 -my_order 2 \
392 -my_max_post_proc_ref_level 3 \
393 -amplitude_of_incident_wave 1 \
394 -duffy true \
395 -lobatto true \
396 -adaptivity false \
397 -error_level 2 \
398 -error_type l2 \
399 -add_incident_wave false \
400 -frequency 10000000 \
401 -velocity 1500 \
402 -complex_wave_number 199.5968
403 
404 mbconvert fe_solution_mesh_post_proc.h5m fe_solution_mesh_post_proc.vtk
405 mv fe_solution_mesh_post_proc.vtk fe_solution_mesh_post_proc_h_10mhz.vtk
406 
407 
408 
409 #
410 #
411 #
412 #
413 #
414 #
415 #
416 #
417 #
418 #
419 #
420 #
421 #
422 #
423 #
424 #
425 #
426 #
427 #
428 #
429 #
430 #
431 #
432 #
433 # ### END ###