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

brad at geodynamics.org brad at geodynamics.org
Wed May 19 16:48:26 PDT 2010


Author: brad
Date: 2010-05-19 16:48:26 -0700 (Wed, 19 May 2010)
New Revision: 16752

Modified:
   short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc
Log:
Fixed build errors.

Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc	2010-05-19 22:55:08 UTC (rev 16751)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc	2010-05-19 23:48:26 UTC (rev 16752)
@@ -626,10 +626,10 @@
   const PetscMat jacobianMatrix = jacobian->matrix();
   assert(0 != jacobianMatrix);
 
-  const ALE::Obj<SieveMesh::order_type>& faultGlobalOrder =
+  const ALE::Obj<SieveMesh::order_type>& lagrangeGlobalOrder =
       sieveMesh->getFactory()->getGlobalOrder(sieveMesh, "faultDefault",
-        solutionSection, numKSP-1);
-  assert(!faultGlobalOrder.isNull());
+        solutionSection, spaceDim);
+  assert(!lagrangeGlobalOrder.isNull());
 
   _logger->eventEnd(setupEvent);
 #if !defined(DETAILED_EVENT_LOGGING)
@@ -692,11 +692,11 @@
 #endif
 
     // Set global preconditioner index associated with Lagrange constraint vertex.
-    const int indexLprecond = faultGlobalOrder->getIndex(v_lagrange);
+    const int indexLprecond = lagrangeGlobalOrder->getIndex(v_lagrange);
     
     // Set diagonal entries in preconditioned matrix.
     for (int iDim=0; iDim < spaceDim; ++iDim)
-      MatSetValue(precondMatrix,
+      MatSetValue(*precondMatrix,
 		  indexLprecond + iDim,
 		  indexLprecond + iDim,
 		  precondVertexL[iDim],



More information about the CIG-COMMITS mailing list