[cig-commits] r14708 - in short/3D/PyLith/branches/pylith-swig: . tests/1d/line2 unittests/libtests/meshio

[knepley at geodynamics.org]

Author: knepley
Date: 2009-04-14 17:19:58 -0700 (Tue, 14 Apr 2009)
New Revision: 14708

Modified:
   short/3D/PyLith/branches/pylith-swig/TODO
   short/3D/PyLith/branches/pylith-swig/tests/1d/line2/axialextension.cfg
   short/3D/PyLith/branches/pylith-swig/tests/1d/line2/testpylith.py
   short/3D/PyLith/branches/pylith-swig/unittests/libtests/meshio/TestDataWriterVTKBCMesh.cc
Log:
Fixed BC test, small changes


Modified: short/3D/PyLith/branches/pylith-swig/TODO
===================================================================
--- short/3D/PyLith/branches/pylith-swig/TODO	2009-04-14 23:58:42 UTC (rev 14707)
+++ short/3D/PyLith/branches/pylith-swig/TODO	2009-04-15 00:19:58 UTC (rev 14708)
@@ -77,6 +77,7 @@
   pytests.meshio.TestOutputNeumann
 
 1. Memory model [Matt]
+   playpen/memcheck/*
 
 2. Nondimensionalize [Brad]
     Not working
@@ -114,12 +115,10 @@
 
 10. PETSc SNES nonlinear solvers [Matt]
 
-11. Create spatialdata/units/UnitsParser (call Pyre units parser).
+11. Add ability to use symmetric matrix, set block size (default is 1)
 
-12. Add ability to use symmetric matrix, set block size (default is 1)
+12. Update manual for version 1.4 [Brad and Charles]
 
-13. Update manual for version 1.4 [Brad and Charles]
-
   * Reduce memory use with ordering elements by material
 
   * Nondimensionalization
@@ -137,6 +136,7 @@
 -- Release (v1.4) --
 
 Modularize output? [Matt]
+  "field_dimensioned" is prefered to "field"
 
 Fault friction
 

Modified: short/3D/PyLith/branches/pylith-swig/tests/1d/line2/axialextension.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-swig/tests/1d/line2/axialextension.cfg	2009-04-14 23:58:42 UTC (rev 14707)
+++ short/3D/PyLith/branches/pylith-swig/tests/1d/line2/axialextension.cfg	2009-04-15 00:19:58 UTC (rev 14708)
@@ -1,6 +1,9 @@
 # -*- Python -*-
 [axialextension]
 
+[axialextension.launcher]
+command = mpiexec -n ${nodes}
+
 # ----------------------------------------------------------------------
 # journal
 # ----------------------------------------------------------------------
@@ -31,7 +34,7 @@
 dimension = 1
 bc = [bc]
 
-#formulation.solver = pylith.problems.SolverNonlinear
+formulation.solver = pylith.problems.SolverNonlinear
 
 [axialextension.timedependent.formulation.time_step]
 total_time = 0.0*s

Modified: short/3D/PyLith/branches/pylith-swig/tests/1d/line2/testpylith.py
===================================================================
--- short/3D/PyLith/branches/pylith-swig/tests/1d/line2/testpylith.py	2009-04-14 23:58:42 UTC (rev 14707)
+++ short/3D/PyLith/branches/pylith-swig/tests/1d/line2/testpylith.py	2009-04-15 00:19:58 UTC (rev 14708)
@@ -1,4 +1,4 @@
-#!/usr/bin/env python
+#!/usr/bin/env nemesis
 #
 # ======================================================================
 #

Modified: short/3D/PyLith/branches/pylith-swig/unittests/libtests/meshio/TestDataWriterVTKBCMesh.cc
===================================================================
--- short/3D/PyLith/branches/pylith-swig/unittests/libtests/meshio/TestDataWriterVTKBCMesh.cc	2009-04-14 23:58:42 UTC (rev 14707)
+++ short/3D/PyLith/branches/pylith-swig/unittests/libtests/meshio/TestDataWriterVTKBCMesh.cc	2009-04-15 00:19:58 UTC (rev 14708)
@@ -288,10 +288,11 @@
     const ALE::Obj<topology::SubMesh::SieveMesh>& sieveSubMesh =
       _submesh->sieveMesh();
     CPPUNIT_ASSERT(!sieveSubMesh.isNull());
-    const int cellDepth = sieveSubMesh->depth()-1;
     const ALE::Obj<topology::SubMesh::SieveMesh::label_sequence>& cells = 
-      sieveSubMesh->depthStratum(cellDepth);
+      sieveSubMesh->heightStratum(1);
     assert(!cells.isNull());
+    const topology::SubMesh::SieveMesh::label_sequence::iterator cellsBegin = 
+      cells->begin();
     const topology::SubMesh::SieveMesh::label_sequence::iterator cellsEnd = 
       cells->end();
 
@@ -308,7 +309,7 @@
       const ALE::Obj<topology::SubMesh::RealSection>& section = field.section();
       CPPUNIT_ASSERT(!section.isNull());
       int icell = 0;
-      for (topology::SubMesh::SieveMesh::label_sequence::iterator c_iter=cells->begin();
+      for (topology::SubMesh::SieveMesh::label_sequence::iterator c_iter=cellsBegin;
 	   c_iter != cellsEnd;
 	   ++c_iter, ++icell) {
 	const double* values = &_data->cellFields[i][icell*fiberDim];