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

brad at geodynamics.org brad at geodynamics.org
Fri Mar 8 15:45:30 PST 2013


Author: brad
Date: 2013-03-08 15:45:29 -0800 (Fri, 08 Mar 2013)
New Revision: 21477

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc
Log:
Code cleanup.

Modified: short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc	2013-03-08 23:43:50 UTC (rev 21476)
+++ short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc	2013-03-08 23:45:29 UTC (rev 21477)
@@ -613,11 +613,6 @@
   } // switch
 
 
-#if 1 // DEBUGGING
-   _fields->get("relative disp").view("BEFORE RELATIVE DISPLACEMENT");
-   fields->get("dispIncr(t->t+dt)").view("BEFORE DISP INCR (t->t+dt)");
-#endif
-
   err = VecGetArray(dispTVec, &dispTArray);CHECK_PETSC_ERROR(err);
   err = VecGetArray(dispTIncrVec, &dispTIncrArray);CHECK_PETSC_ERROR(err);
   err = VecGetArray(dispTIncrAdjVec, &dispTIncrAdjArray);CHECK_PETSC_ERROR(err);
@@ -853,7 +848,7 @@
       // if residual is very small, we prefer the full step
       break;
 
-#if 0
+#if 0 // DEBUGGING
     const int rank = _faultMesh->sieveMesh()->commRank();
     std::cout << "["<<rank<<"] alphaL: " << pow(10.0, logAlphaL)
 	      << ", residuaL: " << residualL
@@ -2089,13 +2084,6 @@
   PetscVec residualVec = residual.localVector();assert(residualVec);
   residual.zero();
 
-#if 0
-  std::cout << "dLagrangeVec" << std::endl;
-  VecView(dLagrangeVec, PETSC_VIEWER_STDOUT_WORLD);
-  std::cout << "SENSITIVITY mesh: " << faultDMMesh << ", coordinates: " << coordVec << std::endl;
-  VecView(coordVec, PETSC_VIEWER_STDOUT_WORLD);
-#endif
-
   // Loop over cells
   for(PetscInt c = cStart; c < cEnd; ++c) {
     // Compute geometry
@@ -2178,7 +2166,7 @@
   // Update section view of field.
   solution.scatterVectorToSection();
 
-#if 1 // DEBUGGING
+#if 0 // DEBUGGING
   residual.view("SENSITIVITY RESIDUAL");
   solution.view("SENSITIVITY SOLUTION");
 #endif



More information about the CIG-COMMITS mailing list