[cig-commits] r19921 - seismo/3D/SPECFEM3D_GLOBE/trunk

dkomati1 at geodynamics.org dkomati1 at geodynamics.org
Wed Apr 4 15:24:31 PDT 2012

Author: dkomati1
Date: 2012-04-04 15:24:31 -0700 (Wed, 04 Apr 2012)
New Revision: 19921

added more details regarding 3D attenuation to the todo list

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/todo_list_please_dont_remove.txt
--- seismo/3D/SPECFEM3D_GLOBE/trunk/todo_list_please_dont_remove.txt	2012-04-04 22:22:43 UTC (rev 19920)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/todo_list_please_dont_remove.txt	2012-04-04 22:24:31 UTC (rev 19921)
@@ -7,6 +7,7 @@
 - fix the 3D attenuation memory size problem by making it optional and off by default.
 To do that, the following command will be very useful:
 nm --print-size --size-sort --reverse-sort --radix=d ./bin/xspecfem3D | cut -f 2- -d ' ' | less 
+I think switching from (NGLLX,NGLLY,NGLLZ) to (1,1,1) i.e. using the lower left corner of each element in the case of 1D attenuation would be easy to implement and sufficient to fix the problem (since array sizes would decrease by a factor of 125).
 - check if sigma_xz is inverted with sigma_zx and so on in compute_forces; usually this does not matter because the stress tensor is symmetric, but for Cosser
 at or for some C-PML formulations this could matter

More information about the CIG-COMMITS mailing list