[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