[cig-commits] r21335 - in short/3D/PyLith/trunk: libsrc/pylith/faults unittests/libtests/faults unittests/libtests/meshio
knepley at geodynamics.org
knepley at geodynamics.org
Tue Feb 5 05:15:49 PST 2013
Author: knepley
Date: 2013-02-05 05:15:49 -0800 (Tue, 05 Feb 2013)
New Revision: 21335
Modified:
short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc
short/3D/PyLith/trunk/unittests/libtests/faults/TestBruneSlipFn.cc
short/3D/PyLith/trunk/unittests/libtests/faults/TestConstRateSlipFn.cc
short/3D/PyLith/trunk/unittests/libtests/faults/TestEqKinSrc.cc
short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveDyn.cc
short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveImpulses.cc
short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc
short/3D/PyLith/trunk/unittests/libtests/faults/TestLiuCosSlipFn.cc
short/3D/PyLith/trunk/unittests/libtests/faults/TestStepSlipFn.cc
short/3D/PyLith/trunk/unittests/libtests/faults/TestTimeHistorySlipFn.cc
short/3D/PyLith/trunk/unittests/libtests/faults/TestTractPerturbation.cc
short/3D/PyLith/trunk/unittests/libtests/meshio/TestDataWriterSubMesh.cc
Log:
Replaced subpointMap with subpointIS in tests
Modified: short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc 2013-02-05 13:15:19 UTC (rev 21334)
+++ short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc 2013-02-05 13:15:49 UTC (rev 21335)
@@ -1007,10 +1007,13 @@
faultMesh->coordsys(mesh);
+ DM dmMesh = mesh.dmMesh();
+ assert(dmMesh);
+
+
+
const ALE::Obj<SieveMesh>& sieveMesh = mesh.sieveMesh();
assert(!sieveMesh.isNull());
- DM dmMesh = mesh.dmMesh();
- assert(dmMesh);
ALE::Obj<SieveSubMesh>& faultSieveMesh = faultMesh->sieveMesh();
const ALE::Obj<SieveMesh::sieve_type>& sieve = sieveMesh->getSieve();
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestBruneSlipFn.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestBruneSlipFn.cc 2013-02-05 13:15:19 UTC (rev 21334)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestBruneSlipFn.cc 2013-02-05 13:15:49 UTC (rev 21335)
@@ -366,7 +366,7 @@
const ALE::Obj<RealSection>& oldCoordSection = sieveMesh->getRealSection("coordinates");
faultSieveMesh->setRealSection("coordinates", oldCoordSection);
DM dmMesh = faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscSection coordSection;
PetscInt vStart, vEnd;
@@ -374,7 +374,7 @@
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
err = DMPlexGetCoordinateSection(dmMesh, &coordSection);CHECK_PETSC_ERROR(err);
err = PetscSectionSetChart(coordSection, vStart, vEnd);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
@@ -389,7 +389,7 @@
err = VecCreate(mesh->comm(), &coordVec);CHECK_PETSC_ERROR(err);
err = VecSetSizes(coordVec, coordSize, PETSC_DETERMINE);CHECK_PETSC_ERROR(err);
err = VecSetFromOptions(coordVec);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt off;
@@ -400,7 +400,8 @@
coords[off+d] = oldCoords[d];
}
}
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
err = DMSetCoordinatesLocal(dmMesh, coordVec);CHECK_PETSC_ERROR(err);
@@ -476,7 +477,7 @@
const ALE::Obj<RealSection>& oldCoordSection = sieveMesh->getRealSection("coordinates");
faultSieveMesh->setRealSection("coordinates", oldCoordSection);
DM dmMesh = faultMesh.dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscSection coordSection;
PetscInt vStart, vEnd;
@@ -484,7 +485,7 @@
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
err = DMPlexGetCoordinateSection(dmMesh, &coordSection);CHECK_PETSC_ERROR(err);
err = PetscSectionSetChart(coordSection, vStart, vEnd);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
@@ -499,7 +500,7 @@
err = VecCreate(mesh.comm(), &coordVec);CHECK_PETSC_ERROR(err);
err = VecSetSizes(coordVec, coordSize, PETSC_DETERMINE);CHECK_PETSC_ERROR(err);
err = VecSetFromOptions(coordVec);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt off;
@@ -510,7 +511,8 @@
coords[off+d] = oldCoords[d];
}
}
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
err = DMSetCoordinatesLocal(dmMesh, coordVec);CHECK_PETSC_ERROR(err);
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestConstRateSlipFn.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestConstRateSlipFn.cc 2013-02-05 13:15:19 UTC (rev 21334)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestConstRateSlipFn.cc 2013-02-05 13:15:49 UTC (rev 21335)
@@ -297,7 +297,7 @@
const ALE::Obj<RealSection>& oldCoordSection = sieveMesh->getRealSection("coordinates");
faultSieveMesh->setRealSection("coordinates", oldCoordSection);
DM dmMesh = faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscSection coordSection;
PetscInt vStart, vEnd;
@@ -305,7 +305,7 @@
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
err = DMPlexGetCoordinateSection(dmMesh, &coordSection);CHECK_PETSC_ERROR(err);
err = PetscSectionSetChart(coordSection, vStart, vEnd);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
@@ -320,7 +320,7 @@
err = VecCreate(mesh->comm(), &coordVec);CHECK_PETSC_ERROR(err);
err = VecSetSizes(coordVec, coordSize, PETSC_DETERMINE);CHECK_PETSC_ERROR(err);
err = VecSetFromOptions(coordVec);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt off;
@@ -331,7 +331,8 @@
coords[off+d] = oldCoords[d];
}
}
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
err = DMSetCoordinatesLocal(dmMesh, coordVec);CHECK_PETSC_ERROR(err);
@@ -401,7 +402,7 @@
const ALE::Obj<RealSection>& oldCoordSection = sieveMesh->getRealSection("coordinates");
faultSieveMesh->setRealSection("coordinates", oldCoordSection);
DM dmMesh = faultMesh.dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscSection coordSection;
PetscInt vStart, vEnd;
@@ -409,7 +410,7 @@
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
err = DMPlexGetCoordinateSection(dmMesh, &coordSection);CHECK_PETSC_ERROR(err);
err = PetscSectionSetChart(coordSection, vStart, vEnd);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
@@ -424,7 +425,7 @@
err = VecCreate(mesh.comm(), &coordVec);CHECK_PETSC_ERROR(err);
err = VecSetSizes(coordVec, coordSize, PETSC_DETERMINE);CHECK_PETSC_ERROR(err);
err = VecSetFromOptions(coordVec);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt off;
@@ -435,7 +436,8 @@
coords[off+d] = oldCoords[d];
}
}
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
err = DMSetCoordinatesLocal(dmMesh, coordVec);CHECK_PETSC_ERROR(err);
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestEqKinSrc.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestEqKinSrc.cc 2013-02-05 13:15:19 UTC (rev 21334)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestEqKinSrc.cc 2013-02-05 13:15:49 UTC (rev 21335)
@@ -203,7 +203,7 @@
const ALE::Obj<RealSection>& oldCoordSection = sieveMesh->getRealSection("coordinates");
faultSieveMesh->setRealSection("coordinates", oldCoordSection);
DM dmMesh = faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscSection coordSection;
PetscInt vStart, vEnd;
@@ -211,7 +211,7 @@
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
err = DMPlexGetCoordinateSection(dmMesh, &coordSection);CHECK_PETSC_ERROR(err);
err = PetscSectionSetChart(coordSection, vStart, vEnd);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
@@ -226,7 +226,7 @@
err = VecCreate(mesh->comm(), &coordVec);CHECK_PETSC_ERROR(err);
err = VecSetSizes(coordVec, coordSize, PETSC_DETERMINE);CHECK_PETSC_ERROR(err);
err = VecSetFromOptions(coordVec);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt off;
@@ -237,7 +237,8 @@
coords[off+d] = oldCoords[d];
}
}
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
err = DMSetCoordinatesLocal(dmMesh, coordVec);CHECK_PETSC_ERROR(err);
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveDyn.cc 2013-02-05 13:15:19 UTC (rev 21334)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveDyn.cc 2013-02-05 13:15:49 UTC (rev 21335)
@@ -138,17 +138,17 @@
_initialize(&mesh, &fault, &fields);
DM dmMesh = fault._faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscInt vStart, vEnd, numPoints;
PetscErrorCode err;
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
- CPPUNIT_ASSERT(subpointMap);
- err = ISGetSize(subpointMap, &numPoints);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
+ CPPUNIT_ASSERT(subpointIS);
+ err = ISGetSize(subpointIS, &numPoints);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt faultPoint;
@@ -156,7 +156,8 @@
CPPUNIT_ASSERT(faultPoint >= 0);
CPPUNIT_ASSERT_EQUAL(faultPoint, v);
} // for
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
CPPUNIT_ASSERT_EQUAL(_data->numConstraintVert, vEnd-vStart);
// Check orientation
@@ -595,16 +596,16 @@
fault._calcTractions(&tractions, fields.get("disp(t)"));
DM faultDMMesh = fault._faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscInt vStart, vEnd, numPoints;
PetscErrorCode err;
CPPUNIT_ASSERT(faultDMMesh);
err = DMPlexGetDepthStratum(faultDMMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(faultDMMesh, &subpointMap);CHECK_PETSC_ERROR(err);
- CPPUNIT_ASSERT(subpointMap);
- err = ISGetSize(subpointMap, &numPoints);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(faultDMMesh, &subpointIS);CHECK_PETSC_ERROR(err);
+ CPPUNIT_ASSERT(subpointIS);
+ err = ISGetSize(subpointIS, &numPoints);CHECK_PETSC_ERROR(err);
PetscSection dispSection = fields.get("disp(t)").petscSection();
Vec dispVec = fields.get("disp(t)").localVector();
@@ -613,7 +614,7 @@
int iVertex = 0;
const PylithScalar tolerance = 1.0e-06;
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(tractionsVec, &tractionsArray);CHECK_PETSC_ERROR(err);
err = VecGetArray(dispVec, &dispArray);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v, ++iVertex) {
@@ -642,7 +643,8 @@
} // for
err = VecRestoreArray(tractionsVec, &tractionsArray);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(dispVec, &dispArray);CHECK_PETSC_ERROR(err);
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
} // testCalcTractions
// ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveImpulses.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveImpulses.cc 2013-02-05 13:15:19 UTC (rev 21334)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveImpulses.cc 2013-02-05 13:15:49 UTC (rev 21335)
@@ -154,17 +154,17 @@
_initialize(&mesh, &fault, &fields);
DM dmMesh = fault._faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscInt vStart, vEnd, numPoints;
PetscErrorCode err;
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
- CPPUNIT_ASSERT(subpointMap);
- err = ISGetSize(subpointMap, &numPoints);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
+ CPPUNIT_ASSERT(subpointIS);
+ err = ISGetSize(subpointIS, &numPoints);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt faultPoint;
@@ -172,7 +172,8 @@
CPPUNIT_ASSERT(faultPoint >= 0);
CPPUNIT_ASSERT_EQUAL(faultPoint, v);
} // for
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
CPPUNIT_ASSERT_EQUAL(_data->numConstraintVert, vEnd-vStart);
// Check orientation
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc 2013-02-05 13:15:19 UTC (rev 21334)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc 2013-02-05 13:15:49 UTC (rev 21335)
@@ -140,17 +140,17 @@
_initialize(&mesh, &fault, &fields);
DM dmMesh = fault._faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscInt vStart, vEnd, numPoints;
PetscErrorCode err;
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
- CPPUNIT_ASSERT(subpointMap);
- err = ISGetSize(subpointMap, &numPoints);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
+ CPPUNIT_ASSERT(subpointIS);
+ err = ISGetSize(subpointIS, &numPoints);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt faultPoint;
@@ -158,7 +158,8 @@
CPPUNIT_ASSERT(faultPoint >= 0);
CPPUNIT_ASSERT_EQUAL(faultPoint, _data->verticesFault[v-vStart]);
} // for
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
CPPUNIT_ASSERT_EQUAL(_data->numFaultVertices, vEnd-vStart);
// Check cohesive vertex info
@@ -479,19 +480,19 @@
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
DM faultDMMesh = fault._faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
CPPUNIT_ASSERT(faultDMMesh);
- err = DMPlexGetSubpointMap(faultDMMesh, &subpointMap);CHECK_PETSC_ERROR(err);
- CPPUNIT_ASSERT(subpointMap);
+ err = DMPlexCreateSubpointIS(faultDMMesh, &subpointIS);CHECK_PETSC_ERROR(err);
+ CPPUNIT_ASSERT(subpointIS);
const PylithScalar tolerance = 1.0e-06;
int iVertex = 0;
const PetscInt *points;
PetscInt numPoints;
- err = ISGetSize(subpointMap, &numPoints);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetSize(subpointIS, &numPoints);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(jacobianVec, &jacobianArray);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v, ++iVertex) {
PetscInt faultPoint;
@@ -519,7 +520,8 @@
CPPUNIT_ASSERT_DOUBLES_EQUAL(valE, jacobianArray[off+d], tolerance);
} // for
} // for
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(jacobianVec, &jacobianArray);CHECK_PETSC_ERROR(err);
} // testIntegrateJacobianLumped
@@ -712,16 +714,16 @@
int iVertex = 0;
const PylithScalar tolerance = 1.0e-06;
DM faultDMMesh = fault._faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscInt numPoints, vStart, vEnd;
CPPUNIT_ASSERT(faultDMMesh);
err = DMPlexGetDepthStratum(faultDMMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(faultDMMesh, &subpointMap);CHECK_PETSC_ERROR(err);
- CPPUNIT_ASSERT(subpointMap);
- err = ISGetSize(subpointMap, &numPoints);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(faultDMMesh, &subpointIS);CHECK_PETSC_ERROR(err);
+ CPPUNIT_ASSERT(subpointIS);
+ err = ISGetSize(subpointIS, &numPoints);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(tractionsVec, &tractionsArray);CHECK_PETSC_ERROR(err);
err = VecGetArray(dispVec, &dispArray);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v, ++iVertex) {
@@ -746,7 +748,8 @@
CPPUNIT_ASSERT_DOUBLES_EQUAL(tractionE, tractionsArray[off+d], tolerance);
} // for
} // for
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(tractionsVec, &tractionsArray);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(dispVec, &dispArray);CHECK_PETSC_ERROR(err);
} // testCalcTractionsChange
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestLiuCosSlipFn.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestLiuCosSlipFn.cc 2013-02-05 13:15:19 UTC (rev 21334)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestLiuCosSlipFn.cc 2013-02-05 13:15:49 UTC (rev 21335)
@@ -360,7 +360,7 @@
const ALE::Obj<RealSection>& oldCoordSection = sieveMesh->getRealSection("coordinates");
faultSieveMesh->setRealSection("coordinates", oldCoordSection);
DM dmMesh = faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscSection coordSection;
PetscInt vStart, vEnd;
@@ -368,7 +368,7 @@
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
err = DMPlexGetCoordinateSection(dmMesh, &coordSection);CHECK_PETSC_ERROR(err);
err = PetscSectionSetChart(coordSection, vStart, vEnd);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
@@ -383,7 +383,7 @@
err = VecCreate(mesh->comm(), &coordVec);CHECK_PETSC_ERROR(err);
err = VecSetSizes(coordVec, coordSize, PETSC_DETERMINE);CHECK_PETSC_ERROR(err);
err = VecSetFromOptions(coordVec);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt off;
@@ -394,7 +394,8 @@
coords[off+d] = oldCoords[d];
}
}
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
err = DMSetCoordinatesLocal(dmMesh, coordVec);CHECK_PETSC_ERROR(err);
@@ -469,7 +470,7 @@
const ALE::Obj<RealSection>& oldCoordSection = sieveMesh->getRealSection("coordinates");
faultSieveMesh->setRealSection("coordinates", oldCoordSection);
DM dmMesh = faultMesh.dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscSection coordSection;
PetscInt vStart, vEnd;
@@ -477,7 +478,7 @@
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
err = DMPlexGetCoordinateSection(dmMesh, &coordSection);CHECK_PETSC_ERROR(err);
err = PetscSectionSetChart(coordSection, vStart, vEnd);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
@@ -492,7 +493,7 @@
err = VecCreate(mesh.comm(), &coordVec);CHECK_PETSC_ERROR(err);
err = VecSetSizes(coordVec, coordSize, PETSC_DETERMINE);CHECK_PETSC_ERROR(err);
err = VecSetFromOptions(coordVec);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt off;
@@ -503,7 +504,8 @@
coords[off+d] = oldCoords[d];
}
}
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
err = DMSetCoordinatesLocal(dmMesh, coordVec);CHECK_PETSC_ERROR(err);
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestStepSlipFn.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestStepSlipFn.cc 2013-02-05 13:15:19 UTC (rev 21334)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestStepSlipFn.cc 2013-02-05 13:15:49 UTC (rev 21335)
@@ -293,7 +293,7 @@
const ALE::Obj<RealSection>& oldCoordSection = sieveMesh->getRealSection("coordinates");
faultSieveMesh->setRealSection("coordinates", oldCoordSection);
DM dmMesh = faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscSection coordSection;
PetscInt vStart, vEnd;
@@ -301,7 +301,7 @@
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
err = DMPlexGetCoordinateSection(dmMesh, &coordSection);CHECK_PETSC_ERROR(err);
err = PetscSectionSetChart(coordSection, vStart, vEnd);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
@@ -316,7 +316,7 @@
err = VecCreate(mesh->comm(), &coordVec);CHECK_PETSC_ERROR(err);
err = VecSetSizes(coordVec, coordSize, PETSC_DETERMINE);CHECK_PETSC_ERROR(err);
err = VecSetFromOptions(coordVec);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt off;
@@ -327,7 +327,8 @@
coords[off+d] = oldCoords[d];
}
}
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
err = DMSetCoordinatesLocal(dmMesh, coordVec);CHECK_PETSC_ERROR(err);
@@ -397,7 +398,7 @@
const ALE::Obj<RealSection>& oldCoordSection = sieveMesh->getRealSection("coordinates");
faultSieveMesh->setRealSection("coordinates", oldCoordSection);
DM dmMesh = faultMesh.dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscSection coordSection;
PetscInt vStart, vEnd;
@@ -405,7 +406,7 @@
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
err = DMPlexGetCoordinateSection(dmMesh, &coordSection);CHECK_PETSC_ERROR(err);
err = PetscSectionSetChart(coordSection, vStart, vEnd);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
@@ -420,7 +421,7 @@
err = VecCreate(mesh.comm(), &coordVec);CHECK_PETSC_ERROR(err);
err = VecSetSizes(coordVec, coordSize, PETSC_DETERMINE);CHECK_PETSC_ERROR(err);
err = VecSetFromOptions(coordVec);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt off;
@@ -431,7 +432,8 @@
coords[off+d] = oldCoords[d];
}
}
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
err = DMSetCoordinatesLocal(dmMesh, coordVec);CHECK_PETSC_ERROR(err);
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestTimeHistorySlipFn.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestTimeHistorySlipFn.cc 2013-02-05 13:15:19 UTC (rev 21334)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestTimeHistorySlipFn.cc 2013-02-05 13:15:49 UTC (rev 21335)
@@ -322,7 +322,7 @@
const ALE::Obj<RealSection>& oldCoordSection = sieveMesh->getRealSection("coordinates");
faultSieveMesh->setRealSection("coordinates", oldCoordSection);
DM dmMesh = faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscSection coordSection;
PetscInt vStart, vEnd;
@@ -330,7 +330,7 @@
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
err = DMPlexGetCoordinateSection(dmMesh, &coordSection);CHECK_PETSC_ERROR(err);
err = PetscSectionSetChart(coordSection, vStart, vEnd);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
@@ -345,7 +345,7 @@
err = VecCreate(mesh->comm(), &coordVec);CHECK_PETSC_ERROR(err);
err = VecSetSizes(coordVec, coordSize, PETSC_DETERMINE);CHECK_PETSC_ERROR(err);
err = VecSetFromOptions(coordVec);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt off;
@@ -356,7 +356,8 @@
coords[off+d] = oldCoords[d];
}
}
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
err = DMSetCoordinatesLocal(dmMesh, coordVec);CHECK_PETSC_ERROR(err);
@@ -429,7 +430,7 @@
const ALE::Obj<RealSection>& oldCoordSection = sieveMesh->getRealSection("coordinates");
faultSieveMesh->setRealSection("coordinates", oldCoordSection);
DM dmMesh = faultMesh.dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscSection coordSection;
PetscInt vStart, vEnd;
@@ -437,7 +438,7 @@
CPPUNIT_ASSERT(dmMesh);
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(dmMesh, &subpointMap);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(dmMesh, &subpointIS);CHECK_PETSC_ERROR(err);
err = DMPlexGetCoordinateSection(dmMesh, &coordSection);CHECK_PETSC_ERROR(err);
err = PetscSectionSetChart(coordSection, vStart, vEnd);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
@@ -452,7 +453,7 @@
err = VecCreate(mesh.comm(), &coordVec);CHECK_PETSC_ERROR(err);
err = VecSetSizes(coordVec, coordSize, PETSC_DETERMINE);CHECK_PETSC_ERROR(err);
err = VecSetFromOptions(coordVec);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt off;
@@ -463,7 +464,8 @@
coords[off+d] = oldCoords[d];
}
}
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
err = DMSetCoordinatesLocal(dmMesh, coordVec);CHECK_PETSC_ERROR(err);
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestTractPerturbation.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestTractPerturbation.cc 2013-02-05 13:15:19 UTC (rev 21334)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestTractPerturbation.cc 2013-02-05 13:15:49 UTC (rev 21335)
@@ -376,14 +376,14 @@
const ALE::Obj<RealSection>& oldCoordSection = sieveMesh->getRealSection("coordinates");
faultSieveMesh->setRealSection("coordinates", oldCoordSection);
DM faultDMMesh = faultMesh->dmMesh();
- IS subpointMap;
+ IS subpointIS;
const PetscInt *points;
PetscSection coordSection;
PetscInt vStart, vEnd;
CPPUNIT_ASSERT(faultDMMesh);
err = DMPlexGetDepthStratum(faultDMMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
- err = DMPlexGetSubpointMap(faultDMMesh, &subpointMap);CHECK_PETSC_ERROR(err);
+ err = DMPlexCreateSubpointIS(faultDMMesh, &subpointIS);CHECK_PETSC_ERROR(err);
err = DMPlexGetCoordinateSection(faultDMMesh, &coordSection);CHECK_PETSC_ERROR(err);
err = PetscSectionSetChart(coordSection, vStart, vEnd);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
@@ -398,7 +398,7 @@
err = VecCreate(mesh->comm(), &coordVec);CHECK_PETSC_ERROR(err);
err = VecSetSizes(coordVec, coordSize, PETSC_DETERMINE);CHECK_PETSC_ERROR(err);
err = VecSetFromOptions(coordVec);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
err = VecGetArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
for(PetscInt v = vStart; v < vEnd; ++v) {
PetscInt off;
@@ -409,7 +409,8 @@
coords[off+d] = oldCoords[d];
}
}
- err = ISRestoreIndices(subpointMap, &points);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &points);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
err = VecRestoreArray(coordVec, &coords);CHECK_PETSC_ERROR(err);
err = DMSetCoordinatesLocal(faultDMMesh, coordVec);CHECK_PETSC_ERROR(err);
Modified: short/3D/PyLith/trunk/unittests/libtests/meshio/TestDataWriterSubMesh.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/meshio/TestDataWriterSubMesh.cc 2013-02-05 13:15:19 UTC (rev 21334)
+++ short/3D/PyLith/trunk/unittests/libtests/meshio/TestDataWriterSubMesh.cc 2013-02-05 13:15:49 UTC (rev 21335)
@@ -120,10 +120,11 @@
err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, &vEnd);CHECK_PETSC_ERROR(err);
#if 0
DM dmSubmesh = _submesh->dmMesh();
- IS subpointMap;
- err = DMPlexGetSubpointMap(dmSubmesh, &subpointMap);CHECK_PETSC_ERROR(err);
- err = ISGetIndices(subpointMap, &ind);CHECK_PETSC_ERROR(err);
- err = ISRestoreIndices(subpointMap, &ind);CHECK_PETSC_ERROR(err);
+ IS subpointIS;
+ err = DMPlexCreateSubpointIS(dmSubmesh, &subpointIS);CHECK_PETSC_ERROR(err);
+ err = ISGetIndices(subpointIS, &ind);CHECK_PETSC_ERROR(err);
+ err = ISRestoreIndices(subpointIS, &ind);CHECK_PETSC_ERROR(err);
+ err = ISDestroy(&subpointIS);CHECK_PETSC_ERROR(err);
#endif
// Set vertex fields
More information about the CIG-COMMITS
mailing list