[cig-commits] r21921 - in short/3D/PyLith/trunk/libsrc/pylith: faults topology

[knepley at geodynamics.org]

Author: knepley
Date: 2013-04-23 07:13:40 -0700 (Tue, 23 Apr 2013)
New Revision: 21921

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc
   short/3D/PyLith/trunk/libsrc/pylith/topology/SubMesh.cc
Log:
Fixed DMPlexCreateSubmesh() calls

Modified: short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc	2013-04-23 00:40:28 UTC (rev 21920)
+++ short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc	2013-04-23 14:13:40 UTC (rev 21921)
@@ -77,11 +77,11 @@
     std::ostringstream tmp;
 
     err = DMLabelGetName(groupField, &groupName);CHECK_PETSC_ERROR(err);
-    err = DMPlexCreateSubmesh(dmMesh, groupName, &subdm);CHECK_PETSC_ERROR(err);
+    err = DMPlexCreateSubmesh(dmMesh, groupName, 1, &subdm);CHECK_PETSC_ERROR(err);
     err = DMPlexCreateLabel(subdm, labelName);CHECK_PETSC_ERROR(err);
     err = DMPlexGetLabel(subdm, labelName, &label);CHECK_PETSC_ERROR(err);
     err = DMPlexMarkBoundaryFaces(subdm, label);CHECK_PETSC_ERROR(err);
-    err = DMPlexCreateSubmesh(subdm, labelName, &faultBoundaryDM);CHECK_PETSC_ERROR(err);
+    err = DMPlexCreateSubmesh(subdm, labelName, 1, &faultBoundaryDM);CHECK_PETSC_ERROR(err);
     faultMesh->setDMMesh(subdm);
   } else {
     // TODO: This leg will be unnecessary
@@ -191,7 +191,7 @@
     err = DMPlexCreateLabel(faultDMMesh, labelName);CHECK_PETSC_ERROR(err);
     err = DMPlexGetLabel(faultDMMesh, labelName, &label);CHECK_PETSC_ERROR(err);
     err = DMPlexMarkBoundaryFaces(faultDMMesh, label);CHECK_PETSC_ERROR(err);
-    err = DMPlexCreateSubmesh(faultDMMesh, labelName, &faultBoundaryDM);CHECK_PETSC_ERROR(err);
+    err = DMPlexCreateSubmesh(faultDMMesh, labelName, 1, &faultBoundaryDM);CHECK_PETSC_ERROR(err);
 
     err = ISRestoreIndices(pointIS, &points);CHECK_PETSC_ERROR(err);
     err = ISDestroy(&pointIS);CHECK_PETSC_ERROR(err);

Modified: short/3D/PyLith/trunk/libsrc/pylith/topology/SubMesh.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/topology/SubMesh.cc	2013-04-23 00:40:28 UTC (rev 21920)
+++ short/3D/PyLith/trunk/libsrc/pylith/topology/SubMesh.cc	2013-04-23 14:13:40 UTC (rev 21921)
@@ -1,3 +1,4 @@
+
 // -*- C++ -*-
 //
 // ======================================================================
@@ -166,7 +167,7 @@
 
   /* TODO: Add creation of pointSF for submesh */
   err = DMDestroy(&_newMesh);CHECK_PETSC_ERROR(err);
-  err = DMPlexCreateSubmesh(dmMesh, label, &_newMesh);CHECK_PETSC_ERROR(err);
+  err = DMPlexCreateSubmesh(dmMesh, label, 1, &_newMesh);CHECK_PETSC_ERROR(err);
 
   // Set data from mesh.
   coordsys(mesh);