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

[knepley at geodynamics.org]

Author: knepley
Date: 2009-06-10 17:37:29 -0700 (Wed, 10 Jun 2009)
New Revision: 15190

Modified:
   short/3D/PyLith/trunk/libsrc/problems/SolverLinear.cc
Log:
FieldSplit seems to be working


Modified: short/3D/PyLith/trunk/libsrc/problems/SolverLinear.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/problems/SolverLinear.cc	2009-06-10 23:54:19 UTC (rev 15189)
+++ short/3D/PyLith/trunk/libsrc/problems/SolverLinear.cc	2009-06-11 00:37:29 UTC (rev 15190)
@@ -83,6 +83,8 @@
 
     err = KSPGetPC(_ksp, &pc); CHECK_PETSC_ERROR(err);
     err = PCSetType(pc, PCFIELDSPLIT); CHECK_PETSC_ERROR(err);
+    err = PCSetOptionsPrefix(pc, "fs_"); CHECK_PETSC_ERROR(err);
+    err = PCSetFromOptions(pc); CHECK_PETSC_ERROR(err);
 #if defined(FIELD_SPLIT)
     constructFieldSplit(residual.section(), sieveMesh->getFactory()->getGlobalOrder(sieveMesh, "default", residual.section()), residual.vector(), pc);
 #endif