[cig-commits] [commit] knepley/upgrade-petsc-interface: Removed flip fault flag in creating fault mesh. (d60fd06)

cig_noreply at geodynamics.org cig_noreply at geodynamics.org
Tue Nov 5 15:54:36 PST 2013


Repository : ssh://geoshell/pylith

On branch  : knepley/upgrade-petsc-interface
Link       : https://github.com/geodynamics/pylith/compare/d80fbaf3bf5f3a2bcf921fa4fccec1f4f10cb3dd...f4fc4fbe8988e18f1e97a95b95ecae5a8ad9b7e9

>---------------------------------------------------------------

commit d60fd068198db06dd4d46a320cfb0f5df891975a
Author: Brad Aagaard <baagaard at usgs.gov>
Date:   Tue Nov 5 15:53:01 2013 -0800

    Removed flip fault flag in creating fault mesh.


>---------------------------------------------------------------

d60fd068198db06dd4d46a320cfb0f5df891975a
 libsrc/pylith/faults/CohesiveTopology.cc           | 70 +---------------------
 libsrc/pylith/faults/CohesiveTopology.hh           |  5 +-
 libsrc/pylith/faults/Fault.hh                      |  6 +-
 libsrc/pylith/faults/FaultCohesive.cc              |  5 +-
 libsrc/pylith/faults/FaultCohesive.hh              |  4 +-
 modulesrc/faults/Fault.i                           |  3 +-
 modulesrc/faults/FaultCohesive.i                   |  4 +-
 unittests/libtests/faults/TestFaultCohesiveDyn.cc  |  3 +-
 unittests/libtests/faults/TestFaultCohesiveDyn.hh  |  1 -
 .../libtests/faults/TestFaultCohesiveDynCases.cc   | 10 ----
 .../libtests/faults/TestFaultCohesiveImpulses.cc   |  3 +-
 .../libtests/faults/TestFaultCohesiveImpulses.hh   |  1 -
 .../faults/TestFaultCohesiveImpulsesCases.cc       |  4 --
 unittests/libtests/faults/TestFaultCohesiveKin.cc  |  4 +-
 unittests/libtests/faults/TestFaultCohesiveKin.hh  |  1 -
 .../libtests/faults/TestFaultCohesiveKinCases.cc   |  8 ---
 .../libtests/faults/TestFaultCohesiveKinSrcs.cc    |  2 -
 .../faults/TestFaultCohesiveKinSrcsCases.cc        |  4 --
 unittests/libtests/friction/TestFrictionModel.cc   |  3 +-
 unittests/libtests/meshio/TestDataWriterBCMesh.cc  |  3 +-
 unittests/libtests/meshio/TestDataWriterBCMesh.hh  |  1 -
 .../libtests/meshio/TestDataWriterFaultMesh.cc     |  3 +-
 .../libtests/meshio/TestDataWriterFaultMesh.hh     |  1 -
 .../meshio/TestDataWriterHDF5BCMeshCases.cc        |  4 --
 .../meshio/TestDataWriterHDF5ExtBCMeshCases.cc     |  4 --
 .../meshio/TestDataWriterHDF5ExtFaultMeshCases.cc  |  4 --
 .../meshio/TestDataWriterHDF5ExtMatMeshCases.cc    |  4 --
 .../meshio/TestDataWriterHDF5ExtMeshCases.cc       |  4 --
 .../meshio/TestDataWriterHDF5ExtSubMeshCases.cc    |  4 --
 .../meshio/TestDataWriterHDF5FaultMeshCases.cc     |  4 --
 .../meshio/TestDataWriterHDF5MatMeshCases.cc       |  4 --
 .../libtests/meshio/TestDataWriterHDF5MeshCases.cc |  4 --
 .../meshio/TestDataWriterHDF5SubMeshCases.cc       |  4 --
 unittests/libtests/meshio/TestDataWriterMesh.cc    |  3 +-
 unittests/libtests/meshio/TestDataWriterMesh.hh    |  3 +-
 unittests/libtests/meshio/TestDataWriterPoints.cc  |  3 +-
 unittests/libtests/meshio/TestDataWriterPoints.hh  |  1 -
 unittests/libtests/meshio/TestDataWriterSubMesh.cc |  3 +-
 unittests/libtests/meshio/TestDataWriterSubMesh.hh |  1 -
 .../meshio/TestDataWriterVTKBCMeshCases.cc         |  4 --
 .../meshio/TestDataWriterVTKFaultMeshCases.cc      |  4 --
 .../meshio/TestDataWriterVTKMatMeshCases.cc        |  4 --
 .../libtests/meshio/TestDataWriterVTKMeshCases.cc  |  4 --
 .../meshio/TestDataWriterVTKPointsCases.cc         |  4 --
 .../meshio/TestDataWriterVTKSubMeshCases.cc        |  4 --
 45 files changed, 21 insertions(+), 206 deletions(-)

diff --git a/libsrc/pylith/faults/CohesiveTopology.cc b/libsrc/pylith/faults/CohesiveTopology.cc
index 5d7a105..3791bc2 100644
--- a/libsrc/pylith/faults/CohesiveTopology.cc
+++ b/libsrc/pylith/faults/CohesiveTopology.cc
@@ -32,8 +32,7 @@ void
 pylith::faults::CohesiveTopology::createFault(topology::Mesh* faultMesh,
                                               PetscDM& faultBoundary,
                                               const topology::Mesh& mesh,
-                                              PetscDMLabel groupField,
-                                              const bool flipFault)
+                                              PetscDMLabel groupField)
 { // createFault
   PYLITH_METHOD_BEGIN;
 
@@ -57,47 +56,6 @@ pylith::faults::CohesiveTopology::createFault(topology::Mesh* faultMesh,
     err = DMPlexCreateSubmesh(dmMesh, groupField, 1, &subdm);PYLITH_CHECK_ERROR(err);
     err = DMPlexOrient(subdm);PYLITH_CHECK_ERROR(err);
 
-    if (flipFault) {
-      PetscInt maxConeSize, *revcone, *revconeO;
-      PetscInt pStart, pEnd, h;
-
-      err = DMPlexGetVTKCellHeight(subdm, &h);PYLITH_CHECK_ERROR(err);
-      err = DMPlexGetHeightStratum(subdm, h, &pStart, &pEnd);PYLITH_CHECK_ERROR(err);
-      err = DMPlexGetMaxSizes(subdm, &maxConeSize, NULL);PYLITH_CHECK_ERROR(err);
-      err = DMGetWorkArray(subdm, maxConeSize, PETSC_INT, &revcone);PYLITH_CHECK_ERROR(err);
-      err = DMGetWorkArray(subdm, maxConeSize, PETSC_INT, &revconeO);PYLITH_CHECK_ERROR(err);
-      for (PetscInt p = pStart; p < pEnd; ++p) {
-        const PetscInt *cone, *coneO, *support;
-        PetscInt        coneSize, faceSize, supportSize, cp, sp;
-
-        err = DMPlexGetConeSize(subdm, p, &coneSize);PYLITH_CHECK_ERROR(err);
-        err = DMPlexGetCone(subdm, p, &cone);PYLITH_CHECK_ERROR(err);
-        err = DMPlexGetConeOrientation(subdm, p, &coneO);PYLITH_CHECK_ERROR(err);
-        for (cp = 0; cp < coneSize; ++cp) {
-          err = DMPlexGetConeSize(subdm, cone[coneSize-1-cp], &faceSize);PYLITH_CHECK_ERROR(err);
-          revcone[cp]  = cone[coneSize-1-cp];
-          revconeO[cp] = coneO[coneSize-1-cp] >= 0 ? -(faceSize-coneO[coneSize-1-cp]) : faceSize+coneO[coneSize-1-cp];
-        }
-        err = DMPlexSetCone(subdm, p, revcone);PYLITH_CHECK_ERROR(err);
-        err = DMPlexSetConeOrientation(subdm, p, revconeO);PYLITH_CHECK_ERROR(err);
-        /* Reverse orientations of support */
-        faceSize = coneSize;
-        err = DMPlexGetSupportSize(subdm, p, &supportSize);PYLITH_CHECK_ERROR(err);
-        err = DMPlexGetSupport(subdm, p, &support);PYLITH_CHECK_ERROR(err);
-        for (sp = 0; sp < supportSize; ++sp) {
-          err = DMPlexGetConeSize(subdm, support[sp], &coneSize);PYLITH_CHECK_ERROR(err);
-          err = DMPlexGetCone(subdm, support[sp], &cone);PYLITH_CHECK_ERROR(err);
-          err = DMPlexGetConeOrientation(subdm, support[sp], &coneO);PYLITH_CHECK_ERROR(err);
-          for (cp = 0; cp < coneSize; ++cp) {
-            if (cone[cp] != p) continue;
-            err = DMPlexInsertConeOrientation(subdm, support[sp], cp, coneO[cp] >= 0 ? -(faceSize-coneO[cp]) : faceSize+coneO[cp]);PYLITH_CHECK_ERROR(err);
-          }
-        }
-      }
-      err = DMRestoreWorkArray(subdm, maxConeSize, PETSC_INT, &revcone);PYLITH_CHECK_ERROR(err);
-      err = DMRestoreWorkArray(subdm, maxConeSize, PETSC_INT, &revconeO);PYLITH_CHECK_ERROR(err);
-    }
-
     err = DMPlexCreateLabel(subdm, labelName);PYLITH_CHECK_ERROR(err);
     err = DMPlexGetLabel(subdm, labelName, &label);PYLITH_CHECK_ERROR(err);
     err = DMPlexMarkBoundaryFaces(subdm, label);PYLITH_CHECK_ERROR(err);
@@ -146,32 +104,6 @@ pylith::faults::CohesiveTopology::createFault(topology::Mesh* faultMesh,
     err = DMPlexSetSubpointMap(faultDMMesh, subpointMap);PYLITH_CHECK_ERROR(err);
     err = DMLabelDestroy(&subpointMap);PYLITH_CHECK_ERROR(err);
     err = DMDestroy(&faultDMMeshTmp);PYLITH_CHECK_ERROR(err);
-    if (flipFault) {
-      PetscInt maxConeSize, *revcone, *revconeO;
-      PetscInt pStart, pEnd;
-
-      err = DMPlexGetHeightStratum(faultDMMesh, h, &pStart, &pEnd);PYLITH_CHECK_ERROR(err);
-      err = DMPlexGetMaxSizes(faultDMMesh, &maxConeSize, NULL);PYLITH_CHECK_ERROR(err);
-      err = DMGetWorkArray(faultDMMesh, maxConeSize, PETSC_INT, &revcone);PYLITH_CHECK_ERROR(err);
-      err = DMGetWorkArray(faultDMMesh, maxConeSize, PETSC_INT, &revconeO);PYLITH_CHECK_ERROR(err);
-      for (PetscInt p = pStart; p < pEnd; ++p) {
-        const PetscInt *cone, *coneO;
-        PetscInt        coneSize, faceSize, c;
-
-        err = DMPlexGetConeSize(faultDMMesh, p, &coneSize);PYLITH_CHECK_ERROR(err);
-        err = DMPlexGetCone(faultDMMesh, p, &cone);PYLITH_CHECK_ERROR(err);
-        err = DMPlexGetConeOrientation(faultDMMesh, p, &coneO);PYLITH_CHECK_ERROR(err);
-        for (c = 0; c < coneSize; ++c) {
-          err = DMPlexGetConeSize(faultDMMesh, cone[coneSize-1-c], &faceSize);PYLITH_CHECK_ERROR(err);
-          revcone[c]  = cone[coneSize-1-c];
-          revconeO[c] = coneO[coneSize-1-c] >= 0 ? -(faceSize-coneO[coneSize-1-c]) : faceSize+coneO[coneSize-1-c];
-        }
-        err = DMPlexSetCone(faultDMMesh, p, revcone);PYLITH_CHECK_ERROR(err);
-        err = DMPlexSetConeOrientation(faultDMMesh, p, revconeO);PYLITH_CHECK_ERROR(err);
-      }
-      err = DMRestoreWorkArray(faultDMMesh, maxConeSize, PETSC_INT, &revcone);PYLITH_CHECK_ERROR(err);
-      err = DMRestoreWorkArray(faultDMMesh, maxConeSize, PETSC_INT, &revconeO);PYLITH_CHECK_ERROR(err);
-    }
     
     std::string submeshLabel = "fault_" + std::string(groupName);
     faultMesh->dmMesh(faultDMMesh, submeshLabel.c_str());
diff --git a/libsrc/pylith/faults/CohesiveTopology.hh b/libsrc/pylith/faults/CohesiveTopology.hh
index 5efa6f2..6a84813 100644
--- a/libsrc/pylith/faults/CohesiveTopology.hh
+++ b/libsrc/pylith/faults/CohesiveTopology.hh
@@ -49,15 +49,12 @@ public :
    * @param mesh Finite-element mesh of domain.
    * @param groupdField Group of vertices assocated with faces of
    *   cells defining fault surface
-   * @param flipFault Flag indicating to flip positive/negative sides
-   * of the fault.
    */
   static
   void createFault(topology::Mesh* faultMesh,
 		   DM& faultBoundary,
 		   const topology::Mesh& mesh,
-		   DMLabel groupField,
-		   const bool flipFault =false);
+		   DMLabel groupField);
 
   /** Create cohesive cells.
    *
diff --git a/libsrc/pylith/faults/Fault.hh b/libsrc/pylith/faults/Fault.hh
index 35d0d7f..178292b 100644
--- a/libsrc/pylith/faults/Fault.hh
+++ b/libsrc/pylith/faults/Fault.hh
@@ -133,13 +133,15 @@ public :
   /** Adjust mesh topology for fault implementation.
    *
    * @param mesh PETSc mesh
+   * @param firstFaultVertex The first point eligible to become a new fault vertex
+   * @param firstLagrangeVertex The first point eligible to become a new Lagrange vertex
+   * @param firstFaultCell The first point eligible to become a new fault cell
    */
   virtual
   void adjustTopology(topology::Mesh* const mesh,
                       int *firstFaultVertex,
                       int *firstLagrangeVertex,
-                      int *firstFaultCell,
-                      const bool flipFault = false) = 0;
+                      int *firstFaultCell) = 0;
 
   /** Initialize fault. Determine orientation and setup boundary
    * condition parameters.
diff --git a/libsrc/pylith/faults/FaultCohesive.cc b/libsrc/pylith/faults/FaultCohesive.cc
index 31a73e1..77a8e16 100644
--- a/libsrc/pylith/faults/FaultCohesive.cc
+++ b/libsrc/pylith/faults/FaultCohesive.cc
@@ -109,8 +109,7 @@ void
 pylith::faults::FaultCohesive::adjustTopology(topology::Mesh* const mesh,
                                               int *firstFaultVertex,
                                               int *firstLagrangeVertex,
-                                              int *firstFaultCell,
-                                              const bool flipFault)
+                                              int *firstFaultCell)
 { // adjustTopology
   PYLITH_METHOD_BEGIN;
 
@@ -144,7 +143,7 @@ pylith::faults::FaultCohesive::adjustTopology(topology::Mesh* const mesh,
       err = DMPlexGetDimension(dmMesh, &dim);PYLITH_CHECK_ERROR(err);
       err = DMPlexGetDepth(dmMesh, &depth);PYLITH_CHECK_ERROR(err);
       err = DMPlexGetLabel(dmMesh, charlabel, &groupField);PYLITH_CHECK_ERROR(err);
-      CohesiveTopology::createFault(&faultMesh, faultBoundary, *mesh, groupField, flipFault);
+      CohesiveTopology::createFault(&faultMesh, faultBoundary, *mesh, groupField);
 
       if (dim > 1 && dim == depth) {
         CohesiveTopology::createInterpolated(mesh, faultMesh, faultBoundary, groupField, id(), *firstFaultVertex, *firstLagrangeVertex, *firstFaultCell, useLagrangeConstraints());
diff --git a/libsrc/pylith/faults/FaultCohesive.hh b/libsrc/pylith/faults/FaultCohesive.hh
index cbb55af..a7ab31e 100644
--- a/libsrc/pylith/faults/FaultCohesive.hh
+++ b/libsrc/pylith/faults/FaultCohesive.hh
@@ -83,13 +83,11 @@ public :
    * @param firstFaultVertex The first point eligible to become a new fault vertex
    * @param firstLagrangeVertex The first point eligible to become a new Lagrange vertex
    * @param firstFaultCell The first point eligible to become a new fault cell
-   * @param flipFault Flip fault orientation.
    */
   void adjustTopology(topology::Mesh* const mesh,
                       int *firstFaultVertex,
                       int *firstLagrangeVertex,
-                      int *firstFaultCell,
-                      const bool flipFault = false);
+                      int *firstFaultCell);
 
   /** Cohesive cells use Lagrange multiplier constraints?
    *
diff --git a/modulesrc/faults/Fault.i b/modulesrc/faults/Fault.i
index 6fbc03b..1b20a70 100644
--- a/modulesrc/faults/Fault.i
+++ b/modulesrc/faults/Fault.i
@@ -119,8 +119,7 @@ namespace pylith {
       void adjustTopology(pylith::topology::Mesh* const mesh,
                           int *firstFaultVertex,
                           int *firstLagrangeVertex,
-                          int *firstFaultCell,
-                          const bool flipFault = false) = 0;
+                          int *firstFaultCell) = 0;
       %clear int *firstFaultVertex, int *firstFaultCell;
       
       /** Initialize fault. Determine orientation and setup boundary
diff --git a/modulesrc/faults/FaultCohesive.i b/modulesrc/faults/FaultCohesive.i
index 74d4301..19d267c 100644
--- a/modulesrc/faults/FaultCohesive.i
+++ b/modulesrc/faults/FaultCohesive.i
@@ -64,14 +64,12 @@ namespace pylith {
       /** Adjust mesh topology for fault implementation.
        *
        * @param mesh PETSc mesh.
-       * @param flipFault Flip fault orientation.
        */
       %apply int *INOUT {int *firstFaultVertex, int *firstLagrangeVertex, int *firstFaultCell};
       void adjustTopology(pylith::topology::Mesh* const mesh,
                           int *firstFaultVertex,
                           int *firstLagrangeVertex,
-                          int *firstFaultCell,
-                          const bool flipFault = false);
+                          int *firstFaultCell);
       %clear int *firstFaultVertex, int *firstLagrangeVertex, int *firstFaultCell;
       
       /** Cohesive cells use Lagrange multiplier constraints?
diff --git a/unittests/libtests/faults/TestFaultCohesiveDyn.cc b/unittests/libtests/faults/TestFaultCohesiveDyn.cc
index 00027ff..6ca49ef 100644
--- a/unittests/libtests/faults/TestFaultCohesiveDyn.cc
+++ b/unittests/libtests/faults/TestFaultCohesiveDyn.cc
@@ -59,7 +59,6 @@ pylith::faults::TestFaultCohesiveDyn::setUp(void)
   _dbInitialTract = 0;
   _friction = 0;
   _dbFriction = 0;
-  _flipFault = false;
 
   PYLITH_METHOD_END;
 } // setUp
@@ -679,7 +678,7 @@ pylith::faults::TestFaultCohesiveDyn::_initialize(topology::Mesh* const mesh,
   fault->label(_data->label);
   fault->quadrature(_quadrature);
   
-  fault->adjustTopology(mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
+  fault->adjustTopology(mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell);
   
   const PylithScalar upDir[3] = { 0.0, 0.0, 1.0 };
   
diff --git a/unittests/libtests/faults/TestFaultCohesiveDyn.hh b/unittests/libtests/faults/TestFaultCohesiveDyn.hh
index cdc77c0..3f23a4f 100644
--- a/unittests/libtests/faults/TestFaultCohesiveDyn.hh
+++ b/unittests/libtests/faults/TestFaultCohesiveDyn.hh
@@ -75,7 +75,6 @@ protected:
   spatialdata::spatialdb::SpatialDB* _dbInitialTract; ///< Initial tractions.
   friction::FrictionModel* _friction; ///< Friction model
   spatialdata::spatialdb::SpatialDB* _dbFriction; ///< Friction parameters.
-  bool _flipFault; ///< If true, flip fault orientation.
 
   // PUBLIC METHODS /////////////////////////////////////////////////////
 public:
diff --git a/unittests/libtests/faults/TestFaultCohesiveDynCases.cc b/unittests/libtests/faults/TestFaultCohesiveDynCases.cc
index 7c675c1..5382942 100644
--- a/unittests/libtests/faults/TestFaultCohesiveDynCases.cc
+++ b/unittests/libtests/faults/TestFaultCohesiveDynCases.cc
@@ -43,8 +43,6 @@ pylith::faults::TestFaultCohesiveDynTri3::setUp(void)
   feassemble::GeometryLine2D geometry;
   _quadrature->refGeometry(&geometry);
   
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
@@ -66,8 +64,6 @@ pylith::faults::TestFaultCohesiveDynTri3d::setUp(void)
   feassemble::GeometryLine2D geometry;
   _quadrature->refGeometry(&geometry);
   
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
@@ -89,8 +85,6 @@ pylith::faults::TestFaultCohesiveDynQuad4::setUp(void)
   feassemble::GeometryLine2D geometry;
   _quadrature->refGeometry(&geometry);
   
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
@@ -112,8 +106,6 @@ pylith::faults::TestFaultCohesiveDynTet4::setUp(void)
   feassemble::GeometryTri3D geometry;
   _quadrature->refGeometry(&geometry);
   
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
@@ -135,8 +127,6 @@ pylith::faults::TestFaultCohesiveDynHex8::setUp(void)
   feassemble::GeometryQuad3D geometry;
   _quadrature->refGeometry(&geometry);
   
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
diff --git a/unittests/libtests/faults/TestFaultCohesiveImpulses.cc b/unittests/libtests/faults/TestFaultCohesiveImpulses.cc
index 3a8a74f..026e983 100644
--- a/unittests/libtests/faults/TestFaultCohesiveImpulses.cc
+++ b/unittests/libtests/faults/TestFaultCohesiveImpulses.cc
@@ -52,7 +52,6 @@ pylith::faults::TestFaultCohesiveImpulses::setUp(void)
   _data = 0;
   _quadrature = new feassemble::Quadrature();CPPUNIT_ASSERT(_quadrature);
   _dbImpulseAmp = 0;
-  _flipFault = false;
 
   PYLITH_METHOD_END;
 } // setUp
@@ -371,7 +370,7 @@ pylith::faults::TestFaultCohesiveImpulses::_initialize(topology::Mesh* const mes
   fault->label(_data->label);
   fault->quadrature(_quadrature);
   
-  fault->adjustTopology(mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
+  fault->adjustTopology(mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell);
   
   const PylithScalar upDir[3] = { 0.0, 0.0, 1.0 };
   
diff --git a/unittests/libtests/faults/TestFaultCohesiveImpulses.hh b/unittests/libtests/faults/TestFaultCohesiveImpulses.hh
index 16015c6..2382fde 100644
--- a/unittests/libtests/faults/TestFaultCohesiveImpulses.hh
+++ b/unittests/libtests/faults/TestFaultCohesiveImpulses.hh
@@ -69,7 +69,6 @@ protected:
   CohesiveImpulsesData* _data; ///< Data for testing
   feassemble::Quadrature* _quadrature; ///< Fault quadrature.
   spatialdata::spatialdb::SpatialDB* _dbImpulseAmp; ///< Initial tractions.
-  bool _flipFault; ///< If true, flip fault orientation.
 
   // PUBLIC METHODS /////////////////////////////////////////////////////
 public:
diff --git a/unittests/libtests/faults/TestFaultCohesiveImpulsesCases.cc b/unittests/libtests/faults/TestFaultCohesiveImpulsesCases.cc
index 211a149..06c74bf 100644
--- a/unittests/libtests/faults/TestFaultCohesiveImpulsesCases.cc
+++ b/unittests/libtests/faults/TestFaultCohesiveImpulsesCases.cc
@@ -45,8 +45,6 @@ pylith::faults::TestFaultCohesiveImpulsesTri3::setUp(void)
   feassemble::GeometryLine2D geometry;
   _quadrature->refGeometry(&geometry); 
 
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
@@ -110,8 +108,6 @@ pylith::faults::TestFaultCohesiveImpulsesHex8::setUp(void)
   feassemble::GeometryQuad3D geometry;
   _quadrature->refGeometry(&geometry); 
 
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
diff --git a/unittests/libtests/faults/TestFaultCohesiveKin.cc b/unittests/libtests/faults/TestFaultCohesiveKin.cc
index e51799d..3a313c7 100644
--- a/unittests/libtests/faults/TestFaultCohesiveKin.cc
+++ b/unittests/libtests/faults/TestFaultCohesiveKin.cc
@@ -61,8 +61,6 @@ pylith::faults::TestFaultCohesiveKin::setUp(void)
   _slipfns.resize(nsrcs);
   _slipfns[0] = new BruneSlipFn();
 
-  _flipFault = false;
-
   PYLITH_METHOD_END;
 } // setUp
 
@@ -715,7 +713,7 @@ pylith::faults::TestFaultCohesiveKin::_initialize(topology::Mesh* const mesh,
   PetscErrorCode err = DMPlexGetStratumSize(dmMesh, _data->label, 1, &firstLagrangeVertex);PYLITH_CHECK_ERROR(err);
   PetscInt firstFaultCell = firstLagrangeVertex + firstLagrangeVertex;
   
-  fault->adjustTopology(mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
+  fault->adjustTopology(mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell);
   
   const PylithScalar upDir[3] = { 0.0, 0.0, 1.0 };
   fault->normalizer(normalizer);
diff --git a/unittests/libtests/faults/TestFaultCohesiveKin.hh b/unittests/libtests/faults/TestFaultCohesiveKin.hh
index 4000235..98579c1 100644
--- a/unittests/libtests/faults/TestFaultCohesiveKin.hh
+++ b/unittests/libtests/faults/TestFaultCohesiveKin.hh
@@ -65,7 +65,6 @@ protected :
   feassemble::Quadrature* _quadrature; ///< Fault quadrature.a
   std::vector<EqKinSrc*> _eqsrcs; ///< Array of Kinematic earthquake sources.
   std::vector<BruneSlipFn*> _slipfns; ///< Slip time function.
-  bool _flipFault; ///< If true, flip fault orientation.
 
   // PUBLIC METHODS /////////////////////////////////////////////////////
 public :
diff --git a/unittests/libtests/faults/TestFaultCohesiveKinCases.cc b/unittests/libtests/faults/TestFaultCohesiveKinCases.cc
index f1714d5..15f3612 100644
--- a/unittests/libtests/faults/TestFaultCohesiveKinCases.cc
+++ b/unittests/libtests/faults/TestFaultCohesiveKinCases.cc
@@ -44,8 +44,6 @@ pylith::faults::TestFaultCohesiveKinTri3::setUp(void)
   feassemble::GeometryLine2D geometry;
   _quadrature->refGeometry(&geometry);
   
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
@@ -67,8 +65,6 @@ pylith::faults::TestFaultCohesiveKinTri3d::setUp(void)
   feassemble::GeometryLine2D geometry;
   _quadrature->refGeometry(&geometry);
   
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
@@ -111,8 +107,6 @@ pylith::faults::TestFaultCohesiveKinQuad4e::setUp(void)
   feassemble::GeometryLine2D geometry;
   _quadrature->refGeometry(&geometry);
   
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
@@ -197,8 +191,6 @@ pylith::faults::TestFaultCohesiveKinHex8::setUp(void)
   feassemble::GeometryQuad3D geometry;
   _quadrature->refGeometry(&geometry);
   
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
diff --git a/unittests/libtests/faults/TestFaultCohesiveKinSrcs.cc b/unittests/libtests/faults/TestFaultCohesiveKinSrcs.cc
index dc1ed80..3d94c7a 100644
--- a/unittests/libtests/faults/TestFaultCohesiveKinSrcs.cc
+++ b/unittests/libtests/faults/TestFaultCohesiveKinSrcs.cc
@@ -57,8 +57,6 @@ pylith::faults::TestFaultCohesiveKinSrcs::setUp(void)
   _slipfns[0] = new BruneSlipFn();
   _slipfns[1] = new BruneSlipFn();
 
-  _flipFault = false;
-
   PYLITH_METHOD_END;
 } // setUp
 
diff --git a/unittests/libtests/faults/TestFaultCohesiveKinSrcsCases.cc b/unittests/libtests/faults/TestFaultCohesiveKinSrcsCases.cc
index 919bcba..5821ccc 100644
--- a/unittests/libtests/faults/TestFaultCohesiveKinSrcsCases.cc
+++ b/unittests/libtests/faults/TestFaultCohesiveKinSrcsCases.cc
@@ -45,8 +45,6 @@ pylith::faults::TestFaultCohesiveKinSrcsTri3::setUp(void)
   feassemble::GeometryLine2D geometry;
   _quadrature->refGeometry(&geometry);
   
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
@@ -110,8 +108,6 @@ pylith::faults::TestFaultCohesiveKinSrcsHex8::setUp(void)
   feassemble::GeometryQuad3D geometry;
   _quadrature->refGeometry(&geometry);
   
-  _flipFault = true;
-
   PYLITH_METHOD_END;
 } // setUp
 
diff --git a/unittests/libtests/friction/TestFrictionModel.cc b/unittests/libtests/friction/TestFrictionModel.cc
index 751a5a0..919c8c7 100644
--- a/unittests/libtests/friction/TestFrictionModel.cc
+++ b/unittests/libtests/friction/TestFrictionModel.cc
@@ -924,7 +924,6 @@ pylith::friction::TestFrictionModel::_initialize(topology::Mesh* mesh,
 			quadWts, numQuadPts,
 			spaceDim);
 
-  const bool flipFault = false;
   const char* label = "fault";
 
   PetscInt labelSize;
@@ -938,7 +937,7 @@ pylith::friction::TestFrictionModel::_initialize(topology::Mesh* mesh,
   fault->id(100);
   fault->label(label);
   fault->quadrature(&quadrature);
-  fault->adjustTopology(mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, flipFault);
+  fault->adjustTopology(mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell);
   fault->normalizer(normalizer);
 
   spatialdata::spatialdb::SimpleDB db;
diff --git a/unittests/libtests/meshio/TestDataWriterBCMesh.cc b/unittests/libtests/meshio/TestDataWriterBCMesh.cc
index 3e692cd..019c0ee 100644
--- a/unittests/libtests/meshio/TestDataWriterBCMesh.cc
+++ b/unittests/libtests/meshio/TestDataWriterBCMesh.cc
@@ -47,7 +47,6 @@ pylith::meshio::TestDataWriterBCMesh::setUp(void)
   _data = 0;
   _mesh = 0;
   _submesh = 0;
-  _flipFault = false;
 
   PYLITH_METHOD_END;
 } // setUp
@@ -100,7 +99,7 @@ pylith::meshio::TestDataWriterBCMesh::_initialize(void)
     } // if
     fault.label(_data->faultLabel);
     fault.id(_data->faultId);
-    fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
+    fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell);
   } // if
 
   delete _submesh; _submesh = new topology::Mesh(*_mesh, _data->bcLabel);CPPUNIT_ASSERT(_submesh);
diff --git a/unittests/libtests/meshio/TestDataWriterBCMesh.hh b/unittests/libtests/meshio/TestDataWriterBCMesh.hh
index 3868d15..1dff799 100644
--- a/unittests/libtests/meshio/TestDataWriterBCMesh.hh
+++ b/unittests/libtests/meshio/TestDataWriterBCMesh.hh
@@ -74,7 +74,6 @@ protected :
   DataWriterData* _data; ///< Data for testing
   topology::Mesh* _mesh; ///< Mesh for domain
   topology::Mesh* _submesh; ///< Mesh for subdomain.
-  bool _flipFault; ///< If true, flip fault orientation.
 
 }; // class TestDataWriterBCMesh
 
diff --git a/unittests/libtests/meshio/TestDataWriterFaultMesh.cc b/unittests/libtests/meshio/TestDataWriterFaultMesh.cc
index d9cc5cf..a573690 100644
--- a/unittests/libtests/meshio/TestDataWriterFaultMesh.cc
+++ b/unittests/libtests/meshio/TestDataWriterFaultMesh.cc
@@ -51,7 +51,6 @@ pylith::meshio::TestDataWriterFaultMesh::setUp(void)
   _mesh = new topology::Mesh();
   const bool isSubMesh = true;
   _faultMesh = new topology::Mesh(isSubMesh);
-  _flipFault = false;
 
   PYLITH_METHOD_END;
 } // setUp
@@ -104,7 +103,7 @@ pylith::meshio::TestDataWriterFaultMesh::_initialize(void)
   } // if
   fault.label(_data->faultLabel);
   fault.id(_data->faultId);
-  fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
+  fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell);
   faults::CohesiveTopology::createFaultParallel(_faultMesh, *_mesh, _data->faultId, _data->faultLabel, useLagrangeConstraints);
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterFaultMesh.hh b/unittests/libtests/meshio/TestDataWriterFaultMesh.hh
index e5c2ca9..4429ea7 100644
--- a/unittests/libtests/meshio/TestDataWriterFaultMesh.hh
+++ b/unittests/libtests/meshio/TestDataWriterFaultMesh.hh
@@ -74,7 +74,6 @@ protected :
   DataWriterData* _data; ///< Data for testing
   topology::Mesh* _mesh; ///< Mesh for domain
   topology::Mesh* _faultMesh; ///< Fault mesh.
-  bool _flipFault; ///< If true, flip fault orientation.
 
 }; // class TestDataWriterFaultMesh
 
diff --git a/unittests/libtests/meshio/TestDataWriterHDF5BCMeshCases.cc b/unittests/libtests/meshio/TestDataWriterHDF5BCMeshCases.cc
index ca6566a..19cb05e 100644
--- a/unittests/libtests/meshio/TestDataWriterHDF5BCMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterHDF5BCMeshCases.cc
@@ -34,7 +34,6 @@ pylith::meshio::TestDataWriterHDF5BCMeshTri3::setUp(void)
 
   TestDataWriterHDF5BCMesh::setUp();
   _data = new DataWriterHDF5DataBCMeshTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -53,7 +52,6 @@ pylith::meshio::TestDataWriterHDF5BCMeshQuad4::setUp(void)
 
   TestDataWriterHDF5BCMesh::setUp();
   _data = new DataWriterHDF5DataBCMeshQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -72,7 +70,6 @@ pylith::meshio::TestDataWriterHDF5BCMeshTet4::setUp(void)
 
   TestDataWriterHDF5BCMesh::setUp();
   _data = new DataWriterHDF5DataBCMeshTet4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -91,7 +88,6 @@ pylith::meshio::TestDataWriterHDF5BCMeshHex8::setUp(void)
 
   TestDataWriterHDF5BCMesh::setUp();
   _data = new DataWriterHDF5DataBCMeshHex8;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterHDF5ExtBCMeshCases.cc b/unittests/libtests/meshio/TestDataWriterHDF5ExtBCMeshCases.cc
index 7725243..be4adf5 100644
--- a/unittests/libtests/meshio/TestDataWriterHDF5ExtBCMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterHDF5ExtBCMeshCases.cc
@@ -35,7 +35,6 @@ pylith::meshio::TestDataWriterHDF5ExtBCMeshTri3::setUp(void)
 
   TestDataWriterHDF5ExtBCMesh::setUp();
   _data = new DataWriterHDF5DataBCMeshTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -54,7 +53,6 @@ pylith::meshio::TestDataWriterHDF5ExtBCMeshQuad4::setUp(void)
 
   TestDataWriterHDF5ExtBCMesh::setUp();
   _data = new DataWriterHDF5DataBCMeshQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -73,7 +71,6 @@ pylith::meshio::TestDataWriterHDF5ExtBCMeshTet4::setUp(void)
 
   TestDataWriterHDF5ExtBCMesh::setUp();
   _data = new DataWriterHDF5DataBCMeshTet4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -92,7 +89,6 @@ pylith::meshio::TestDataWriterHDF5ExtBCMeshHex8::setUp(void)
 
   TestDataWriterHDF5ExtBCMesh::setUp();
   _data = new DataWriterHDF5DataBCMeshHex8;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterHDF5ExtFaultMeshCases.cc b/unittests/libtests/meshio/TestDataWriterHDF5ExtFaultMeshCases.cc
index ccd7a62..2e7bf77 100644
--- a/unittests/libtests/meshio/TestDataWriterHDF5ExtFaultMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterHDF5ExtFaultMeshCases.cc
@@ -30,7 +30,6 @@ pylith::meshio::TestDataWriterHDF5ExtFaultMeshTri3::setUp(void)
 { // setUp
   TestDataWriterHDF5ExtFaultMesh::setUp();
   _data = new DataWriterHDF5DataFaultMeshTri3;
-  _flipFault = true;
 
   _initialize();
 } // setUp
@@ -46,7 +45,6 @@ pylith::meshio::TestDataWriterHDF5ExtFaultMeshQuad4::setUp(void)
 { // setUp
   TestDataWriterHDF5ExtFaultMesh::setUp();
   _data = new DataWriterHDF5DataFaultMeshQuad4;
-  _flipFault = true;
 
   _initialize();
 } // setUp
@@ -62,7 +60,6 @@ pylith::meshio::TestDataWriterHDF5ExtFaultMeshTet4::setUp(void)
 { // setUp
   TestDataWriterHDF5ExtFaultMesh::setUp();
   _data = new DataWriterHDF5DataFaultMeshTet4;
-  _flipFault = false;
 
   _initialize();
 } // setUp
@@ -78,7 +75,6 @@ pylith::meshio::TestDataWriterHDF5ExtFaultMeshHex8::setUp(void)
 { // setUp
   TestDataWriterHDF5ExtFaultMesh::setUp();
   _data = new DataWriterHDF5DataFaultMeshHex8;
-  _flipFault = true;
 
   _initialize();
 } // setUp
diff --git a/unittests/libtests/meshio/TestDataWriterHDF5ExtMatMeshCases.cc b/unittests/libtests/meshio/TestDataWriterHDF5ExtMatMeshCases.cc
index ee2f784..f5f2142 100644
--- a/unittests/libtests/meshio/TestDataWriterHDF5ExtMatMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterHDF5ExtMatMeshCases.cc
@@ -43,7 +43,6 @@ pylith::meshio::TestDataWriterHDF5ExtMatMeshTri3::setUp(void)
 
   TestDataWriterHDF5ExtMesh::setUp();
   _data = new DataWriterHDF5DataMatMeshTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -59,7 +58,6 @@ pylith::meshio::TestDataWriterHDF5ExtMatMeshQuad4::setUp(void)
 
   TestDataWriterHDF5ExtMesh::setUp();
   _data = new DataWriterHDF5DataMatMeshQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -75,7 +73,6 @@ pylith::meshio::TestDataWriterHDF5ExtMatMeshTet4::setUp(void)
 
   TestDataWriterHDF5ExtMesh::setUp();
   _data = new DataWriterHDF5DataMatMeshTet4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -91,7 +88,6 @@ pylith::meshio::TestDataWriterHDF5ExtMatMeshHex8::setUp(void)
 
   TestDataWriterHDF5ExtMesh::setUp();
   _data = new DataWriterHDF5DataMatMeshHex8;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterHDF5ExtMeshCases.cc b/unittests/libtests/meshio/TestDataWriterHDF5ExtMeshCases.cc
index a23faf5..ec660a7 100644
--- a/unittests/libtests/meshio/TestDataWriterHDF5ExtMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterHDF5ExtMeshCases.cc
@@ -43,7 +43,6 @@ pylith::meshio::TestDataWriterHDF5ExtMeshTri3::setUp(void)
 
   TestDataWriterHDF5ExtMesh::setUp();
   _data = new DataWriterHDF5DataMeshTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -59,7 +58,6 @@ pylith::meshio::TestDataWriterHDF5ExtMeshQuad4::setUp(void)
 
   TestDataWriterHDF5ExtMesh::setUp();
   _data = new DataWriterHDF5DataMeshQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -75,7 +73,6 @@ pylith::meshio::TestDataWriterHDF5ExtMeshTet4::setUp(void)
 
   TestDataWriterHDF5ExtMesh::setUp();
   _data = new DataWriterHDF5DataMeshTet4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -91,7 +88,6 @@ pylith::meshio::TestDataWriterHDF5ExtMeshHex8::setUp(void)
 
   TestDataWriterHDF5ExtMesh::setUp();
   _data = new DataWriterHDF5DataMeshHex8;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterHDF5ExtSubMeshCases.cc b/unittests/libtests/meshio/TestDataWriterHDF5ExtSubMeshCases.cc
index db2cd90..5f98c6d 100644
--- a/unittests/libtests/meshio/TestDataWriterHDF5ExtSubMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterHDF5ExtSubMeshCases.cc
@@ -34,7 +34,6 @@ pylith::meshio::TestDataWriterHDF5ExtSubMeshTri3::setUp(void)
 
   TestDataWriterHDF5ExtSubMesh::setUp();
   _data = new DataWriterHDF5DataSubMeshTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -53,7 +52,6 @@ pylith::meshio::TestDataWriterHDF5ExtSubMeshQuad4::setUp(void)
 
   TestDataWriterHDF5ExtSubMesh::setUp();
   _data = new DataWriterHDF5DataSubMeshQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -72,7 +70,6 @@ pylith::meshio::TestDataWriterHDF5ExtSubMeshTet4::setUp(void)
 
   TestDataWriterHDF5ExtSubMesh::setUp();
   _data = new DataWriterHDF5DataSubMeshTet4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -91,7 +88,6 @@ pylith::meshio::TestDataWriterHDF5ExtSubMeshHex8::setUp(void)
 
   TestDataWriterHDF5ExtSubMesh::setUp();
   _data = new DataWriterHDF5DataSubMeshHex8;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterHDF5FaultMeshCases.cc b/unittests/libtests/meshio/TestDataWriterHDF5FaultMeshCases.cc
index 3b8b86f..2d62d05 100644
--- a/unittests/libtests/meshio/TestDataWriterHDF5FaultMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterHDF5FaultMeshCases.cc
@@ -35,7 +35,6 @@ pylith::meshio::TestDataWriterHDF5FaultMeshTri3::setUp(void)
 
   TestDataWriterHDF5FaultMesh::setUp();
   _data = new DataWriterHDF5DataFaultMeshTri3;
-  _flipFault = true;
 
   _initialize();
 
@@ -55,7 +54,6 @@ pylith::meshio::TestDataWriterHDF5FaultMeshQuad4::setUp(void)
 
   TestDataWriterHDF5FaultMesh::setUp();
   _data = new DataWriterHDF5DataFaultMeshQuad4;
-  _flipFault = true;
 
   _initialize();
 
@@ -75,7 +73,6 @@ pylith::meshio::TestDataWriterHDF5FaultMeshTet4::setUp(void)
 
   TestDataWriterHDF5FaultMesh::setUp();
   _data = new DataWriterHDF5DataFaultMeshTet4;
-  _flipFault = false;
 
   _initialize();
 
@@ -95,7 +92,6 @@ pylith::meshio::TestDataWriterHDF5FaultMeshHex8::setUp(void)
 
   TestDataWriterHDF5FaultMesh::setUp();
   _data = new DataWriterHDF5DataFaultMeshHex8;
-  _flipFault = true;
 
   _initialize();
 
diff --git a/unittests/libtests/meshio/TestDataWriterHDF5MatMeshCases.cc b/unittests/libtests/meshio/TestDataWriterHDF5MatMeshCases.cc
index db7f533..2bcb000 100644
--- a/unittests/libtests/meshio/TestDataWriterHDF5MatMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterHDF5MatMeshCases.cc
@@ -43,7 +43,6 @@ pylith::meshio::TestDataWriterHDF5MatMeshTri3::setUp(void)
 
   TestDataWriterHDF5Mesh::setUp();
   _data = new DataWriterHDF5DataMatMeshTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -59,7 +58,6 @@ pylith::meshio::TestDataWriterHDF5MatMeshQuad4::setUp(void)
 
   TestDataWriterHDF5Mesh::setUp();
   _data = new DataWriterHDF5DataMatMeshQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -75,7 +73,6 @@ pylith::meshio::TestDataWriterHDF5MatMeshTet4::setUp(void)
 
   TestDataWriterHDF5Mesh::setUp();
   _data = new DataWriterHDF5DataMatMeshTet4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -91,7 +88,6 @@ pylith::meshio::TestDataWriterHDF5MatMeshHex8::setUp(void)
 
   TestDataWriterHDF5Mesh::setUp();
   _data = new DataWriterHDF5DataMatMeshHex8;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterHDF5MeshCases.cc b/unittests/libtests/meshio/TestDataWriterHDF5MeshCases.cc
index 5c65349..35e5aa4 100644
--- a/unittests/libtests/meshio/TestDataWriterHDF5MeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterHDF5MeshCases.cc
@@ -43,7 +43,6 @@ pylith::meshio::TestDataWriterHDF5MeshTri3::setUp(void)
 
   TestDataWriterHDF5Mesh::setUp();
   _data = new DataWriterHDF5DataMeshTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -59,7 +58,6 @@ pylith::meshio::TestDataWriterHDF5MeshQuad4::setUp(void)
 
   TestDataWriterHDF5Mesh::setUp();
   _data = new DataWriterHDF5DataMeshQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -75,7 +73,6 @@ pylith::meshio::TestDataWriterHDF5MeshTet4::setUp(void)
 
   TestDataWriterHDF5Mesh::setUp();
   _data = new DataWriterHDF5DataMeshTet4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -91,7 +88,6 @@ pylith::meshio::TestDataWriterHDF5MeshHex8::setUp(void)
 
   TestDataWriterHDF5Mesh::setUp();
   _data = new DataWriterHDF5DataMeshHex8;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterHDF5SubMeshCases.cc b/unittests/libtests/meshio/TestDataWriterHDF5SubMeshCases.cc
index 137ddbe..9fe2de6 100644
--- a/unittests/libtests/meshio/TestDataWriterHDF5SubMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterHDF5SubMeshCases.cc
@@ -35,7 +35,6 @@ pylith::meshio::TestDataWriterHDF5SubMeshTri3::setUp(void)
 
   TestDataWriterHDF5SubMesh::setUp();
   _data = new DataWriterHDF5DataSubMeshTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -54,7 +53,6 @@ pylith::meshio::TestDataWriterHDF5SubMeshQuad4::setUp(void)
 
   TestDataWriterHDF5SubMesh::setUp();
   _data = new DataWriterHDF5DataSubMeshQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -73,7 +71,6 @@ pylith::meshio::TestDataWriterHDF5SubMeshTet4::setUp(void)
 
   TestDataWriterHDF5SubMesh::setUp();
   _data = new DataWriterHDF5DataSubMeshTet4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -92,7 +89,6 @@ pylith::meshio::TestDataWriterHDF5SubMeshHex8::setUp(void)
 
   TestDataWriterHDF5SubMesh::setUp();
   _data = new DataWriterHDF5DataSubMeshHex8;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterMesh.cc b/unittests/libtests/meshio/TestDataWriterMesh.cc
index 94ae302..2289814 100644
--- a/unittests/libtests/meshio/TestDataWriterMesh.cc
+++ b/unittests/libtests/meshio/TestDataWriterMesh.cc
@@ -46,7 +46,6 @@ pylith::meshio::TestDataWriterMesh::setUp(void)
 
   _data = 0;
   _mesh = 0;
-  _flipFault = false;
 
   PYLITH_METHOD_END;
 } // setUp
@@ -98,7 +97,7 @@ pylith::meshio::TestDataWriterMesh::_initialize(void)
     } // if
     fault.label(_data->faultLabel);
     fault.id(_data->faultId);
-    fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
+    fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell);
   } // if
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterMesh.hh b/unittests/libtests/meshio/TestDataWriterMesh.hh
index 7c1dae6..1c7d561 100644
--- a/unittests/libtests/meshio/TestDataWriterMesh.hh
+++ b/unittests/libtests/meshio/TestDataWriterMesh.hh
@@ -27,7 +27,7 @@
 #if !defined(pylith_meshio_testdatawritermesh_hh)
 #define pylith_meshio_testdatawritermesh_hh
 
-#include "pylith/topology/topologyfwd.hh" // USES Mesh, Field
+#include "pylith/topology/topologyfwd.hh" // HOLDSA Mesh, Field
 
 /// Namespace for pylith package
 namespace pylith {
@@ -73,7 +73,6 @@ protected :
 
   DataWriterData* _data; ///< Data for testing
   topology::Mesh* _mesh; ///< Mesh for data
-  bool _flipFault; ///< If true, flip fault orientation.
 
 }; // class TestDataWriterMesh
 
diff --git a/unittests/libtests/meshio/TestDataWriterPoints.cc b/unittests/libtests/meshio/TestDataWriterPoints.cc
index 7dc9d88..0146dd3 100644
--- a/unittests/libtests/meshio/TestDataWriterPoints.cc
+++ b/unittests/libtests/meshio/TestDataWriterPoints.cc
@@ -46,7 +46,6 @@ pylith::meshio::TestDataWriterPoints::setUp(void)
 
   _data = 0;
   _mesh = 0;
-  _flipFault = false;
 
   PYLITH_METHOD_END;
 } // setUp
@@ -98,7 +97,7 @@ pylith::meshio::TestDataWriterPoints::_initialize(void)
     } // if
     fault.label(_data->faultLabel);
     fault.id(_data->faultId);
-    fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
+    fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell);
   } // if
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterPoints.hh b/unittests/libtests/meshio/TestDataWriterPoints.hh
index b2f5b07..9552cad 100644
--- a/unittests/libtests/meshio/TestDataWriterPoints.hh
+++ b/unittests/libtests/meshio/TestDataWriterPoints.hh
@@ -66,7 +66,6 @@ protected :
 
   DataWriterDataPoints* _data; ///< Data for testing
   topology::Mesh* _mesh; ///< Mesh for data
-  bool _flipFault; ///< If true, flip fault orientation.
 
 }; // class TestDataWriterPoints
 
diff --git a/unittests/libtests/meshio/TestDataWriterSubMesh.cc b/unittests/libtests/meshio/TestDataWriterSubMesh.cc
index 86eb87c..fbeb7c2 100644
--- a/unittests/libtests/meshio/TestDataWriterSubMesh.cc
+++ b/unittests/libtests/meshio/TestDataWriterSubMesh.cc
@@ -47,7 +47,6 @@ pylith::meshio::TestDataWriterSubMesh::setUp(void)
   _data = 0;
   _mesh = 0;
   _submesh = 0;
-  _flipFault = false;
 
   PYLITH_METHOD_END;
 } // setUp
@@ -100,7 +99,7 @@ pylith::meshio::TestDataWriterSubMesh::_initialize(void)
     } // if
     fault.label(_data->faultLabel);
     fault.id(_data->faultId);
-    fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
+    fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell);
   } // if
 
   CPPUNIT_ASSERT(_data->bcLabel);
diff --git a/unittests/libtests/meshio/TestDataWriterSubMesh.hh b/unittests/libtests/meshio/TestDataWriterSubMesh.hh
index c3d5ad4..ca9d1e4 100644
--- a/unittests/libtests/meshio/TestDataWriterSubMesh.hh
+++ b/unittests/libtests/meshio/TestDataWriterSubMesh.hh
@@ -74,7 +74,6 @@ protected :
   DataWriterData* _data; ///< Data for testing
   topology::Mesh* _mesh; ///< Mesh for domain
   topology::Mesh* _submesh; ///< Mesh for subdomain.
-  bool _flipFault; ///< If true, flip fault orientation.
 
 }; // class TestDataWriterSubMesh
 
diff --git a/unittests/libtests/meshio/TestDataWriterVTKBCMeshCases.cc b/unittests/libtests/meshio/TestDataWriterVTKBCMeshCases.cc
index 7402953..c5ebf88 100644
--- a/unittests/libtests/meshio/TestDataWriterVTKBCMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterVTKBCMeshCases.cc
@@ -35,7 +35,6 @@ pylith::meshio::TestDataWriterVTKBCMeshTri3::setUp(void)
 
   TestDataWriterVTKBCMesh::setUp();
   _data = new DataWriterVTKDataBCMeshTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -54,7 +53,6 @@ pylith::meshio::TestDataWriterVTKBCMeshQuad4::setUp(void)
 
   TestDataWriterVTKBCMesh::setUp();
   _data = new DataWriterVTKDataBCMeshQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -73,7 +71,6 @@ pylith::meshio::TestDataWriterVTKBCMeshTet4::setUp(void)
 
   TestDataWriterVTKBCMesh::setUp();
   _data = new DataWriterVTKDataBCMeshTet4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -92,7 +89,6 @@ pylith::meshio::TestDataWriterVTKBCMeshHex8::setUp(void)
 
   TestDataWriterVTKBCMesh::setUp();
   _data = new DataWriterVTKDataBCMeshHex8;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterVTKFaultMeshCases.cc b/unittests/libtests/meshio/TestDataWriterVTKFaultMeshCases.cc
index fa7b63f..007e3b4 100644
--- a/unittests/libtests/meshio/TestDataWriterVTKFaultMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterVTKFaultMeshCases.cc
@@ -34,7 +34,6 @@ pylith::meshio::TestDataWriterVTKFaultMeshTri3::setUp(void)
 
   TestDataWriterVTKFaultMesh::setUp();
   _data = new DataWriterVTKDataFaultMeshTri3;
-  _flipFault = true;
 
   _initialize();
 
@@ -54,7 +53,6 @@ pylith::meshio::TestDataWriterVTKFaultMeshQuad4::setUp(void)
 
   TestDataWriterVTKFaultMesh::setUp();
   _data = new DataWriterVTKDataFaultMeshQuad4;
-  _flipFault = true;
 
   _initialize();
 
@@ -74,7 +72,6 @@ pylith::meshio::TestDataWriterVTKFaultMeshTet4::setUp(void)
 
   TestDataWriterVTKFaultMesh::setUp();
   _data = new DataWriterVTKDataFaultMeshTet4;
-  _flipFault = false;
 
   _initialize();
 
@@ -94,7 +91,6 @@ pylith::meshio::TestDataWriterVTKFaultMeshHex8::setUp(void)
 
   TestDataWriterVTKFaultMesh::setUp();
   _data = new DataWriterVTKDataFaultMeshHex8;
-  _flipFault = true;
 
   _initialize();
 
diff --git a/unittests/libtests/meshio/TestDataWriterVTKMatMeshCases.cc b/unittests/libtests/meshio/TestDataWriterVTKMatMeshCases.cc
index cde222c..eac8d71 100644
--- a/unittests/libtests/meshio/TestDataWriterVTKMatMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterVTKMatMeshCases.cc
@@ -44,7 +44,6 @@ pylith::meshio::TestDataWriterVTKMatMeshTri3::setUp(void)
 
   TestDataWriterVTKMesh::setUp();
   _data = new DataWriterVTKDataMatMeshTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -60,7 +59,6 @@ pylith::meshio::TestDataWriterVTKMatMeshQuad4::setUp(void)
 
   TestDataWriterVTKMesh::setUp();
   _data = new DataWriterVTKDataMatMeshQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -76,7 +74,6 @@ pylith::meshio::TestDataWriterVTKMatMeshTet4::setUp(void)
 
   TestDataWriterVTKMesh::setUp();
   _data = new DataWriterVTKDataMatMeshTet4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -92,7 +89,6 @@ pylith::meshio::TestDataWriterVTKMatMeshHex8::setUp(void)
 
   TestDataWriterVTKMesh::setUp();
   _data = new DataWriterVTKDataMatMeshHex8;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterVTKMeshCases.cc b/unittests/libtests/meshio/TestDataWriterVTKMeshCases.cc
index 2738e57..1c770f2 100644
--- a/unittests/libtests/meshio/TestDataWriterVTKMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterVTKMeshCases.cc
@@ -44,7 +44,6 @@ pylith::meshio::TestDataWriterVTKMeshTri3::setUp(void)
 
   TestDataWriterVTKMesh::setUp();
   _data = new DataWriterVTKDataMeshTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -60,7 +59,6 @@ pylith::meshio::TestDataWriterVTKMeshQuad4::setUp(void)
 
   TestDataWriterVTKMesh::setUp();
   _data = new DataWriterVTKDataMeshQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -76,7 +74,6 @@ pylith::meshio::TestDataWriterVTKMeshTet4::setUp(void)
 
   TestDataWriterVTKMesh::setUp();
   _data = new DataWriterVTKDataMeshTet4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -92,7 +89,6 @@ pylith::meshio::TestDataWriterVTKMeshHex8::setUp(void)
 
   TestDataWriterVTKMesh::setUp();
   _data = new DataWriterVTKDataMeshHex8;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterVTKPointsCases.cc b/unittests/libtests/meshio/TestDataWriterVTKPointsCases.cc
index 6bb0c44..d15ebae 100644
--- a/unittests/libtests/meshio/TestDataWriterVTKPointsCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterVTKPointsCases.cc
@@ -44,7 +44,6 @@ pylith::meshio::TestDataWriterVTKPointsTri3::setUp(void)
 
   TestDataWriterVTKPoints::setUp();
   _data = new DataWriterVTKDataPointsTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -60,7 +59,6 @@ pylith::meshio::TestDataWriterVTKPointsQuad4::setUp(void)
 
   TestDataWriterVTKPoints::setUp();
   _data = new DataWriterVTKDataPointsQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -76,7 +74,6 @@ pylith::meshio::TestDataWriterVTKPointsTet4::setUp(void)
 
   TestDataWriterVTKPoints::setUp();
   _data = new DataWriterVTKDataPointsTet4;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -92,7 +89,6 @@ pylith::meshio::TestDataWriterVTKPointsHex8::setUp(void)
 
   TestDataWriterVTKPoints::setUp();
   _data = new DataWriterVTKDataPointsHex8;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
diff --git a/unittests/libtests/meshio/TestDataWriterVTKSubMeshCases.cc b/unittests/libtests/meshio/TestDataWriterVTKSubMeshCases.cc
index bc6a5bd..6238451 100644
--- a/unittests/libtests/meshio/TestDataWriterVTKSubMeshCases.cc
+++ b/unittests/libtests/meshio/TestDataWriterVTKSubMeshCases.cc
@@ -35,7 +35,6 @@ pylith::meshio::TestDataWriterVTKSubMeshTri3::setUp(void)
 
   TestDataWriterVTKSubMesh::setUp();
   _data = new DataWriterVTKDataSubMeshTri3;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -54,7 +53,6 @@ pylith::meshio::TestDataWriterVTKSubMeshQuad4::setUp(void)
 
   TestDataWriterVTKSubMesh::setUp();
   _data = new DataWriterVTKDataSubMeshQuad4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -73,7 +71,6 @@ pylith::meshio::TestDataWriterVTKSubMeshTet4::setUp(void)
 
   TestDataWriterVTKSubMesh::setUp();
   _data = new DataWriterVTKDataSubMeshTet4;
-  _flipFault = false;
   _initialize();
 
   PYLITH_METHOD_END;
@@ -92,7 +89,6 @@ pylith::meshio::TestDataWriterVTKSubMeshHex8::setUp(void)
 
   TestDataWriterVTKSubMesh::setUp();
   _data = new DataWriterVTKDataSubMeshHex8;
-  _flipFault = true;
   _initialize();
 
   PYLITH_METHOD_END;



More information about the CIG-COMMITS mailing list