[cig-commits] r18841 - short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults

[brad at geodynamics.org]

Author: brad
Date: 2011-08-19 18:10:07 -0700 (Fri, 19 Aug 2011)
New Revision: 18841

Modified:
   short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults/FaultCohesiveDyn.cc
Log:
Raise tolerances for friction solve to more reasonable levels.

Modified: short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults/FaultCohesiveDyn.cc	2011-08-20 01:09:00 UTC (rev 18840)
+++ short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults/FaultCohesiveDyn.cc	2011-08-20 01:10:07 UTC (rev 18841)
@@ -438,10 +438,10 @@
 
 #if 0 // DEBUGGING
   slipSection->view("SLIP");
-  slipRateSection->view("SLIP RATE");
-  areaSection->view("AREA");
-  dispTSection->view("DISP (t)");
-  dispTIncrSection->view("DISP INCR (t->t+dt)");
+  //slipRateSection->view("SLIP RATE");
+  //areaSection->view("AREA");
+  //dispTSection->view("DISP (t)");
+  //dispTIncrSection->view("DISP INCR (t->t+dt)");
 #endif
 
   const int numVertices = _cohesiveVertices.size();
@@ -596,9 +596,9 @@
 
 #if 0 // DEBUGGING
   dLagrangeTpdtSection->view("AFTER dLagrange");
-  dispTIncrSection->view("AFTER DISP INCR (t->t+dt)");
+  //dispTIncrSection->view("AFTER DISP INCR (t->t+dt)");
   slipSection->view("AFTER SLIP");
-  slipRateSection->view("AFTER SLIP RATE");
+  //slipRateSection->view("AFTER SLIP RATE");
 #endif
 } // constrainSolnSpace
 
@@ -1581,8 +1581,8 @@
     int maxIters = 0;
     err = KSPGetTolerances(_ksp, &rtol, &atol, &dtol, &maxIters); 
     CHECK_PETSC_ERROR(err);
-    rtol = 1.0e-15;
-    atol = 1.0e-25;
+    rtol = 1.0e-12;
+    atol = 1.0e-15;
     err = KSPSetTolerances(_ksp, rtol, atol, dtol, maxIters);
     CHECK_PETSC_ERROR(err);