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

[knepley at geodynamics.org]

Author: knepley
Date: 2013-05-09 18:56:48 -0700 (Thu, 09 May 2013)
New Revision: 22022

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc
Log:
Fixed fault coordinate problems

Modified: short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc	2013-05-10 00:51:32 UTC (rev 22021)
+++ short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc	2013-05-10 01:56:48 UTC (rev 22022)
@@ -914,6 +914,7 @@
     PetscInt dof;
     err = PetscSectionGetDof(coordSection, v, &dof);PYLITH_CHECK_ERROR(err);
     err = PetscSectionSetDof(newCoordSection, v+extraCells, dof);PYLITH_CHECK_ERROR(err);
+    err = PetscSectionSetFieldDof(newCoordSection, v+extraCells, 0, dof);PYLITH_CHECK_ERROR(err);
   }
 
   if (debug) coordinates->view("Coordinates without shadow vertices");
@@ -924,7 +925,11 @@
 
     err = PetscSectionGetDof(coordSection, v, &dof);PYLITH_CHECK_ERROR(err);
     err = PetscSectionSetDof(newCoordSection, vertexRenumberDM[vnew], dof);PYLITH_CHECK_ERROR(err);
-    if (constraintCell) {err = PetscSectionSetDof(newCoordSection, vertexLagrangeRenumberDM[vnew], dof);PYLITH_CHECK_ERROR(err);}
+    err = PetscSectionSetFieldDof(newCoordSection, vertexRenumberDM[vnew], 0, dof);PYLITH_CHECK_ERROR(err);
+    if (constraintCell) {
+      err = PetscSectionSetDof(newCoordSection, vertexLagrangeRenumberDM[vnew], dof);PYLITH_CHECK_ERROR(err);
+      err = PetscSectionSetFieldDof(newCoordSection, vertexLagrangeRenumberDM[vnew], 0, dof);PYLITH_CHECK_ERROR(err);
+    }
 
     coordinates->addPoint(vertexRenumber[v], coordinates->getFiberDimension(v));
     if (constraintCell) coordinates->addPoint(vertexLagrangeRenumber[v], coordinates->getFiberDimension(v));