[cig-commits] r16482 - short/3D/PyLith/trunk/libsrc/faults

[brad at geodynamics.org]

Author: brad
Date: 2010-04-01 11:53:41 -0700 (Thu, 01 Apr 2010)
New Revision: 16482

Modified:
   short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc
Log:
Set tighter default tolerances for sensitivity solve.

Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc	2010-04-01 17:59:12 UTC (rev 16481)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc	2010-04-01 18:53:41 UTC (rev 16482)
@@ -1514,6 +1514,16 @@
     } // if
     err = KSPCreate(_faultMesh->comm(), &_ksp); CHECK_PETSC_ERROR(err);
     err = KSPSetInitialGuessNonzero(_ksp, PETSC_FALSE); CHECK_PETSC_ERROR(err);
+    double rtol = 0.0;
+    double atol = 0.0;
+    double dtol = 0.0;
+    int maxIters = 0;
+    err = KSPGetTolerances(_ksp, &rtol, &atol, &dtol, &maxIters); 
+    CHECK_PETSC_ERROR(err);
+    rtol = 1.0e-15;
+    atol = 1.0e-25;
+    err = KSPSetTolerances(_ksp, rtol, atol, dtol, maxIters);
+    CHECK_PETSC_ERROR(err);
 
     PC pc;
     err = KSPGetPC(_ksp, &pc); CHECK_PETSC_ERROR(err);
@@ -1639,6 +1649,8 @@
   } // for
 
   _jacobian->assemble("final_assembly");
+
+  //_jacobian->view(); // DEBUGGING
 } // _sensitivityUpdateJacobian
 
 // ----------------------------------------------------------------------