[cig-commits] r15489 - in short/3D/PyLith/trunk: libsrc/materials pylith/problems

willic3 at geodynamics.org willic3 at geodynamics.org
Thu Jul 30 15:55:21 PDT 2009


Author: willic3
Date: 2009-07-30 15:55:21 -0700 (Thu, 30 Jul 2009)
New Revision: 15489

Modified:
   short/3D/PyLith/trunk/libsrc/materials/GenMaxwellIsotropic3D.cc
   short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc
   short/3D/PyLith/trunk/pylith/problems/TimeStepAdapt.py
Log:
Changed stable time step to 0.2 * maxwellTime for Maxwell and generalized
Maxwell viscoelastic models, and changed default stability factor for
adaptive time stepping to 2.0.  This will give a default time step size
of 0.1 * maxwellTime for adaptive time stepping.



Modified: short/3D/PyLith/trunk/libsrc/materials/GenMaxwellIsotropic3D.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/materials/GenMaxwellIsotropic3D.cc	2009-07-30 21:33:24 UTC (rev 15488)
+++ short/3D/PyLith/trunk/libsrc/materials/GenMaxwellIsotropic3D.cc	2009-07-30 22:55:21 UTC (rev 15489)
@@ -830,7 +830,7 @@
   const int numMaxwellModels = _GenMaxwellIsotropic3D::numMaxwellModels;
   for (int i=0; i < numMaxwellModels; ++i) {
     const double maxwellTime = properties[p_maxwellTime+i];
-    const double dt = 0.1*maxwellTime;
+    const double dt = 0.2*maxwellTime;
     if (dt < dtStable)
       dtStable = dt;
   } // for

Modified: short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc	2009-07-30 21:33:24 UTC (rev 15488)
+++ short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc	2009-07-30 22:55:21 UTC (rev 15489)
@@ -695,7 +695,7 @@
   assert(_numVarsQuadPt == numStateVars);
 
   const double maxwellTime = properties[p_maxwellTime];
-  const double dtStable = 0.1 * maxwellTime;
+  const double dtStable = 0.2 * maxwellTime;
 
   return dtStable;
 } // _stableTimeStepImplicit

Modified: short/3D/PyLith/trunk/pylith/problems/TimeStepAdapt.py
===================================================================
--- short/3D/PyLith/trunk/pylith/problems/TimeStepAdapt.py	2009-07-30 21:33:24 UTC (rev 15488)
+++ short/3D/PyLith/trunk/pylith/problems/TimeStepAdapt.py	2009-07-30 22:55:21 UTC (rev 15489)
@@ -65,7 +65,7 @@
     adaptSkip.meta['tip'] = "Number of time steps to skip between " \
         "adjusting value."
 
-    stabilityFactor = pyre.inventory.float("stability_factor", default=1.2,
+    stabilityFactor = pyre.inventory.float("stability_factor", default=2.0,
                                     validator=pyre.inventory.greater(0.0))
     stabilityFactor.meta['tip'] = "'Safety factor' for stable time step."
 



More information about the CIG-COMMITS mailing list