[cig-commits] [commit] devel: fixed confusing sentences about the number of elements and points in the mesh (problem detected by Vadim Monteiller) (afa264d)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem3d

On branch  : devel
Link       : https://github.com/geodynamics/specfem3d/compare/aad4eb97227c36e1cd6293a2f0109990b4e5b30b...afa264da6cf450a0e4775808bd59cc84d9b4fcc8

>---------------------------------------------------------------

commit afa264da6cf450a0e4775808bd59cc84d9b4fcc8
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Tue Sep 16 22:08:03 2014 +0200

    fixed confusing sentences about the number of elements and points in the mesh (problem detected by Vadim Monteiller)


>---------------------------------------------------------------

afa264da6cf450a0e4775808bd59cc84d9b4fcc8
 src/generate_databases/finalize_databases.f90 |  8 ++++----
 src/meshfem3D/meshfem3D.f90                   | 11 ++++++-----
 2 files changed, 10 insertions(+), 9 deletions(-)

diff --git a/src/generate_databases/finalize_databases.f90 b/src/generate_databases/finalize_databases.f90
index ac84660..2912a61 100644
--- a/src/generate_databases/finalize_databases.f90
+++ b/src/generate_databases/finalize_databases.f90
@@ -48,16 +48,16 @@
     write(IMAIN,*) 'Repartition of elements:'
     write(IMAIN,*) '-----------------------'
     write(IMAIN,*)
-    write(IMAIN,*) 'total number of elements in each slice: ',NSPEC_AB
-    write(IMAIN,*) 'total number of points in each slice: ',NGLOB_AB
+    write(IMAIN,*) 'total number of elements in mesh slice 0: ',NSPEC_AB
+    write(IMAIN,*) 'total number of points in mesh slice 0: ',NGLOB_AB
     write(IMAIN,*)
     write(IMAIN,*) 'total number of elements in entire mesh: ',nspec_total
 ! the float() statement below are for the case of more than 2 Gigapoints per mesh, in which
 ! case and integer(kind=4) counter would overflow and display an incorrect negative value;
 ! converting it to float fixes the problem (but then prints some extra decimals equal to zero).
 ! Another option would be to declare the sum as integer(kind=8) and then print it.
-    write(IMAIN,*) 'total number of points in entire mesh: ',nglob_total
-    write(IMAIN,*) 'total number of DOFs in entire mesh: ',nglob_total*dble(NDIM)
+    write(IMAIN,*) 'approximate total number of points in entire mesh (with duplicates on MPI edges): ',nglob_total
+    write(IMAIN,*) 'approximate total number of DOFs in entire mesh (with duplicates on MPI edges): ',nglob_total*dble(NDIM)
     write(IMAIN,*)
     write(IMAIN,*) 'total number of time steps in the solver will be: ',NSTEP
     write(IMAIN,*)
diff --git a/src/meshfem3D/meshfem3D.f90 b/src/meshfem3D/meshfem3D.f90
index 4f0ad1b..56c08a7 100644
--- a/src/meshfem3D/meshfem3D.f90
+++ b/src/meshfem3D/meshfem3D.f90
@@ -857,17 +857,18 @@
     write(IMAIN,*) 'Repartition of elements:'
     write(IMAIN,*) '-----------------------'
     write(IMAIN,*)
-    write(IMAIN,*) 'total number of elements in each slice: ',NSPEC_AB
-    write(IMAIN,*)
-    write(IMAIN,*) 'total number of points in each slice: ',NGLOB_AB
+    write(IMAIN,*) 'total number of elements in mesh slice 0: ',NSPEC_AB
+    write(IMAIN,*) 'total number of points in mesh slice 0: ',NGLOB_AB
     write(IMAIN,*)
     write(IMAIN,*) 'total number of elements in entire mesh: ',NSPEC_AB*NPROC
 ! the float() statements here are for the case of more than 2 Gigapoints per mesh, in which
 ! case and integer(kind=4) counter would overflow and display an incorrect negative value;
 ! converting it to float fixes the problem (but then prints some extra decimals equal to zero).
 ! Another option would be to declare the sum as integer(kind=8) and then print it.
-    write(IMAIN,*) 'total number of points in entire mesh: ',dble(NGLOB_AB)*dble(NPROC)
-    write(IMAIN,*) 'total number of DOFs in entire mesh: ',dble(NGLOB_AB)*dble(NPROC*NDIM)
+    write(IMAIN,*) 'approximate total number of points in entire mesh (with duplicates on MPI edges): ', &
+                                             dble(NGLOB_AB)*dble(NPROC)
+    write(IMAIN,*) 'approximate total number of DOFs in entire mesh (with duplicates on MPI edges): ', &
+                                             dble(NGLOB_AB)*dble(NPROC*NDIM)
     write(IMAIN,*)
     ! write information about precision used for floating-point operations
     if(CUSTOM_REAL == SIZE_REAL) then