[cig-commits] r20961 - in seismo/3D/SPECFEM3D/trunk: examples/homogeneous_halfspace_HEX8 src/decompose_mesh src/shared

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-10-27 15:47:29 -0700 (Sat, 27 Oct 2012)
New Revision: 20961

Added:
   seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8/plotall_comparison_seismograms.gnu
Modified:
   seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/Makefile.in
   seismo/3D/SPECFEM3D/trunk/src/shared/check_mesh_resolution.f90
Log:
now mention the fact that element imbalance can be normal if fluid, elastic and poro elements are mixed (also with CPML elements).
also added a gnuplot script for the HEX27 example.


Added: seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8/plotall_comparison_seismograms.gnu
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8/plotall_comparison_seismograms.gnu	                        (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8/plotall_comparison_seismograms.gnu	2012-10-27 22:47:29 UTC (rev 20961)
@@ -0,0 +1,22 @@
+
+set xlabel "time (s)"
+set ylabel "displacement (m)"
+
+plot 'REF_SEIS/X1.DB.BXZ.semd' w l lc 1, 'OUTPUT_FILES/X1.DB.BXZ.semd' w l lc 3
+pause -1 "hit key..."
+
+plot 'REF_SEIS/X10.DB.BXZ.semd' w l lc 1, 'OUTPUT_FILES/X10.DB.BXZ.semd' w l lc 3
+pause -1 "hit key..."
+
+plot 'REF_SEIS/X20.DB.BXZ.semd' w l lc 1, 'OUTPUT_FILES/X20.DB.BXZ.semd' w l lc 3
+pause -1 "hit key..."
+
+plot 'REF_SEIS/X30.DB.BXZ.semd' w l lc 1, 'OUTPUT_FILES/X30.DB.BXZ.semd' w l lc 3
+pause -1 "hit key..."
+
+plot 'REF_SEIS/X40.DB.BXZ.semd' w l lc 1, 'OUTPUT_FILES/X40.DB.BXZ.semd' w l lc 3
+pause -1 "hit key..."
+
+plot 'REF_SEIS/X50.DB.BXZ.semd' w l lc 1, 'OUTPUT_FILES/X50.DB.BXZ.semd' w l lc 3
+pause -1 "hit key..."
+

Modified: seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/Makefile.in
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/Makefile.in	2012-10-27 15:51:11 UTC (rev 20960)
+++ seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/Makefile.in	2012-10-27 22:47:29 UTC (rev 20961)
@@ -54,7 +54,7 @@
 @COND_MPI_FALSE at COND_MPI_CPPFLAGS =
 
 SCOTCH_LIBS = -L"@SCOTCH_LIBDIR@" -lscotch -lscotcherr
-### DK DK added support for METIS as well
+### added support for METIS as well, thus uncomment the line below and compile METIS if you want to use it instead of SCOTCH
 #SCOTCH_LIBS = -L"@SCOTCH_LIBDIR@" -lscotch -lscotcherr -L./src/decompose_mesh/metis-4.0.3 -L./metis-4.0.3 -lmetis
 SCOTCH_INC = -I"@SCOTCH_INCLUDEDIR@"
 

Modified: seismo/3D/SPECFEM3D/trunk/src/shared/check_mesh_resolution.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/shared/check_mesh_resolution.f90	2012-10-27 15:51:11 UTC (rev 20960)
+++ seismo/3D/SPECFEM3D/trunk/src/shared/check_mesh_resolution.f90	2012-10-27 22:47:29 UTC (rev 20961)
@@ -287,6 +287,11 @@
                       ' = ',sngl((dble(NGLOB_AB_global_max) / dble(NGLOB_AB_global_min) - 1.d0) * 100.d0),' %'
     write(IMAIN,*) 'NGLOB_global_sum = ',NGLOB_AB_global_sum
     write(IMAIN,*)
+    write(IMAIN,*) 'If you have elements of a single type (all acoustic, all elastic, all poroelastic, and without CPML)'
+    write(IMAIN,*) 'in the whole mesh, then there should be no significant imbalance in the above numbers.'
+    write(IMAIN,*) 'Otherwise, it is normal to have imbalance in elements and points because the domain decomposer'
+    write(IMAIN,*) 'compensates for the different cost of different elements by partitioning them unevenly among processes.'
+    write(IMAIN,*)
     write(IMAIN,*) '********'
     write(IMAIN,*) 'Model: P velocity min,max = ',vpmin_glob,vpmax_glob
     write(IMAIN,*) 'Model: S velocity min,max = ',vsmin_glob,vsmax_glob