[cig-commits] r12988 - short/3D/PyLith/trunk/libsrc/faults
brad at geodynamics.org
brad at geodynamics.org
Thu Oct 2 21:04:11 PDT 2008
Author: brad
Date: 2008-10-02 21:04:11 -0700 (Thu, 02 Oct 2008)
New Revision: 12988
Modified:
short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.cc
Log:
Commented out debugging stuff.
Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.cc 2008-10-03 04:01:22 UTC (rev 12987)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.cc 2008-10-03 04:04:11 UTC (rev 12988)
@@ -622,14 +622,17 @@
ncV.clear();
} // for
+#if 0
_orientation->view("ORIENTATION Before complete");
_orientation->setDebug(2);
+#endif
// Assemble orientation information
ALE::Completion::completeSectionAdd(_faultMesh->getSendOverlap(), _faultMesh->getRecvOverlap(), _orientation, _orientation);
+
+#if 0
_orientation->view("ORIENTATION After complete");
+#endif
- _orientation->view("ORIENTATION");
-
// Loop over vertices, make orientation information unit magnitude
double_array vertexDir(orientationSize);
int count = 0;
@@ -697,7 +700,7 @@
PetscLogFlops(5 + count * 3);
} // if
- _orientation->view("ORIENTATION");
+ //_orientation->view("ORIENTATION");
} // _calcOrientation
// ----------------------------------------------------------------------
@@ -833,7 +836,7 @@
// Assemble area information
ALE::Completion::completeSectionAdd(_faultMesh->getSendOverlap(), _faultMesh->getRecvOverlap(), _area, _area);
- _area->view("AREA");
+ //_area->view("AREA");
} // _calcArea
// ----------------------------------------------------------------------
@@ -851,8 +854,10 @@
assert(!_pseudoStiffness.isNull());
assert(!_area.isNull());
+#if 0
_pseudoStiffness->view("PSEUDOSTIFFNESS");
_area->view("AREA");
+#endif
const ALE::Obj<Mesh::label_sequence>& vertices = mesh->depthStratum(0);
assert(!vertices.isNull());
@@ -912,10 +917,13 @@
PetscLogFlops(numVertices * (1 + fiberDim) );
+#if 0
_faultMesh->view("FAULT MESH");
solution->view("SOLUTION");
_area->view("AREA");
+ _pseudoStiffness->view("CONDITIONING");
(*tractions)->view("TRACTIONS");
+#endif
} // _calcTractionsChange
// ----------------------------------------------------------------------
More information about the cig-commits
mailing list