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

[brad at geodynamics.org]

Author: brad
Date: 2010-03-09 16:03:42 -0800 (Tue, 09 Mar 2010)
New Revision: 16396

Modified:
   short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc
Log:
Turned off some debugging output.

Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc	2010-03-09 22:08:59 UTC (rev 16395)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc	2010-03-10 00:03:42 UTC (rev 16396)
@@ -199,8 +199,6 @@
   } // for
 
   PetscLogFlops(numVertices*spaceDim);
-
-  residual.view("RESIDUAL FAULT");
 } // integrateResidualAssembled
 
 // ----------------------------------------------------------------------
@@ -1250,14 +1248,14 @@
       for (int kDim = 0; kDim < spaceDim; ++kDim)
         forcesInitialVertexGlobal[iDim] +=
           forcesInitialVertexFault[kDim] * 
-	  orientationVertex[iDim*spaceDim+kDim];
+	  orientationVertex[kDim*spaceDim+iDim];
 
     assert(forcesInitialVertexGlobal.size() == 
 	   forcesInitialSection->getFiberDimension(v_fault));
     forcesInitialSection->updatePoint(v_fault, &forcesInitialVertexGlobal[0]);
   } // for
 
-  forcesInitial.view("INITIAL FORCES"); // DEBUGGING
+  //forcesInitial.view("INITIAL FORCES"); // DEBUGGING
 } // _setupInitialTractions
 
 // ----------------------------------------------------------------------
@@ -1394,7 +1392,7 @@
     for (int iDim = 0; iDim < spaceDim; ++iDim)
       for (int kDim = 0; kDim < spaceDim; ++kDim)
         tractionsVertexFault[iDim] +=
-          tractionsVertexGlobal[kDim] * orientationVertex[kDim*spaceDim+iDim];
+          tractionsVertexGlobal[kDim] * orientationVertex[iDim*spaceDim+kDim];
     
     assert(tractionsVertexFault.size() == 
 	   tractionsSection->getFiberDimension(v_fault));