[cig-commits] r20702 - short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/feassemble

[brad at geodynamics.org]

Author: brad
Date: 2012-09-07 08:52:41 -0700 (Fri, 07 Sep 2012)
New Revision: 20702

Modified:
   short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/feassemble/GeometryTet3D.cc
Log:
Setup empirical scaling of radius of inscribed sphere of tet for stable time step.

Modified: short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/feassemble/GeometryTet3D.cc
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/feassemble/GeometryTet3D.cc	2012-09-07 15:42:19 UTC (rev 20701)
+++ short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/feassemble/GeometryTet3D.cc	2012-09-07 15:52:41 UTC (rev 20702)
@@ -284,9 +284,10 @@
     faceArea(coordinatesCell, 2) +
     faceArea(coordinatesCell, 3);
     
-  const PylithScalar r = 3.0* v / a;
-  if (r < minWidth) {
-    minWidth = r;
+  const PylithScalar r = 3.0 * v / a;
+  const PylithScalar rwidth = 3.0*r; // based on empirical tests
+  if (rwidth < minWidth) {
+    minWidth = rwidth;
   } // if
 
   PetscLogFlops(5);