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

[brad at geodynamics.org]

Author: brad
Date: 2011-04-19 17:10:43 -0700 (Tue, 19 Apr 2011)
New Revision: 18252

Modified:
   short/3D/PyLith/trunk/libsrc/problems/Solver.cc
Log:
Added missing return statement.

Modified: short/3D/PyLith/trunk/libsrc/problems/Solver.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/problems/Solver.cc	2011-04-19 23:55:43 UTC (rev 18251)
+++ short/3D/PyLith/trunk/libsrc/problems/Solver.cc	2011-04-20 00:10:43 UTC (rev 18252)
@@ -38,11 +38,11 @@
 typedef pylith::topology::Mesh::RealSection RealSection;
 
 EXTERN_C_BEGIN
-PetscErrorCode  MyMatGetSubMatrix(Mat mat, IS isrow, IS iscol, MatReuse reuse, Mat *newmat) {
+PetscErrorCode  MyMatGetSubMatrix(Mat mat, IS isrow, IS iscol, MatReuse reuse, Mat *newmat) { // MyMatGetSubMatrix
   FaultPreconCtx *ctx;
   IS              faultIS;
   PetscBool       isFaultRow, isFaultCol;
-  PetscErrorCode  ierr;
+  PetscErrorCode  ierr = 0;
 
   ierr = MatShellGetContext(mat, (void **) &ctx);CHKERRQ(ierr);
   ierr = PCFieldSplitGetIS(ctx->pc, ctx->faultFieldName, &faultIS);CHKERRQ(ierr);
@@ -52,11 +52,13 @@
     if (reuse == MAT_INITIAL_MATRIX) {
       ierr = PetscObjectReference((PetscObject) ctx->faultA);CHKERRQ(ierr);
       *newmat = ctx->faultA;
-    }
+    } // if
   } else {
     ierr = MatGetSubMatrix(ctx->A, isrow, iscol, reuse, newmat);CHKERRQ(ierr);
-  }
-}
+  } // if/else
+
+  return 0;
+} // MyMatGetSubMatrix
 EXTERN_C_END
 
 // ----------------------------------------------------------------------