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

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-02-19 15:14:17 -0800 (Sun, 19 Feb 2012)
New Revision: 19653

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/todo_list_please_dont_remove.txt
Log:
updated 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-02-19 23:10:59 UTC (rev 19652)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/todo_list_please_dont_remove.txt	2012-02-19 23:14:17 UTC (rev 19653)
@@ -4,7 +4,8 @@
 
 Things that could be done in a future version:
 
-- use a better checkpointing/restarting system with CRC-32 as in SPECFEM3D_GLOBE/tags/v4.1.0_beta_merged_mesher_solver_non_blocking_MPI/src
+- 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
 
 - use a more accurate time scheme: either symplectic schemes (but there could be a problem in the viscoelastic case, see below) or else Runge Kutta RK4, which is non staggered and would thus solve the approximate implementation of the displacement gradient that I used to implement attenuation back in 1999. It is not clear if symplectic schemes could be suitable for viscoelastic media (Tarje implemented them for elastic media only), but RK4 is definitely suitable (Roland Martin used it successfully in 2D viscoelastic).
 
@@ -107,13 +108,13 @@
 It would be nice having both improved formulations included for those cases.
 This seems to be a nice outlook for a future version update.
 
-- use a potential of (rho * u) instead of u in the fluid, in case of
-  fluid-fluid discontinuities
+- use a better checkpointing/restarting system with CRC-32 as in SPECFEM3D_GLOBE/tags/v4.1.0_beta_merged_mesher_solver_non_blocking_MPI/src
 
+- use a potential of (rho * u) instead of u in the fluid, in case of fluid-fluid discontinuities
+
 - merge SPECFEM3D_GLOBE into SPECFEM3D, i.e. create a single package and see the mesh of the Earth as a general mesh that we would decompose in parallel using PT-Scotch. Problem: would this scale for huge global Earth simulations, e.g. at a seismic period of 1 second on 200,000 processor cores?
 
-- make the code compatible with helioseismology / general Cowling formulation (will be done
-  with Tarje Nissen-Meyer)
+- make the code compatible with helioseismology / general Cowling formulation (could be done with Tarje Nissen-Meyer)
 
 - could be done by Vala:
 we could use heuristic rules to make source and receiver detection much