[cig-commits] [commit] knepley/upgrade-petsc-interface: Cleanup of names for Lagrange edges. (086b7b9)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Fri Nov 1 10:20:23 PDT 2013
Repository : ssh://geoshell/pylith
On branch : knepley/upgrade-petsc-interface
Link : https://github.com/geodynamics/pylith/compare/4fd816f765cf48e0d10c86110d8fa9fcba507d4e...086b7b94b37b71e66c6947911ac1b4bf11480589
>---------------------------------------------------------------
commit 086b7b94b37b71e66c6947911ac1b4bf11480589
Author: Brad Aagaard <baagaard at usgs.gov>
Date: Fri Nov 1 10:22:44 2013 -0700
Cleanup of names for Lagrange edges.
Also turned off debugging output.
>---------------------------------------------------------------
086b7b94b37b71e66c6947911ac1b4bf11480589
unittests/libtests/faults/TestFaultCohesiveKin.cc | 20 ++++++++++----------
unittests/libtests/faults/TestFaultCohesiveKin.hh | 8 ++++----
2 files changed, 14 insertions(+), 14 deletions(-)
diff --git a/unittests/libtests/faults/TestFaultCohesiveKin.cc b/unittests/libtests/faults/TestFaultCohesiveKin.cc
index ff40595..602ae42 100644
--- a/unittests/libtests/faults/TestFaultCohesiveKin.cc
+++ b/unittests/libtests/faults/TestFaultCohesiveKin.cc
@@ -217,7 +217,7 @@ pylith::faults::TestFaultCohesiveKin::testInitialize(void)
} // for
// Check orientation
- fault._fields->get("orientation").view("ORIENTATION"); // DEBUGGING
+ //fault._fields->get("orientation").view("ORIENTATION"); // DEBUGGING
topology::VecVisitorMesh orientationVisitor(fault._fields->get("orientation"));
const PetscScalar* orientationArray = orientationVisitor.localArray();CPPUNIT_ASSERT(orientationArray);
@@ -271,7 +271,7 @@ pylith::faults::TestFaultCohesiveKin::testIntegrateResidual(void)
topology::Field& residual = fields.get("residual");
fault.integrateResidual(residual, t, &fields);
- residual.view("RESIDUAL"); // DEBUGGING
+ //residual.view("RESIDUAL"); // DEBUGGING
// Check values
CPPUNIT_ASSERT(_data->residual);
@@ -364,7 +364,7 @@ pylith::faults::TestFaultCohesiveKin::testIntegrateJacobian(void)
for (int iCol=0; iCol < ncols; ++iCol) {
const int index = ncols*iRow+iCol;
const PylithScalar valE = valsE[index];
-#if 1 // DEBUGGING
+#if 0 // DEBUGGING
if (fabs(valE-vals[index]*jacobianScale) > tolerance)
std::cout << "ERROR: iRow: " << iRow << ", iCol: " << iCol
<< "valE: " << valE
@@ -406,7 +406,7 @@ pylith::faults::TestFaultCohesiveKin::testIntegrateJacobianLumped(void)
CPPUNIT_ASSERT_EQUAL(false, fault.needNewJacobian());
jacobian.complete();
- jacobian.view("JACOBIAN"); // DEBUGGING
+ //jacobian.view("JACOBIAN"); // DEBUGGING
const int spaceDim = _data->spaceDim;
@@ -482,7 +482,7 @@ pylith::faults::TestFaultCohesiveKin::testAdjustSolnLumped(void)
const topology::Field& dispIncrAdj = fields.get("dispIncr adjust");
solution += dispIncrAdj;
- solution.view("SOLUTION AFTER ADJUSTMENT"); // DEBUGGING
+ //solution.view("SOLUTION AFTER ADJUSTMENT"); // DEBUGGING
PetscInt pStart, pEnd;
PetscErrorCode err = PetscSectionGetChart(jacobian.petscSection(), &pStart, &pEnd);PYLITH_CHECK_ERROR(err);
@@ -740,10 +740,10 @@ pylith::faults::TestFaultCohesiveKin::_initialize(topology::Mesh* const mesh,
} // _initialize
// ----------------------------------------------------------------------
-// Determine if vertex is a Lagrange multiplier constraint vertex.
+// Determine if point is a Lagrange multiplier constraint point.
bool
-pylith::faults::TestFaultCohesiveKin::_isConstraintVertex(const int vertex) const
-{ // _isConstraintVertex
+pylith::faults::TestFaultCohesiveKin::_isConstraintEdge(const int point) const
+{ // _isConstraintEdge
PYLITH_METHOD_BEGIN;
CPPUNIT_ASSERT(_data);
@@ -751,12 +751,12 @@ pylith::faults::TestFaultCohesiveKin::_isConstraintVertex(const int vertex) cons
const int numFaultVertices = _data->numFaultVertices;
bool isFound = false;
for (int i=0; i < _data->numFaultVertices; ++i)
- if (_data->edgesLagrange[i] == vertex) {
+ if (_data->edgesLagrange[i] == point) {
isFound = true;
break;
} // if
PYLITH_METHOD_RETURN(isFound);
-} // _isConstraintVertex
+} // _isConstraintEdge
// End of file
diff --git a/unittests/libtests/faults/TestFaultCohesiveKin.hh b/unittests/libtests/faults/TestFaultCohesiveKin.hh
index f1f07b9..4000235 100644
--- a/unittests/libtests/faults/TestFaultCohesiveKin.hh
+++ b/unittests/libtests/faults/TestFaultCohesiveKin.hh
@@ -129,13 +129,13 @@ private :
FaultCohesiveKin* const fault,
topology::SolutionFields* const fields) const;
- /** Determine if vertex is a Lagrange multiplier constraint vertex.
+ /** Determine if point is a Lagrange multiplier constraint point.
*
- * @param vertex Label of vertex.
+ * @param point Label of point.
*
- * @returns True if vertex is a constraint vertex, false otherwise.
+ * @returns True if point is a constraint point, false otherwise.
*/
- bool _isConstraintVertex(const int vertex) const;
+ bool _isConstraintEdge(const int point) const;
}; // class TestFaultCohesiveKin
More information about the CIG-COMMITS
mailing list