[cig-commits] r13873 - in short/3D/PyLith/trunk: libsrc/faults unittests/libtests/faults unittests/libtests/faults/data

brad at geodynamics.org brad at geodynamics.org
Thu Jan 15 14:51:53 PST 2009


Author: brad
Date: 2009-01-15 14:51:53 -0800 (Thu, 15 Jan 2009)
New Revision: 13873

Modified:
   short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc
   short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.hh
   short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.cc
   short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.hh
   short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc
   short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.hh
   short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinSrcs.cc
   short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinTri3.cc
   short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataTri3.cc
Log:
Worked on fixing flipFault switch (fixed creation of vertices but cohesive cells have wrong orientation). Looks like the flipFault switch won't work and need to change orientation of the fault face which is established in CohesiveTopology::createFaultSieveFromVertices().

Modified: short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc	2009-01-15 22:42:24 UTC (rev 13872)
+++ short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc	2009-01-15 22:51:53 UTC (rev 13873)
@@ -19,6 +19,7 @@
 
 #include <cassert> // USES assert()
 
+// ----------------------------------------------------------------------
 void
 pylith::faults::CohesiveTopology::createFaultSieveFromVertices(const int dim,
                                                                const int firstCell,
@@ -138,6 +139,7 @@
   if (!faultSieve->commRank()) delete [] indices;
 }
 
+// ----------------------------------------------------------------------
 void
 pylith::faults::CohesiveTopology::createFaultSieveFromFaces(const int dim,
                                                             const int firstCell,
@@ -192,6 +194,7 @@
   }
 }
 
+// ----------------------------------------------------------------------
 void
 pylith::faults::CohesiveTopology::orientFaultSieve(const int dim,
                                                    const Obj<Mesh>& mesh,
@@ -553,10 +556,6 @@
     Mesh::point_type cell      = cells[0];
     Mesh::point_type otherCell = cells[1];
 
-    if (flipFault) {
-      otherCell = cells[0];
-      cell      = cells[1];
-    }
     if (debug) std::cout << "  Checking orientation against cell " << cell << std::endl;
     ALE::ISieveTraversal<sieve_type>::orientedClosure(*ifaultSieve, face, cV2);
     const int               coneSize = cV2.getSize();
@@ -613,6 +612,12 @@
         }
       }
     }
+
+    if (flipFault) {
+      Mesh::point_type tmpCell = otherCell;
+      otherCell = cell;
+      cell      = tmpCell;
+    } // if
     if (found) {
       if (debug) std::cout << "  Choosing other cell" << std::endl;
       Mesh::point_type tmpCell = otherCell;
@@ -630,18 +635,18 @@
           if (faceVertices[v%numFaultCorners] != faceCone[c]) {
             found = false;
             break;
-          }
-        }
-      }
+          } // if
+        } // for
+      } // if
       if (!found) {
         std::cout << "Considering fault face " << face << std::endl;
         std::cout << "  bordered by cells " << cell << " and " << otherCell << std::endl;
         for(int c = 0; c < coneSize; ++c) {
           std::cout << "    Checking " << faceCone[c] << " against " << faceVertices[c] << std::endl;
-        }
-      }
+        } // for
+      } // if
       assert(found);
-    }
+    } // else
     noReplaceCells.insert(otherCell);
     replaceCells.insert(cell);
     replaceVertices.insert(faceCone, &faceCone[coneSize]);

Modified: short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.hh
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.hh	2009-01-15 22:42:24 UTC (rev 13872)
+++ short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.hh	2009-01-15 22:51:53 UTC (rev 13873)
@@ -186,8 +186,8 @@
               const Obj<Mesh>& mesh,
               const Obj<Mesh::int_section_type>& groupField,
               const int materialId,
-              const bool constraintCell = false,
-              const bool flipFault = false);
+              const bool constraintCell =false,
+              const bool flipFault =false);
 
   /** Create (distributed) fault mesh from cohesive cells.
    *
@@ -217,7 +217,7 @@
   static
   unsigned int _numFaceVertices(const Mesh::point_type& cell,
                                 const ALE::Obj<Mesh>& mesh,
-                                const int depth = -1);
+                                const int depth =-1);
 
   /** Determine a face orientation
    *    We should really have an interpolated mesh, instead of

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.cc	2009-01-15 22:42:24 UTC (rev 13872)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.cc	2009-01-15 22:51:53 UTC (rev 13873)
@@ -96,7 +96,7 @@
 { // testAdjustTopologyLine2
   CohesiveDataLine2 data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyLine2
 
 // ----------------------------------------------------------------------
@@ -106,7 +106,7 @@
 { // testAdjustTopologyTri3
   CohesiveDataTri3 data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, true);
 } // testAdjustTopologyTri3
 
 // ----------------------------------------------------------------------
@@ -116,7 +116,7 @@
 { // testAdjustTopologyTri3b
   CohesiveDataTri3b data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyTri3b
 
 // ----------------------------------------------------------------------
@@ -126,7 +126,7 @@
 { // testAdjustTopologyTri3c
   CohesiveDataTri3c data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyTri3c
 
 // ----------------------------------------------------------------------
@@ -136,7 +136,7 @@
 { // testAdjustTopologyTri3d
   CohesiveDataTri3d data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, true);
 } // testAdjustTopologyTri3d
 
 // ----------------------------------------------------------------------
@@ -146,7 +146,7 @@
 { // testAdjustTopologyTri3e
   CohesiveDataTri3e data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyTri3e
 
 // ----------------------------------------------------------------------
@@ -156,7 +156,7 @@
 { // testAdjustTopologyTri3f
   CohesiveDataTri3f data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, true);
 } // testAdjustTopologyTri3f
 
 // ----------------------------------------------------------------------
@@ -166,7 +166,7 @@
 { // testAdjustTopologyQuad4
   CohesiveDataQuad4 data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyQuad4
 
 // ----------------------------------------------------------------------
@@ -176,7 +176,7 @@
 { // testAdjustTopologyQuad4b
   CohesiveDataQuad4b data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyQuad4b
 
 // ----------------------------------------------------------------------
@@ -186,7 +186,7 @@
 { // testAdjustTopologyQuad4c
   CohesiveDataQuad4c data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, true);
 } // testAdjustTopologyQuad4c
 
 // ----------------------------------------------------------------------
@@ -196,7 +196,7 @@
 { // testAdjustTopologyQuad4d
   CohesiveDataQuad4d data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, true);
 } // testAdjustTopologyQuad4d
 
 // ----------------------------------------------------------------------
@@ -206,7 +206,7 @@
 { // testAdjustTopologyQuad4e
   CohesiveDataQuad4e data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, true);
 } // testAdjustTopologyQuad4e
 
 // ----------------------------------------------------------------------
@@ -216,7 +216,7 @@
 { // testAdjustTopologyQuad4f
   CohesiveDataQuad4f data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, true);
 } // testAdjustTopologyQuad4f
 
 // ----------------------------------------------------------------------
@@ -226,7 +226,7 @@
 { // testAdjustTopologyQuad4g
   CohesiveDataQuad4g data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyQuad4g
 
 // ----------------------------------------------------------------------
@@ -237,7 +237,7 @@
   CohesiveDataQuad4h data;
   FaultCohesiveDyn faultA;
   FaultCohesiveDyn faultB;
-  _testAdjustTopology(&faultA, &faultB, data);
+  _testAdjustTopology(&faultA, &faultB, data, false, false);
 } // testAdjustTopologyQuad4h
 
 // ----------------------------------------------------------------------
@@ -247,7 +247,7 @@
 { // testAdjustTopologyTet4
   CohesiveDataTet4 data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyTet4
 
 // ----------------------------------------------------------------------
@@ -257,7 +257,7 @@
 { // testAdjustTopologyTet4b
   CohesiveDataTet4b data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, true);
 } // testAdjustTopologyTet4b
 
 // ----------------------------------------------------------------------
@@ -267,7 +267,7 @@
 { // testAdjustTopologyTet4c
   CohesiveDataTet4c data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyTet4c
 
 // ----------------------------------------------------------------------
@@ -277,7 +277,7 @@
 { // testAdjustTopologyTet4d
   CohesiveDataTet4d data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyTet4d
 
 // ----------------------------------------------------------------------
@@ -287,7 +287,7 @@
 { // testAdjustTopologyTet4f
   CohesiveDataTet4f data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyTet4f
 
 // ----------------------------------------------------------------------
@@ -297,7 +297,7 @@
 { // testAdjustTopologyTet4g
   CohesiveDataTet4g data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyTet4g
 
 // ----------------------------------------------------------------------
@@ -307,7 +307,7 @@
 { // testAdjustTopologyTet4h
   CohesiveDataTet4h data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyTet4h
 
 // ----------------------------------------------------------------------
@@ -317,7 +317,7 @@
 { // testAdjustTopologyTet4i
   CohesiveDataTet4i data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyTet4i
 
 // ----------------------------------------------------------------------
@@ -327,7 +327,7 @@
 { // testAdjustTopologyTet4j
   CohesiveDataTet4j data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyTet4j
 
 // ----------------------------------------------------------------------
@@ -337,7 +337,7 @@
 { // testAdjustTopologyHex8
   CohesiveDataHex8 data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyHex8
 
 // ----------------------------------------------------------------------
@@ -347,7 +347,7 @@
 { // testAdjustTopologyHex8b
   CohesiveDataHex8b data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyHex8b
 
 // ----------------------------------------------------------------------
@@ -357,7 +357,7 @@
 { // testAdjustTopologyHex8c
   CohesiveDataHex8c data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyHex8c
 
 // ----------------------------------------------------------------------
@@ -367,7 +367,7 @@
 { // testAdjustTopologyHex8d
   CohesiveDataHex8d data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyHex8d
 
 // ----------------------------------------------------------------------
@@ -377,7 +377,7 @@
 { // testAdjustTopologyHex8e
   CohesiveDataHex8e data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyHex8e
 
 // ----------------------------------------------------------------------
@@ -387,7 +387,7 @@
 { // testAdjustTopologyHex8f
   CohesiveDataHex8f data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyHex8f
 
 // ----------------------------------------------------------------------
@@ -397,7 +397,7 @@
 { // testAdjustTopologyHex8g
   CohesiveDataHex8g data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyHex8g
 
 // ----------------------------------------------------------------------
@@ -407,7 +407,7 @@
 { // testAdjustTopologyHex8h
   CohesiveDataHex8h data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyHex8h
 
 // ----------------------------------------------------------------------
@@ -417,7 +417,7 @@
 { // testAdjustTopologyHex8i
   CohesiveDataHex8i data;
   FaultCohesiveDyn fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyHex8i
 
 // ----------------------------------------------------------------------
@@ -428,7 +428,7 @@
 { // testAdjustTopologyLine2Lagrange
   CohesiveDataLine2Lagrange data;
   FaultCohesiveKin fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyLine2Lagrange
 
 // ----------------------------------------------------------------------
@@ -439,7 +439,7 @@
 { // testAdjustTopologyTri3Lagrange
   CohesiveDataTri3Lagrange data;
   FaultCohesiveKin fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, true);
 } // testAdjustTopologyTri3Lagrange
 
 // ----------------------------------------------------------------------
@@ -450,7 +450,7 @@
 { // testAdjustTopologyQuad4Lagrange
   CohesiveDataQuad4Lagrange data;
   FaultCohesiveKin fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyQuad4Lagrange
 
 // ----------------------------------------------------------------------
@@ -461,7 +461,7 @@
 { // testAdjustTopologyTet4Lagrange
   CohesiveDataTet4Lagrange data;
   FaultCohesiveKin fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyTet4Lagrange
 
 // ----------------------------------------------------------------------
@@ -472,14 +472,15 @@
 { // testAdjustTopologyHex8Lagrange
   CohesiveDataHex8Lagrange data;
   FaultCohesiveKin fault;
-  _testAdjustTopology(&fault, data);
+  _testAdjustTopology(&fault, data, false);
 } // testAdjustTopologyHex8Lagrange
 
 // ----------------------------------------------------------------------
 // Test adjustTopology().
 void
 pylith::faults::TestFaultCohesive::_testAdjustTopology(Fault* fault,
-						      const CohesiveData& data)
+						       const CohesiveData& data,
+						       const bool flipFault)
 { // _testAdjustTopology
   ALE::Obj<Mesh> mesh;
   meshio::MeshIOAscii iohandler;
@@ -491,7 +492,8 @@
   CPPUNIT_ASSERT(0 != fault);
   fault->id(1);
   fault->label("fault");
-  fault->adjustTopology(mesh);
+  mesh->setDebug(true);
+  fault->adjustTopology(mesh, flipFault);
 
   CPPUNIT_ASSERT_EQUAL(data.cellDim, mesh->getDimension());
 
@@ -531,7 +533,7 @@
   ALE::ISieveVisitor::PointRetriever<Mesh::sieve_type> pV(sieve->getMaxConeSize());
   int iCell = 0;
   i = 0;
-  //mesh->view(data.filename);
+  mesh->view(data.filename);
   for(Mesh::label_sequence::iterator c_iter = cells->begin();
       c_iter != cells->end();
       ++c_iter) {
@@ -597,7 +599,9 @@
 void
 pylith::faults::TestFaultCohesive::_testAdjustTopology(Fault* faultA,
 						       Fault* faultB,
-						      const CohesiveData& data)
+						       const CohesiveData& data,
+						       const bool flipFaultA,
+						       const bool flipFaultB)
 { // _testAdjustTopology
   ALE::Obj<Mesh> mesh;
   meshio::MeshIOAscii iohandler;
@@ -609,12 +613,12 @@
   CPPUNIT_ASSERT(0 != faultA);
   faultA->id(1);
   faultA->label("faultA");
-  faultA->adjustTopology(mesh);
+  faultA->adjustTopology(mesh, flipFaultA);
 
   CPPUNIT_ASSERT(0 != faultB);
   faultB->id(1);
   faultB->label("faultB");
-  faultB->adjustTopology(mesh);
+  faultB->adjustTopology(mesh, flipFaultB);
 
   CPPUNIT_ASSERT_EQUAL(data.cellDim, mesh->getDimension());
 

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.hh
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.hh	2009-01-15 22:42:24 UTC (rev 13872)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.hh	2009-01-15 22:51:53 UTC (rev 13873)
@@ -221,19 +221,25 @@
    *
    * @param fault Fault for cohesive elements.
    * @param data Cohesive element data.
+   * @param flipFault If true, add vertices to opposite side of fault.
    */
   void _testAdjustTopology(Fault* fault,
-			   const CohesiveData& data);
+			   const CohesiveData& data,
+			   const bool flipFault);
 
   /** Test adjustTopology().
    *
    * @param faultA First fault for cohesive elements.
    * @param faultB Second fault for cohesive elements.
    * @param data Cohesive element data.
+   * @param flipFaultA If true, add vertices to opposite side of fault A.
+   * @param flipFaultB If true, add vertices to opposite side of fault B.
    */
   void _testAdjustTopology(Fault* faultA,
 			   Fault* faultB,
-			   const CohesiveData& data);
+			   const CohesiveData& data,
+			   const bool flipFaultA,
+			   const bool flipFaultB);
 
 }; // class TestFaultCohesive
 

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc	2009-01-15 22:42:24 UTC (rev 13872)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc	2009-01-15 22:51:53 UTC (rev 13873)
@@ -46,6 +46,8 @@
   _eqsrcs[0] = new EqKinSrc();
   _slipfns.resize(nsrcs);
   _slipfns[0] = new BruneSlipFn();
+
+  _flipFault = false;
 } // setUp
 
 // ----------------------------------------------------------------------
@@ -130,6 +132,8 @@
   ALE::Obj<Mesh> mesh;
   FaultCohesiveKin fault;
   _initialize(&mesh, &fault);
+
+  //mesh->view(_data->meshFilename);
   
   Mesh::renumbering_type& renumbering = fault._faultMesh->getRenumbering();
   const ALE::Obj<Mesh::label_sequence>& vertices = 
@@ -909,7 +913,7 @@
     fault->quadrature(_quadrature);
     
     fault->eqsrcs(const_cast<const char**>(names), sources, nsrcs);
-    fault->adjustTopology(*mesh);
+    fault->adjustTopology(*mesh, _flipFault);
 
     const double upDirVals[] = { 0.0, 0.0, 1.0 };
     double_array upDir(upDirVals, 3);

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.hh
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.hh	2009-01-15 22:42:24 UTC (rev 13872)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.hh	2009-01-15 22:51:53 UTC (rev 13873)
@@ -65,6 +65,7 @@
   feassemble::Quadrature* _quadrature; ///< Data used in testing
   std::vector<EqKinSrc*> _eqsrcs; ///< Array of Kinematic earthquake sources.
   std::vector<BruneSlipFn*> _slipfns; ///< Slip time function.
+  bool _flipFault; ///< If true, dd new vertices to opposite side of fault
 
   // PUBLIC METHODS /////////////////////////////////////////////////////
 public :

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinSrcs.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinSrcs.cc	2009-01-15 22:42:24 UTC (rev 13872)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinSrcs.cc	2009-01-15 22:51:53 UTC (rev 13873)
@@ -47,6 +47,8 @@
   _slipfns.resize(nsrcs);
   _slipfns[0] = new BruneSlipFn();
   _slipfns[1] = new BruneSlipFn();
+
+  _flipFault = false;
 } // setUp
 
 

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinTri3.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinTri3.cc	2009-01-15 22:42:24 UTC (rev 13872)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinTri3.cc	2009-01-15 22:51:53 UTC (rev 13873)
@@ -33,6 +33,8 @@
   CPPUNIT_ASSERT(0 != _quadrature);
   feassemble::GeometryLine2D geometry;
   _quadrature->refGeometry(&geometry);
+  
+  _flipFault = true;
 } // setUp
 
 

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataTri3.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataTri3.cc	2009-01-15 22:42:24 UTC (rev 13872)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataTri3.cc	2009-01-15 22:51:53 UTC (rev 13873)
@@ -31,17 +31,6 @@
  *
  * Cells are 0-1, 8, vertices are 2-7.
  *
- *              6 -7- 3
- *             /|     |\
- *            / |     | \
- *           /  |     |  \
- *          /   |     |   \
- *         2    |     |    5
- *          \   |     |   /
- *           \  |     |  /
- *            \ |     | /
- *             \|     |/
- *              8 -9- 4
  *              6 -8- 3
  *             /|     |\
  *            / |     | \



More information about the CIG-COMMITS mailing list