[cig-commits] r18246 - short/3D/PyLith/trunk/libsrc/problems

brad at geodynamics.org brad at geodynamics.org
Tue Apr 19 09:16:56 PDT 2011


Author: brad
Date: 2011-04-19 09:16:56 -0700 (Tue, 19 Apr 2011)
New Revision: 18246

Modified:
   short/3D/PyLith/trunk/libsrc/problems/Solver.cc
Log:
Small fix. Use jacobianPC member for context.

Modified: short/3D/PyLith/trunk/libsrc/problems/Solver.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/problems/Solver.cc	2011-04-19 16:00:47 UTC (rev 18245)
+++ short/3D/PyLith/trunk/libsrc/problems/Solver.cc	2011-04-19 16:16:56 UTC (rev 18246)
@@ -141,8 +141,8 @@
 					   const topology::Jacobian& jacobian,
 					   const topology::SolutionFields& fields)
 { // _setupFieldSplit
-  assert(0 != pc);
-  assert(0 != formulation);
+  assert(pc);
+  assert(formulation);
 
   PetscErrorCode err = 0;
 
@@ -204,8 +204,11 @@
     // Set preconditioning matrix in formulation
     formulation->customPCMatrix(_jacobianPCFault);
 
+    assert(_jacobianPC);
+    assert(_jacobianPCFault);
+
     _ctx.pc = *pc;
-    _ctx.A = jacobian.matrix();
+    _ctx.A = _jacobianPC;
     _ctx.faultFieldName = "3";
     _ctx.faultA = _jacobianPCFault;
   } // if



More information about the CIG-COMMITS mailing list