[cig-commits] r16524 - in short/3D/PyLith/tags/v1.4.3: . libsrc/bc libsrc/faults libsrc/feassemble libsrc/materials libsrc/meshio libsrc/topology libsrc/utils modulesrc/faults modulesrc/meshio pylith/faults unittests/libtests/bc unittests/libtests/faults unittests/libtests/meshio unittests/libtests/meshio/data unittests/libtests/topology unittests/libtests/utils unittests/pytests/faults unittests/pytests/utils

brad at geodynamics.org brad at geodynamics.org
Fri Apr 9 17:35:29 PDT 2010


Author: brad
Date: 2010-04-09 17:35:22 -0700 (Fri, 09 Apr 2010)
New Revision: 16524

Modified:
   short/3D/PyLith/tags/v1.4.3/TODO
   short/3D/PyLith/tags/v1.4.3/configure.ac
   short/3D/PyLith/tags/v1.4.3/libsrc/bc/AbsorbingDampers.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/bc/BoundaryConditionPoints.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/bc/Neumann.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/faults/CohesiveTopology.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/faults/CohesiveTopology.hh
   short/3D/PyLith/tags/v1.4.3/libsrc/faults/Fault.hh
   short/3D/PyLith/tags/v1.4.3/libsrc/faults/FaultCohesive.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/faults/FaultCohesive.hh
   short/3D/PyLith/tags/v1.4.3/libsrc/faults/FaultCohesiveKin.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/faults/TopologyOps.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/faults/TopologyVisitors.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/feassemble/ElasticityExplicit.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/feassemble/ElasticityImplicit.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/feassemble/IntegratorElasticity.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/materials/ElasticMaterial.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/meshio/CellFilterAvg.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/meshio/DataWriterVTK.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/meshio/DataWriterVTK.hh
   short/3D/PyLith/tags/v1.4.3/libsrc/meshio/MeshIO.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/topology/Mesh.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/topology/SubMesh.cc
   short/3D/PyLith/tags/v1.4.3/libsrc/utils/EventLogger.cc
   short/3D/PyLith/tags/v1.4.3/modulesrc/faults/Fault.i
   short/3D/PyLith/tags/v1.4.3/modulesrc/faults/FaultCohesive.i
   short/3D/PyLith/tags/v1.4.3/modulesrc/meshio/DataWriterVTK.i
   short/3D/PyLith/tags/v1.4.3/pylith/faults/LiuCosSlipFn.py
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/bc/TestBoundaryMesh.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestBruneSlipFn.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestConstRateSlipFn.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestEqKinSrc.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestFaultCohesive.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestFaultCohesiveKin.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestLiuCosSlipFn.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestStepSlipFn.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKBCMesh.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKFaultMesh.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKMesh.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKMesh.hh
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKSubMesh.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_bc_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_bc_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_bc_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_fault_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_fault_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_fault_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_mat_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_mat_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_mat_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_surf_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_surf_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_surf_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_mat_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_mat_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_mat_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_surf_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_surf_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_surf_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_cell_filter_t12.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_cell_t12.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_t12.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_vertex_filter_t12.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_vertex_t12.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_bc_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_bc_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_bc_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_fault_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_fault_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_fault_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_mat_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_mat_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_mat_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_surf_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_surf_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_surf_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_bc_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_bc_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_bc_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_fault_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_fault_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_fault_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_mat_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_mat_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_mat_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_surf_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_surf_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_surf_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_bc_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_bc_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_bc_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_fault_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_fault_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_fault_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_mat_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_mat_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_mat_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_surf_cell_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_surf_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_surf_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_vertex_t10.vtk
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/topology/TestMeshOps.cc
   short/3D/PyLith/tags/v1.4.3/unittests/libtests/utils/TestEventLogger.cc
   short/3D/PyLith/tags/v1.4.3/unittests/pytests/faults/TestFaultCohesiveKin.py
   short/3D/PyLith/tags/v1.4.3/unittests/pytests/utils/TestEventLogger.py
Log:
Manually merged bugfixes from trunk. Main bugfix is retrieving physical properties when updating bulk state variables (Roundup Issue 546). Also fixed bug related to dimensioning fault coordinates (Roundup Issue Issues 548) and filtering cell output for submeshes. Also merged some error checking and setting precision of output in VTK files (C++ only).

Modified: short/3D/PyLith/tags/v1.4.3/TODO
===================================================================
--- short/3D/PyLith/tags/v1.4.3/TODO	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/TODO	2010-04-10 00:35:22 UTC (rev 16524)
@@ -2,6 +2,8 @@
 CURRENT ISSUES/PRIORITIES
 ======================================================================
 
+Need to incorporate bugfixes from trunk for Powerlaw3D??
+
 ----------------------------------------------------------------------
 POST RELEASE 1.4
 ----------------------------------------------------------------------

Modified: short/3D/PyLith/tags/v1.4.3/configure.ac
===================================================================
--- short/3D/PyLith/tags/v1.4.3/configure.ac	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/configure.ac	2010-04-10 00:35:22 UTC (rev 16524)
@@ -109,7 +109,7 @@
 #CIT_CHECK_LIB_MPI
 
 # PETSC
-CIT_PATH_PETSC([3.0.0])
+CIT_PATH_PETSC([3.1.0])
 CIT_HEADER_PETSC
 CIT_CHECK_LIB_PETSC
 CIT_CHECK_LIB_PETSC_SIEVE

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/bc/AbsorbingDampers.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/bc/AbsorbingDampers.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/bc/AbsorbingDampers.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -357,7 +357,7 @@
 
     // Assemble cell contribution into field
     residualVisitor.clear();
-    sieveSubMesh->updateAdd(*c_iter, residualVisitor);
+    sieveSubMesh->updateClosure(*c_iter, residualVisitor);
   } // for
 } // integrateResidual
 

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/bc/BoundaryConditionPoints.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/bc/BoundaryConditionPoints.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/bc/BoundaryConditionPoints.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -67,13 +67,14 @@
   const ALE::Obj<SieveMesh>& sieveMesh = mesh.sieveMesh();
   assert(!sieveMesh.isNull());
 
-  const ALE::Obj<SieveMesh::int_section_type>& groupField = 
-    sieveMesh->getIntSection(_label);
-  if (groupField.isNull()) {
+  if (!sieveMesh->hasIntSection(_label)) {
     std::ostringstream msg;
     msg << "Could not find group of points '" << _label << "' in mesh.";
     throw std::runtime_error(msg.str());
   } // if
+
+  const ALE::Obj<SieveMesh::int_section_type>& groupField = 
+    sieveMesh->getIntSection(_label);
   assert(!groupField.isNull());
   const chart_type& chart = groupField->getChart();
   const chart_type::const_iterator& chartEnd = chart.end();

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/bc/Neumann.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/bc/Neumann.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/bc/Neumann.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -158,7 +158,7 @@
     } // for
     // Assemble cell contribution into field
     residualVisitor.clear();
-    subSieveMesh->updateAdd(*c_iter, residualVisitor);
+    subSieveMesh->updateClosure(*c_iter, residualVisitor);
 
     PetscLogFlops(numQuadPts*(1+numBasis*(1+numBasis*(1+2*spaceDim))));
   } // for

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/faults/CohesiveTopology.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/faults/CohesiveTopology.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/faults/CohesiveTopology.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -124,6 +124,7 @@
                                          const ALE::Obj<topology::Mesh::IntSection>& groupField,
                                          const int materialId,
                                          int& firstFaultVertex,
+                                         int& firstLagrangeVertex,
                                          int& firstFaultCell,
                                          const bool constraintCell)
 { // create
@@ -188,10 +189,12 @@
   assert(!groupNames.isNull());
   const int numFaultVertices = fVertices->size();
   std::map<point_type,point_type> vertexRenumber;
+  std::map<point_type,point_type> vertexLagrangeRenumber;
   std::map<point_type,point_type> cellRenumber;
   if (firstFaultVertex == 0) {
-    firstFaultVertex += sieve->getBaseSize() + sieve->getCapSize();
-    firstFaultCell   += firstFaultVertex;
+    firstFaultVertex    += sieve->getBaseSize() + sieve->getCapSize();
+    firstLagrangeVertex += firstFaultVertex;
+    firstFaultCell      += firstFaultVertex;
   }
 
   for(SieveSubMesh::label_sequence::iterator v_iter = fVerticesBegin;
@@ -213,12 +216,14 @@
     sieveMesh->setDepth(firstFaultVertex, 0);
 #endif
     if (constraintCell) {
-      groupField->addPoint(firstFaultVertex+numFaultVertices, 1);
+      vertexLagrangeRenumber[*v_iter] = firstLagrangeVertex;
+      groupField->addPoint(firstLagrangeVertex, 1);
 #if defined(FAST_STRATIFY)
       // OPTIMIZATION
-      sieveMesh->setHeight(firstFaultVertex+numFaultVertices, 1);
-      sieveMesh->setDepth(firstFaultVertex+numFaultVertices, 0);
+      sieveMesh->setHeight(firstLagrangeVertex, 1);
+      sieveMesh->setDepth(firstLagrangeVertex, 0);
 #endif
+      ++firstLagrangeVertex;
     } // if
     logger.stagePop();
     logger.stagePush("FaultCreation");
@@ -242,9 +247,6 @@
       ++name) {
     sieveMesh->reallocate(sieveMesh->getIntSection(*name));
   } // for
-  if (constraintCell) {
-    firstFaultVertex += numFaultVertices;
-  }
 
   // Split the mesh along the fault sieve and create cohesive elements
   const ALE::Obj<SieveSubMesh::label_sequence>& faces =
@@ -395,8 +397,8 @@
     if (constraintCell) {
       for (int c = 0; c < coneSize; ++c) {
         if (debug)
-	  std::cout << "    Lagrange vertex " << vertexRenumber[faceCone[c]]+numFaultVertices << std::endl;
-        sieve->addArrow(vertexRenumber[faceCone[c]]+numFaultVertices, firstFaultCell, true);
+	  std::cout << "    Lagrange vertex " << vertexLagrangeRenumber[faceCone[c]] << std::endl;
+        sieve->addArrow(vertexLagrangeRenumber[faceCone[c]], firstFaultCell, true);
       } // for
     } // if
     // TODO: Need to reform the material label when sieve is reallocated
@@ -609,8 +611,7 @@
     const ALE::Obj<Mesh::label_type>& label = sieveMesh->createLabel(labelName);
     assert(!label.isNull());
 
-    TopologyOps::computeCensoredDepth(label, sieveMesh->getSieve(),
-				      firstCohesiveCell-(constraintCell?numFaultVertices:0));
+    TopologyOps::computeCensoredDepth(label, sieveMesh->getSieve(), firstFaultVertex);
   } else {
     // Insert new shadow vertices into existing label
     const ALE::Obj<Mesh::label_type>& label = sieveMesh->getLabel(labelName);
@@ -645,7 +646,7 @@
     coordinates->addPoint(vertexRenumber[*v_iter],
 			  coordinates->getFiberDimension(*v_iter));
     if (constraintCell)
-      coordinates->addPoint(vertexRenumber[*v_iter]+numFaultVertices,
+      coordinates->addPoint(vertexLagrangeRenumber[*v_iter],
 			    coordinates->getFiberDimension(*v_iter));
   } // for
   sieveMesh->reallocate(coordinates);
@@ -657,7 +658,7 @@
     coordinates->updatePoint(vertexRenumber[*v_iter], 
 			     coordinates->restrictPoint(*v_iter));
     if (constraintCell)
-      coordinates->updatePoint(vertexRenumber[*v_iter]+numFaultVertices,
+      coordinates->updatePoint(vertexLagrangeRenumber[*v_iter],
 			       coordinates->restrictPoint(*v_iter));
   } // for
   if (debug)
@@ -800,14 +801,23 @@
   assert(!fCoordinates.isNull());
   const ALE::Obj<SieveMesh::label_sequence>& vertices =
     sieveMesh->depthStratum(0);
+  assert(!vertices.isNull());
   const SieveMesh::label_sequence::iterator vBegin = vertices->begin();
   const SieveMesh::label_sequence::iterator vEnd = vertices->end();
 
-  fCoordinates->setChart(topology::Mesh::RealSection::chart_type(faultSieveMesh->heightStratum(0)->size(),
-                                                             faultSieveMesh->getSieve()->getChart().max()));
+  const ALE::Obj<SieveMesh::label_sequence>& fvertices =
+    faultSieveMesh->depthStratum(0);
+  assert(!fvertices.isNull());
+  const point_type fvMin = *std::min_element(fvertices->begin(), 
+					     fvertices->end());
+  const point_type fvMax = *std::max_element(fvertices->begin(),
+					     fvertices->end());
+  
+  fCoordinates->setChart(topology::Mesh::RealSection::chart_type(fvMin, 
+								 fvMax+1));
   for (SieveMesh::label_sequence::iterator v_iter = vBegin;
-      v_iter != vEnd;
-      ++v_iter) {
+       v_iter != vEnd;
+       ++v_iter) {
     if (fRenumbering.find(*v_iter) == fRenumberingEnd)
       continue;
     fCoordinates->setFiberDimension(fRenumbering[*v_iter],
@@ -832,13 +842,9 @@
     const ALE::Obj<topology::Mesh::RealSection>& fCoordinatesDim =
       faultSieveMesh->getRealSection("coordinates_dimensioned");
     assert(!fCoordinatesDim.isNull());
-    const ALE::Obj<SieveMesh::label_sequence>& vertices =
-      sieveMesh->depthStratum(0);
-    const SieveMesh::label_sequence::iterator vBegin = vertices->begin();
-    const SieveMesh::label_sequence::iterator vEnd = vertices->end();
 
-    fCoordinatesDim->setChart(topology::Mesh::RealSection::chart_type(faultSieveMesh->heightStratum(0)->size(),
-								   faultSieveMesh->getSieve()->getChart().max()));
+    fCoordinatesDim->setChart(topology::Mesh::RealSection::chart_type(fvMin,
+								      fvMax+1));
     for (SieveMesh::label_sequence::iterator v_iter = vBegin;
 	 v_iter != vEnd;
 	 ++v_iter) {
@@ -853,9 +859,9 @@
 	++v_iter) {
       if (fRenumbering.find(*v_iter) == fRenumberingEnd)
 	continue;
-    fCoordinatesDim->updatePoint(fRenumbering[*v_iter], 
-			      coordinatesDim->restrictPoint(*v_iter));
-    }
+      fCoordinatesDim->updatePoint(fRenumbering[*v_iter], 
+				   coordinatesDim->restrictPoint(*v_iter));
+    } // for
     //faultSieveMesh->view("Parallel fault mesh");
   } // if
 

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/faults/CohesiveTopology.hh
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/faults/CohesiveTopology.hh	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/faults/CohesiveTopology.hh	2010-04-10 00:35:22 UTC (rev 16524)
@@ -39,8 +39,10 @@
    * @param faultMesh Finite-element mesh of fault (output).
    * @param faultBoundary Finite-element mesh of fault boundary (output).
    * @param mesh Finite-element mesh of domain.
-   * @param faultVertices Vertices assocated with faces of cells defining 
-   *   fault surface
+   * @param groupdField Group of vertices assocated with faces of
+   *   cells defining fault surface
+   * @param flipFault Flag indicating to flip positive/negative sides
+   * of the fault.
    */
   static
   void createFault(topology::SubMesh* faultMesh,
@@ -70,6 +72,7 @@
               const ALE::Obj<topology::Mesh::IntSection>& groupField,
               const int materialId,
               int& firstFaultVertex,
+              int& firstLagrangeVertex,
               int& firstFaultCell,
               const bool constraintCell = false);
 

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/faults/Fault.hh
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/faults/Fault.hh	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/faults/Fault.hh	2010-04-10 00:35:22 UTC (rev 16524)
@@ -92,6 +92,7 @@
   virtual
   void adjustTopology(topology::Mesh* const mesh,
                       int *firstFaultVertex,
+                      int *firstLagrangeVertex,
                       int *firstFaultCell,
                       const bool flipFault = false) = 0;
 

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/faults/FaultCohesive.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/faults/FaultCohesive.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/faults/FaultCohesive.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -106,6 +106,7 @@
 void
 pylith::faults::FaultCohesive::adjustTopology(topology::Mesh* const mesh,
                                               int *firstFaultVertex,
+                                              int *firstLagrangeVertex,
                                               int *firstFaultCell,
                                               const bool flipFault)
 { // adjustTopology
@@ -133,7 +134,7 @@
 				  flipFault);
 
     CohesiveTopology::create(mesh, faultMesh, faultBoundary, groupField, id(), 
-			     *firstFaultVertex, *firstFaultCell, useLagrangeConstraints());
+                             *firstFaultVertex, *firstLagrangeVertex, *firstFaultCell, useLagrangeConstraints());
 
   } else {
     //std::cout << "BEFORE ADJUSTING TOPOLOGY FOR FAULT '" << label() << "' firstFaultVertex: " << *firstFaultVertex << ", firstFaultCell: " << *firstFaultCell << std::endl;
@@ -155,7 +156,7 @@
       sieveMesh->getIntSection(label());
     assert(!groupField.isNull());
     CohesiveTopology::create(mesh, faultMesh, faultBoundary, groupField, id(),
-                  *firstFaultVertex, *firstFaultCell, useLagrangeConstraints());
+                             *firstFaultVertex, *firstLagrangeVertex, *firstFaultCell, useLagrangeConstraints());
 
     //std::cout << "AFTER ADJUSTING TOPOLOGY FOR FAULT '" << label() << "' firstFaultVertex: " << *firstFaultVertex << ", firstFaultCell: " << *firstFaultCell << std::endl;
   } // if/else
@@ -430,7 +431,7 @@
       } // for
     } // for
     areaVisitor.clear();
-    faultSieveMesh->updateAdd(*c_iter, areaVisitor);
+    faultSieveMesh->updateClosure(*c_iter, areaVisitor);
 
     PetscLogFlops( numQuadPts*(1+numBasis*2) );
   } // for

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/faults/FaultCohesive.hh
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/faults/FaultCohesive.hh	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/faults/FaultCohesive.hh	2010-04-10 00:35:22 UTC (rev 16524)
@@ -71,16 +71,19 @@
   /** Adjust mesh topology for fault implementation.
    *
    * If firstFaultVertex == 0, then firstFaultVertex is set to the first point
-   * not currently used in the mesh, and firstFaultCell is incremented with this
-   * point. These values are updated as new fault vertices and cells are added.
+   * not currently used in the mesh, and firstLagrangeVertex/firstFaultCell are
+   * incremented with this point. These values are updated as new fault vertices
+   * and cells are added.
    *
    * @param mesh PETSc mesh.
    * @param firstFaultVertex The first point eligible to become a new fault vertex
+   * @param firstLagrangeVertex The first point eligible to become a new Lagrange vertex
    * @param firstFaultCell The first point eligible to become a new fault cell
    * @param flipFault Flip fault orientation.
    */
   void adjustTopology(topology::Mesh* const mesh,
                       int *firstFaultVertex,
+                      int *firstLagrangeVertex,
                       int *firstFaultCell,
                       const bool flipFault = false);
 

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/faults/FaultCohesiveKin.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/faults/FaultCohesiveKin.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/faults/FaultCohesiveKin.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -414,7 +414,7 @@
 #endif
 
     residualVisitor.clear();
-    sieveMesh->updateAdd(*c_iter, residualVisitor);
+    sieveMesh->updateClosure(*c_iter, residualVisitor);
   } // for
 
   // FIX THIS

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/faults/TopologyOps.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/faults/TopologyOps.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/faults/TopologyOps.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -302,37 +302,55 @@
                                                    const ALE::Obj<ALE::Mesh::arrow_section_type>& orientation,
                                                    const ALE::Obj<ALE::Mesh>& fault)
 {
-  // Must check the orientation here
+  assert(!mesh.isNull());
+  assert(!orientation.isNull());
+  assert(!fault.isNull());
+
   typedef ALE::Selection<ALE::Mesh> selection;
-  const ALE::Obj<ALE::Mesh::sieve_type>& faultSieve      = fault->getSieve();
-  const SieveMesh::point_type            firstFaultCell  = *fault->heightStratum(1)->begin();
-  const ALE::Obj<ALE::Mesh::label_sequence>& fFaces      = fault->heightStratum(2);
-  const int                         numFaultFaces   = fFaces->size();
-  const int                         faultDepth      = fault->depth()-1; // Depth of fault cells
-  int                               numFaultCorners = 0; // The number of vertices in a fault cell
-  int                               faultFaceSize   = 0; // The number of vertices in a face between fault cells
-  int                               faceSize        = 0; // The number of vertices in a mesh face
-  const int                         debug           = fault->debug();
-  ALE::Obj<PointSet>                     newCells        = new PointSet();
-  ALE::Obj<PointSet>                     loopCells       = new PointSet();
-  PointSet                          flippedCells;   // Incorrectly oriented fault cells
-  PointSet                          facesSeen;      // Fault faces already considered
-  PointSet                          cellsSeen;      // Fault cells already matched
-  PointArray                        faceVertices;
 
+  // Must check the orientation here
+  const ALE::Obj<ALE::Mesh::sieve_type>& faultSieve = fault->getSieve();
+  assert(!faultSieve.isNull());
+  const SieveMesh::point_type firstFaultCell  = 
+    *fault->heightStratum(1)->begin();
+  const ALE::Obj<ALE::Mesh::label_sequence>& fFaces = fault->heightStratum(2);
+  assert(!fFaces.isNull());
+  const int numFaultFaces = fFaces->size();
+  const int faultDepth = fault->depth()-1; // Depth of fault cells
+  int numFaultCorners = 0; // The number of vertices in a fault cell
+  int faultFaceSize = 0; // The number of vertices in a face between fault cells
+  int faceSize = 0; // The number of vertices in a mesh face
+  const int debug = fault->debug();
+  ALE::Obj<PointSet> newCells = new PointSet();
+  assert(!newCells.isNull());
+  ALE::Obj<PointSet> loopCells = new PointSet();
+  assert(!loopCells.isNull());
+  PointSet flippedCells;   // Incorrectly oriented fault cells
+  PointSet facesSeen;      // Fault faces already considered
+  PointSet cellsSeen;      // Fault cells already matched
+  PointArray faceVertices;
+
   if (!fault->commRank()) {
     faceSize        = selection::numFaceVertices(mesh);
     numFaultCorners = faultSieve->nCone(firstFaultCell, faultDepth)->size();
     if (debug) std::cout << "  Fault corners " << numFaultCorners << std::endl;
-    if (dim == 0) {
+    if (dim > 0) {
+      assert(numFaultCorners == faceSize);
+    } else {
+      // dim is 0
       assert(numFaultCorners == faceSize-1);
-    } else {
-      assert(numFaultCorners == faceSize);
     }
     if (faultDepth == 1) {
       faultFaceSize = 1;
     } else {
-      faultFaceSize = faultSieve->nCone(*fFaces->begin(), faultDepth-1)->size();
+      if (dim > 0) {
+	assert(fFaces->size() > 0);
+	assert(faultDepth > 0);
+	faultFaceSize = faultSieve->nCone(*fFaces->begin(), faultDepth-1)->size();
+      } else {
+	// dim is 0
+	faultFaceSize = 1;
+      } // if/else
     }
   }
   if (debug) std::cout << "  Fault face size " << faultFaceSize << std::endl;

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/faults/TopologyVisitors.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/faults/TopologyVisitors.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/faults/TopologyVisitors.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -134,8 +134,15 @@
   } // if
   bool classified = false;
     
-  if (debug)
+  if (debug) {
     std::cout << "Checking neighbor " << point << std::endl;
+    ALE::ISieveVisitor::PointRetriever<Sieve> pR(9);
+
+    sieve.cone(point, pR);
+    for(int i = 0; i < pR.getSize(); ++i) {
+      std::cout << "  cone point " << pR.getPoints()[i] << std::endl;
+    }
+  }
   if (vReplaceCells.find(point) != vReplaceCells.end()) {
     if (debug) 
       std::cout << "  already in replaceCells" << std::endl;

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/feassemble/ElasticityExplicit.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/feassemble/ElasticityExplicit.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/feassemble/ElasticityExplicit.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -282,7 +282,7 @@
     // Assemble cell contribution into field
     _logger->eventBegin(updateEvent);
     residualVisitor.clear();
-    sieveMesh->updateAdd(*c_iter, residualVisitor);
+    sieveMesh->updateClosure(*c_iter, residualVisitor);
     _logger->eventEnd(updateEvent);
   } // for
 } // integrateResidual

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/feassemble/ElasticityImplicit.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/feassemble/ElasticityImplicit.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/feassemble/ElasticityImplicit.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -308,7 +308,7 @@
     // Assemble cell contribution into field
     _logger->eventBegin(updateEvent);
     residualVisitor.clear();
-    sieveMesh->updateAdd(*c_iter, residualVisitor);
+    sieveMesh->updateClosure(*c_iter, residualVisitor);
     // residualSection->view("After stress contribution");
     _logger->eventEnd(updateEvent);
   } // for

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/feassemble/IntegratorElasticity.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/feassemble/IntegratorElasticity.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/feassemble/IntegratorElasticity.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -199,11 +199,11 @@
   } else {
       std::cerr << "Bad cell dimension '" << cellDim << "'." << std::endl;
       assert(0);
-      throw std::logic_error("Bad cell dimension in IntegratorElasticity.");
+      throw std::logic_error("Bad cell dimension in "
+			     "IntegratorElasticity::updateStateVars().");
   } // else
 
   // Allocate arrays for cell data.
-  double_array dispCell(numBasis*spaceDim);
   double_array strainCell(numQuadPts*tensorSize);
   strainCell = 0.0;
 
@@ -218,6 +218,7 @@
   const SieveMesh::label_sequence::iterator cellsEnd = cells->end();
 
   // Get fields
+  double_array dispCell(numBasis*spaceDim);
   const topology::Field<topology::Mesh>& disp = fields->get("disp(t)");
   const ALE::Obj<RealSection>& dispSection = disp.section();
   assert(!dispSection.isNull());
@@ -247,6 +248,9 @@
     _quadrature->computeGeometry(coordinatesCell, *c_iter);
 #endif
 
+    // Get physical properties and state variables for cell.
+    _material->retrievePropsAndVars(*c_iter);
+
     // Restrict input fields to cell
     dispVisitor.clear();
     sieveMesh->restrictClosure(*c_iter, dispVisitor);
@@ -261,7 +265,7 @@
     // Update material state
     _material->updateStateVars(strainCell, *c_iter);
   } // for
-} // updateState
+} // updateStateVars
 
 // ----------------------------------------------------------------------
 // Verify configuration is acceptable.
@@ -625,15 +629,17 @@
 pylith::feassemble::IntegratorElasticity::_elasticityResidual1D(
 				     const double_array& stress)
 { // _elasticityResidual1D
+  const int spaceDim = 1;
+  const int cellDim = 1;
+
   const int numQuadPts = _quadrature->numQuadPts();
   const int numBasis = _quadrature->numBasis();
-  const int spaceDim = _quadrature->spaceDim();
-  const int cellDim = _quadrature->cellDim();
   const double_array& quadWts = _quadrature->quadWts();
   const double_array& jacobianDet = _quadrature->jacobianDet();
   const double_array& basisDeriv = _quadrature->basisDeriv();
 
-  assert(1 == cellDim);
+  assert(_quadrature->spaceDim() == spaceDim);
+  assert(_quadrature->cellDim() == cellDim);
   assert(quadWts.size() == numQuadPts);
 
   for (int iQuad=0; iQuad < numQuadPts; ++iQuad) {
@@ -653,30 +659,35 @@
 pylith::feassemble::IntegratorElasticity::_elasticityResidual2D(
 				     const double_array& stress)
 { // _elasticityResidual2D
+  const int cellDim = 2;
+  const int spaceDim = 2;
+  const int stressSize = 3;
+
   const int numQuadPts = _quadrature->numQuadPts();
   const int numBasis = _quadrature->numBasis();
-  const int spaceDim = _quadrature->spaceDim();
-  const int cellDim = _quadrature->cellDim();
   const double_array& quadWts = _quadrature->quadWts();
   const double_array& jacobianDet = _quadrature->jacobianDet();
   const double_array& basisDeriv = _quadrature->basisDeriv();
   
-  assert(2 == cellDim);
+  assert(_quadrature->spaceDim() == spaceDim);
+  assert(_quadrature->cellDim() == cellDim);
   assert(quadWts.size() == numQuadPts);
-  const int stressSize = 3;
 
   for (int iQuad=0; iQuad < numQuadPts; ++iQuad) {
+    const int iQs = iQuad*stressSize;
     const double wt = quadWts[iQuad] * jacobianDet[iQuad];
-    const double s11 = stress[iQuad*stressSize+0];
-    const double s22 = stress[iQuad*stressSize+1];
-    const double s12 = stress[iQuad*stressSize+2];
+    const double s11 = stress[iQs  ];
+    const double s22 = stress[iQs+1];
+    const double s12 = stress[iQs+2];
     for (int iBasis=0, iQ=iQuad*numBasis*spaceDim;
 	 iBasis < numBasis;
 	 ++iBasis) {
-      const double N1 = wt*basisDeriv[iQ+iBasis*spaceDim  ];
-      const double N2 = wt*basisDeriv[iQ+iBasis*spaceDim+1];
-      _cellVector[iBasis*spaceDim  ] -= N1*s11 + N2*s12;
-      _cellVector[iBasis*spaceDim+1] -= N1*s12 + N2*s22;
+      const int iBlock = iBasis*spaceDim;
+      const double N1 = wt*basisDeriv[iQ+iBlock  ];
+      const double N2 = wt*basisDeriv[iQ+iBlock+1];
+
+      _cellVector[iBlock  ] -= N1*s11 + N2*s12;
+      _cellVector[iBlock+1] -= N1*s12 + N2*s22;
     } // for
   } // for
   PetscLogFlops(numQuadPts*(1+numBasis*(8+2+9)));
@@ -688,34 +699,37 @@
 pylith::feassemble::IntegratorElasticity::_elasticityResidual3D(
 				     const double_array& stress)
 { // _elasticityResidual3D
+  const int spaceDim = 3;
+  const int cellDim = 3;
+  const int stressSize = 6;
+
   const int numQuadPts = _quadrature->numQuadPts();
   const int numBasis = _quadrature->numBasis();
-  const int spaceDim = _quadrature->spaceDim();
-  const int cellDim = _quadrature->cellDim();
   const double_array& quadWts = _quadrature->quadWts();
   const double_array& jacobianDet = _quadrature->jacobianDet();
   const double_array& basisDeriv = _quadrature->basisDeriv();
   
-  assert(3 == cellDim);
+  assert(_quadrature->spaceDim() == spaceDim);
+  assert(_quadrature->cellDim() == cellDim);
   assert(quadWts.size() == numQuadPts);
-  const int stressSize = 6;
   
   for (int iQuad=0; iQuad < numQuadPts; ++iQuad) {
+    const int iQs = iQuad * stressSize;
     const double wt = quadWts[iQuad] * jacobianDet[iQuad];
-    const double s11 = stress[iQuad*stressSize+0];
-    const double s22 = stress[iQuad*stressSize+1];
-    const double s33 = stress[iQuad*stressSize+2];
-    const double s12 = stress[iQuad*stressSize+3];
-    const double s23 = stress[iQuad*stressSize+4];
-    const double s13 = stress[iQuad*stressSize+5];
+    const double s11 = stress[iQs  ];
+    const double s22 = stress[iQs+1];
+    const double s33 = stress[iQs+2];
+    const double s12 = stress[iQs+3];
+    const double s23 = stress[iQs+4];
+    const double s13 = stress[iQs+5];
     
     for (int iBasis=0, iQ=iQuad*numBasis*spaceDim;
-	 iBasis < numBasis;
-	 ++iBasis) {
+        iBasis < numBasis;
+        ++iBasis) {
       const int iBlock = iBasis*spaceDim;
-      const double N1 = wt*basisDeriv[iQ+iBasis*spaceDim+0];
-      const double N2 = wt*basisDeriv[iQ+iBasis*spaceDim+1];
-      const double N3 = wt*basisDeriv[iQ+iBasis*spaceDim+2];
+      const double N1 = wt*basisDeriv[iQ+iBlock+0];
+      const double N2 = wt*basisDeriv[iQ+iBlock+1];
+      const double N3 = wt*basisDeriv[iQ+iBlock+2];
 
       _cellVector[iBlock  ] -= N1*s11 + N2*s12 + N3*s13;
       _cellVector[iBlock+1] -= N1*s12 + N2*s22 + N3*s23;
@@ -731,15 +745,17 @@
 pylith::feassemble::IntegratorElasticity::_elasticityJacobian1D(
 			       const double_array& elasticConsts)
 { // _elasticityJacobian1D
+  const int cellDim = 1;
+  const int spaceDim = 1;
+
   const int numQuadPts = _quadrature->numQuadPts();
   const int numBasis = _quadrature->numBasis();
-  const int spaceDim = _quadrature->spaceDim();
-  const int cellDim = _quadrature->cellDim();
   const double_array& quadWts = _quadrature->quadWts();
   const double_array& jacobianDet = _quadrature->jacobianDet();
   const double_array& basisDeriv = _quadrature->basisDeriv();
   
-  assert(1 == cellDim);
+  assert(_quadrature->spaceDim() == spaceDim);
+  assert(_quadrature->cellDim() == cellDim);
   assert(quadWts.size() == numQuadPts);
   
   for (int iQuad=0; iQuad < numQuadPts; ++iQuad) {
@@ -764,17 +780,19 @@
 pylith::feassemble::IntegratorElasticity::_elasticityJacobian2D(
 			       const double_array& elasticConsts)
 { // _elasticityJacobian2D
+  const int spaceDim = 2;
+  const int cellDim = 2;
+  const int numConsts = 6;
+
   const int numQuadPts = _quadrature->numQuadPts();
   const int numBasis = _quadrature->numBasis();
-  const int spaceDim = _quadrature->spaceDim();
-  const int cellDim = _quadrature->cellDim();
   const double_array& quadWts = _quadrature->quadWts();
   const double_array& jacobianDet = _quadrature->jacobianDet();
   const double_array& basisDeriv = _quadrature->basisDeriv();
   
-  assert(2 == cellDim);
+  assert(_quadrature->spaceDim() == spaceDim);
+  assert(_quadrature->cellDim() == cellDim);
   assert(quadWts.size() == numQuadPts);
-  const int numConsts = 6;
 
   for (int iQuad=0; iQuad < numQuadPts; ++iQuad) {
     const double wt = quadWts[iQuad] * jacobianDet[iQuad];
@@ -828,17 +846,19 @@
 pylith::feassemble::IntegratorElasticity::_elasticityJacobian3D(
 			       const double_array& elasticConsts)
 { // _elasticityJacobian3D
+  const int spaceDim = 3;
+  const int cellDim = 3;
+  const int numConsts = 21;
+
   const int numQuadPts = _quadrature->numQuadPts();
   const int numBasis = _quadrature->numBasis();
-  const int spaceDim = _quadrature->spaceDim();
-  const int cellDim = _quadrature->cellDim();
   const double_array& quadWts = _quadrature->quadWts();
   const double_array& jacobianDet = _quadrature->jacobianDet();
   const double_array& basisDeriv = _quadrature->basisDeriv();
   
-  assert(3 == cellDim);
+  assert(_quadrature->spaceDim() == spaceDim);
+  assert(_quadrature->cellDim() == cellDim);
   assert(quadWts.size() == numQuadPts);
-  const int numConsts = 21;
 
   // Compute Jacobian for consistent tangent matrix
   for (int iQuad=0; iQuad < numQuadPts; ++iQuad) {
@@ -970,12 +990,12 @@
   assert(0 != strain);
   
   const int dim = 2;
+  const int strainSize = 3;
 
   assert(basisDeriv.size() == numQuadPts*numBasis*dim);
   assert(disp.size() == numBasis*dim);
 
   (*strain) = 0.0;
-  const int strainSize = 3;
   for (int iQuad=0; iQuad < numQuadPts; ++iQuad)
     for (int iBasis=0, iQ=iQuad*numBasis*dim; iBasis < numBasis; ++iBasis) {
       (*strain)[iQuad*strainSize+0] += 
@@ -1000,12 +1020,12 @@
   assert(0 != strain);
 
   const int dim = 3;
+  const int strainSize = 6;
 
   assert(basisDeriv.size() == numQuadPts*numBasis*dim);
   assert(disp.size() == numBasis*dim);
 
   (*strain) = 0.0;
-  const int strainSize = 6;
   for (int iQuad=0; iQuad < numQuadPts; ++iQuad)
     for (int iBasis=0, iQ=iQuad*numBasis*dim; iBasis < numBasis; ++iBasis) {
       (*strain)[iQuad*strainSize+0] += 

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/materials/ElasticMaterial.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/materials/ElasticMaterial.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/materials/ElasticMaterial.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -226,11 +226,6 @@
   assert(_initialStrainCell.size() == numQuadPts*_tensorSize);
   assert(totalStrain.size() == numQuadPts*_tensorSize);
 
-  const ALE::Obj<RealSection>& stateVarsSection = _stateVars->section();
-  assert(!stateVarsSection.isNull());
-  stateVarsSection->restrictPoint(cell, &_stateVarsCell[0],
-				  _stateVarsCell.size());
-  
   for (int iQuad=0; iQuad < numQuadPts; ++iQuad)
     _updateStateVars(&_stateVarsCell[iQuad*numVarsQuadPt], numVarsQuadPt,
 		     &_propertiesCell[iQuad*numPropsQuadPt], 
@@ -239,6 +234,8 @@
 		     &_initialStressCell[iQuad*_tensorSize], _tensorSize,
 		     &_initialStrainCell[iQuad*_tensorSize], _tensorSize);
   
+  const ALE::Obj<RealSection>& stateVarsSection = _stateVars->section();
+  assert(!stateVarsSection.isNull());  
   stateVarsSection->updatePoint(cell, &_stateVarsCell[0]);
 } // updateStateVars
 

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/meshio/CellFilterAvg.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/meshio/CellFilterAvg.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/meshio/CellFilterAvg.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -93,10 +93,14 @@
   
   const ALE::Obj<SieveMesh>& sieveMesh = fieldIn.mesh().sieveMesh();
   assert(!sieveMesh.isNull());
+  const int cellDepth = (sieveMesh->depth() == -1) ? -1 : 1;
+  const int depth = (0 == label) ? cellDepth : labelId;
+  const std::string labelName = (0 == label) ?
+    ((sieveMesh->hasLabel("censored depth")) ?
+     "censored depth" : "depth") : label;
 
-  const ALE::Obj<label_sequence>& cells = (0 == label) ?
-    sieveMesh->heightStratum(0) :
-    sieveMesh->getLabelStratum(label, labelId);
+  const ALE::Obj<label_sequence>& cells = 
+    sieveMesh->getLabelStratum(labelName, depth);
   assert(!cells.isNull());
   const typename label_sequence::iterator cellsBegin = cells->begin();
   const typename label_sequence::iterator cellsEnd = cells->end();
@@ -168,6 +172,7 @@
        c_iter != cellsEnd;
        ++c_iter) {
     const double* values = sectionIn->restrictPoint(*c_iter);
+    assert(totalFiberDim == sectionIn->getFiberDimension(*c_iter));
     
     fieldAvgCell = 0.0;
     for (int iQuad=0; iQuad < numQuadPts; ++iQuad)

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/meshio/DataWriterVTK.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/meshio/DataWriterVTK.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/meshio/DataWriterVTK.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -27,7 +27,8 @@
   _timeFormat("%f"),
   _viewer(0),
   _wroteVertexHeader(false),
-  _wroteCellHeader(false)
+  _wroteCellHeader(false),
+  _precision(6)
 { // constructor
 } // constructor
 
@@ -40,6 +41,21 @@
 } // destructor  
 
 // ----------------------------------------------------------------------
+// Set precision of floating point values in output.
+template<typename mesh_type, typename field_type>
+void
+pylith::meshio::DataWriterVTK<mesh_type,field_type>::precision(const int value)
+{ // precision
+  if (value <= 0) {
+    std::ostringstream msg;
+    msg << "Floating point precision (" << value << ") must be positive.";
+    throw std::runtime_error(msg.str());
+  } // if
+
+  _precision = value;
+} // precision
+
+// ----------------------------------------------------------------------
 // Deallocate PETSc and local data structures.
 template<typename mesh_type, typename field_type>
 void
@@ -201,7 +217,7 @@
     } // if
 
     err = VTKViewer::writeField(section, field.label(), fiberDim, numbering,
-				_viewer, enforceDim);
+				_viewer, enforceDim, _precision);
     CHECK_PETSC_ERROR(err);
   } catch (const std::exception& err) {
     std::ostringstream msg;
@@ -269,7 +285,7 @@
     } // if
 
     VTKViewer::writeField(section, field.label(), fiberDim, numbering,
-			  _viewer, enforceDim);
+			  _viewer, enforceDim, _precision);
   } catch (const std::exception& err) {
     std::ostringstream msg;
     msg << "Error while writing field '" << field.label() << "' at time " 

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/meshio/DataWriterVTK.hh
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/meshio/DataWriterVTK.hh	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/meshio/DataWriterVTK.hh	2010-04-10 00:35:22 UTC (rev 16524)
@@ -70,6 +70,12 @@
    */
   void timeConstant(const double value);
 
+  /** Set precision of floating point values in output.
+   *
+   * @param value Precision for floating point values.
+   */
+  void precision(const int value);
+
   /** Prepare file for data at a new time step.
    *
    * @param t Time stamp for new data
@@ -140,6 +146,8 @@
 
   PetscViewer _viewer; ///< Output file
 
+  int _precision; ///< Precision of floating point values in output.
+
   bool _wroteVertexHeader; ///< True if wrote header for vertex data.
   bool _wroteCellHeader; ///< True if wrote header for cell data
 

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/meshio/MeshIO.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/meshio/MeshIO.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/meshio/MeshIO.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -18,7 +18,6 @@
 #include "pylith/utils/array.hh" // USES double_array, int_array
 
 #include "spatialdata/geocoords/CoordSys.hh" // USES CoordSys
-#include "spatialdata/units/Nondimensional.hh" // USES Nondimensional
 
 #include "Selection.hh" // USES boundary()
 
@@ -290,6 +289,13 @@
   const ALE::Obj<SieveMesh>& sieveMesh = _mesh->sieveMesh();
   assert(!sieveMesh.isNull());
 
+  if (sieveMesh->hasIntSection(name)) {
+    std::ostringstream msg;
+    msg << "Could not setup group '" << name
+	<< "'. Group already exists in mesh.";
+    throw std::runtime_error(msg.str());
+  } // if
+
   ALE::MemoryLogger& logger = ALE::MemoryLogger::singleton();
   logger.stagePush("VertexGroups");
   const ALE::Obj<IntSection>& groupField = sieveMesh->getIntSection(name);
@@ -396,6 +402,13 @@
   const ALE::Obj<SieveMesh>& sieveMesh = _mesh->sieveMesh();
   assert(!sieveMesh.isNull());
 
+  if (!sieveMesh->hasIntSection(name)) {
+    std::ostringstream msg;
+    msg << "Could not get group '" << name
+	<< "'. Group missing from mesh.";
+    throw std::runtime_error(msg.str());
+  } // if
+
   const ALE::Obj<IntSection>& groupField = sieveMesh->getIntSection(name);
   assert(!groupField.isNull());
   const IntSection::chart_type& chart = groupField->getChart();

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/topology/Mesh.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/topology/Mesh.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/topology/Mesh.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -11,6 +11,7 @@
 //
 
 #include <portinfo>
+#include <stdexcept>
 
 #include "Mesh.hh" // implementation of class methods
 
@@ -18,6 +19,10 @@
 #include "spatialdata/units/Nondimensional.hh" // USES Nondimensional
 #include "pylith/utils/array.hh" // USES double_array
 
+#include <stdexcept> // USES std::runtime_error
+#include <sstream> // USES std::ostringstream
+#include <cassert> // USES assert()
+
 // ----------------------------------------------------------------------
 // Default constructor
 pylith::topology::Mesh::Mesh(void) :
@@ -116,15 +121,21 @@
       ++v_iter) {
     coordsSection->restrictPoint(*v_iter,
 			       &coordsVertex[0], coordsVertex.size());
-    coordsDimSection->restrictPoint(*v_iter, &coordsDimVertex[0],
-				  coordsDimVertex.size());
-    for (int iDim=0; iDim < spaceDim; ++iDim)
-      coordsDimVertex[iDim] = coordsVertex[iDim];
 
+    // Save dimensioned coordinates in coordsDimVertex
+    coordsDimVertex = coordsVertex;
+
     // Nondimensionalize original coordinates.
     normalizer.nondimensionalize(&coordsVertex[0], spaceDim, lengthScale);
 
+    // Update section with nondimensional coordinates
+    assert(coordsVertex.size() == 
+	   coordsSection->getFiberDimension(*v_iter));
     coordsSection->updatePoint(*v_iter, &coordsVertex[0]);
+    
+    // Update section with dimensioned coordinates
+    assert(coordsDimVertex.size() == 
+	   coordsDimSection->getFiberDimension(*v_iter));
     coordsDimSection->updatePoint(*v_iter, &coordsDimVertex[0]);
   } // for
 } // nondimensionalize
@@ -153,16 +164,25 @@
 int
 pylith::topology::Mesh::groupSize(const char *name)
 { // groupSize
+  if (!_mesh->hasIntSection(name)) {
+    std::ostringstream msg;
+    msg << "Cannot get size of group '" << name
+	<< "'. Group missing from mesh.";
+    throw std::runtime_error(msg.str());
+  } // if
+
   const ALE::Obj<IntSection>&            group    = _mesh->getIntSection(name);
   const IntSection::chart_type&          chart    = group->getChart();
   IntSection::chart_type::const_iterator chartEnd = chart.end();
   int                                    size     = 0;
 
-  for(IntSection::chart_type::const_iterator c_iter = chart.begin(); c_iter != chartEnd; ++c_iter) {
-    if (group->getFiberDimension(*c_iter)) {
+  for(IntSection::chart_type::const_iterator c_iter = chart.begin();
+      c_iter != chartEnd;
+      ++c_iter) {
+    if (group->getFiberDimension(*c_iter))
       size++;
-    }
-  }
+  } // for
+
   return size;
 } // groupSize
 

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/topology/SubMesh.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/topology/SubMesh.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/topology/SubMesh.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -11,6 +11,7 @@
 //
 
 #include <portinfo>
+#include <stdexcept>
 
 #include "SubMesh.hh" // implementation of class methods
 
@@ -18,6 +19,10 @@
 
 #include <Selection.hh> // USES ALE::Selection
 
+#include <stdexcept> // USES std::runtime_error
+#include <sstream> // USES std::ostringstream
+#include <cassert> // USES assert()
+
 // ----------------------------------------------------------------------
 // Default constructor
 pylith::topology::SubMesh::SubMesh(void) :
@@ -62,12 +67,14 @@
   const ALE::Obj<DomainSieveMesh>& meshSieveMesh = mesh.sieveMesh();
   assert(!meshSieveMesh.isNull());
 
-  const ALE::Obj<IntSection>& groupField = meshSieveMesh->getIntSection(label);
-  if (groupField.isNull()) {
+  if (!meshSieveMesh->hasIntSection(label)) {
     std::ostringstream msg;
     msg << "Could not find group of points '" << label << "' in mesh.";
     throw std::runtime_error(msg.str());
   } // if
+
+  const ALE::Obj<IntSection>& groupField = meshSieveMesh->getIntSection(label);
+  assert(!groupField.isNull());
   _mesh = 
     ALE::Selection<DomainSieveMesh>::submeshV<SieveMesh>(meshSieveMesh,
 							 groupField);

Modified: short/3D/PyLith/tags/v1.4.3/libsrc/utils/EventLogger.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/libsrc/utils/EventLogger.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/libsrc/utils/EventLogger.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -41,7 +41,7 @@
 			   "initializaing EventLogger.");
   
   _events.clear();
-  PetscErrorCode err = PetscCookieRegister(_className.c_str(), &_classId);
+  PetscErrorCode err = PetscClassIdRegister(_className.c_str(), &_classId);
   if (err) {
     std::ostringstream msg;
     msg << "Could not register logging class '" << _className << "'.";

Modified: short/3D/PyLith/tags/v1.4.3/modulesrc/faults/Fault.i
===================================================================
--- short/3D/PyLith/tags/v1.4.3/modulesrc/faults/Fault.i	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/modulesrc/faults/Fault.i	2010-04-10 00:35:22 UTC (rev 16524)
@@ -75,6 +75,7 @@
       virtual
       void adjustTopology(pylith::topology::Mesh* const mesh,
                           int *firstFaultVertex,
+                          int *firstLagrangeVertex,
                           int *firstFaultCell,
                           const bool flipFault = false) = 0;
       %clear int *firstFaultVertex, int *firstFaultCell;

Modified: short/3D/PyLith/tags/v1.4.3/modulesrc/faults/FaultCohesive.i
===================================================================
--- short/3D/PyLith/tags/v1.4.3/modulesrc/faults/FaultCohesive.i	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/modulesrc/faults/FaultCohesive.i	2010-04-10 00:35:22 UTC (rev 16524)
@@ -59,15 +59,25 @@
 
       /** Adjust mesh topology for fault implementation.
        *
+       * If firstFaultVertex == 0, then firstFaultVertex is set to the
+       * first point not currently used in the mesh, and
+       * firstLagrangeVertex/firstFaultCell are incremented with this
+       * point. These values are updated as new fault vertices and
+       * cells are added.
+       *
        * @param mesh PETSc mesh.
+       * @param firstFaultVertex The first point eligible to become a new fault vertex
+       * @param firstLagrangeVertex The first point eligible to become a new Lagrange vertex
+       * @param firstFaultCell The first point eligible to become a new fault cell
        * @param flipFault Flip fault orientation.
        */
-      %apply int *INOUT {int *firstFaultVertex, int *firstFaultCell};
+      %apply int *INOUT {int *firstFaultVertex, int *firstLagrangeVertex, int *firstFaultCell};
       void adjustTopology(pylith::topology::Mesh* const mesh,
                           int *firstFaultVertex,
+                          int *firstLagrangeVertex,
                           int *firstFaultCell,
                           const bool flipFault = false);
-      %clear int *firstFaultVertex, int *firstFaultCell;
+      %clear int *firstFaultVertex, int *firstLagrangeVertex, int *firstFaultCell;
       
       /** Cohesive cells use Lagrange multiplier constraints?
        *

Modified: short/3D/PyLith/tags/v1.4.3/modulesrc/meshio/DataWriterVTK.i
===================================================================
--- short/3D/PyLith/tags/v1.4.3/modulesrc/meshio/DataWriterVTK.i	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/modulesrc/meshio/DataWriterVTK.i	2010-04-10 00:35:22 UTC (rev 16524)
@@ -63,6 +63,12 @@
        */
       void timeConstant(const double value);
       
+      /** Set precision of floating point values in output.
+       *
+       * @param value Precision for floating point values.
+       */
+      void precision(const int value);
+
       /** Prepare file for data at a new time step.
        *
        * @param t Time stamp for new data

Modified: short/3D/PyLith/tags/v1.4.3/pylith/faults/LiuCosSlipFn.py
===================================================================
--- short/3D/PyLith/tags/v1.4.3/pylith/faults/LiuCosSlipFn.py	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/pylith/faults/LiuCosSlipFn.py	2010-04-10 00:35:22 UTC (rev 16524)
@@ -10,7 +10,7 @@
 # ----------------------------------------------------------------------
 #
 
-## @file pylith/faults/BruneSlipFn.py
+## @file pylith/faults/LiuCosSlipFn.py
 ##
 ## @brief Sine/cosine slip time function from Liu, Archuleta, and Hartzell,
 ## BSSA, 2006 (doi:10.1785/0120060036) which has a rapid rise and then
@@ -62,13 +62,13 @@
 
   # PUBLIC METHODS /////////////////////////////////////////////////////
 
-  def __init__(self, name="bruneslipfn"):
+  def __init__(self, name="liucosslipfn"):
     """
     Constructor.
     """
     SlipTimeFn.__init__(self, name)
     ModuleLiuCosSlipFn.__init__(self)
-    self._loggingPrefix = "BrSF "
+    self._loggingPrefix = "LCSF "
     return
 
 

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/bc/TestBoundaryMesh.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/bc/TestBoundaryMesh.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/bc/TestBoundaryMesh.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -129,11 +129,12 @@
 
   // Adjust topology
   faults::FaultCohesiveKin fault;
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh.sieveMesh()->getIntSection(_data->faultLabel)->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh.sieveMesh()->getIntSection(_data->faultLabel)->size();
+  int firstFaultCell      = mesh.sieveMesh()->getIntSection(_data->faultLabel)->size();
   fault.label(_data->faultLabel);
   fault.id(_data->faultId);
-  fault.adjustTopology(&mesh, &firstFaultVertex, &firstFaultCell, _flipFault);
+  fault.adjustTopology(&mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
 
   // Create submesh
   topology::SubMesh submesh(mesh, _data->bcLabel);

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestBruneSlipFn.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestBruneSlipFn.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestBruneSlipFn.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -398,8 +398,9 @@
   mesh->coordsys(&cs);
 
   // Create fault mesh
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh->sieveMesh()->getIntSection(faultLabel)->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh->sieveMesh()->getIntSection(faultLabel)->size();
+  int firstFaultCell      = mesh->sieveMesh()->getIntSection(faultLabel)->size();
   const bool useLagrangeConstraints = true;
   if (useLagrangeConstraints) {
     firstFaultCell += mesh->sieveMesh()->getIntSection(faultLabel)->size();
@@ -412,7 +413,7 @@
   CohesiveTopology::create(mesh, *faultMesh, faultBoundary, 
                            sieveMesh->getIntSection(faultLabel),
                            faultId,
-                           firstFaultVertex, firstFaultCell,
+                           firstFaultVertex, firstLagrangeVertex, firstFaultCell,
                            useLagrangeConstraints);
   // Need to copy coordinates from mesh to fault mesh since we are not
   // using create() instead of createParallel().
@@ -471,8 +472,9 @@
 
   // Create fault mesh
   topology::SubMesh faultMesh;
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
+  int firstFaultCell      = mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
   const bool useLagrangeConstraints = true;
   if (useLagrangeConstraints) {
     firstFaultCell += mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
@@ -485,7 +487,7 @@
   CohesiveTopology::create(&mesh, faultMesh, faultBoundary, 
                            sieveMesh->getIntSection(data.faultLabel),
                            data.faultId,
-                           firstFaultVertex, firstFaultCell,
+                           firstFaultVertex, firstLagrangeVertex, firstFaultCell,
                            useLagrangeConstraints);
   // Need to copy coordinates from mesh to fault mesh since we are not
   // using create() instead of createParallel().

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestConstRateSlipFn.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestConstRateSlipFn.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestConstRateSlipFn.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -317,8 +317,9 @@
   mesh->coordsys(&cs);
 
   // Create fault mesh
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh->sieveMesh()->getIntSection(faultLabel)->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh->sieveMesh()->getIntSection(faultLabel)->size();
+  int firstFaultCell      = mesh->sieveMesh()->getIntSection(faultLabel)->size();
   const bool useLagrangeConstraints = true;
   if (useLagrangeConstraints) {
     firstFaultCell += mesh->sieveMesh()->getIntSection(faultLabel)->size();
@@ -331,7 +332,7 @@
   CohesiveTopology::create(mesh, *faultMesh, faultBoundary, 
                            sieveMesh->getIntSection(faultLabel),
                            faultId,
-                           firstFaultVertex, firstFaultCell,
+                           firstFaultVertex, firstLagrangeVertex, firstFaultCell,
                            useLagrangeConstraints);
   // Need to copy coordinates from mesh to fault mesh since we are not
   // using create() instead of createParallel().
@@ -384,8 +385,9 @@
 
   // Create fault mesh
   topology::SubMesh faultMesh;
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
+  int firstFaultCell      = mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
   const bool useLagrangeConstraints = true;
   if (useLagrangeConstraints) {
     firstFaultCell += mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
@@ -398,7 +400,7 @@
   CohesiveTopology::create(&mesh, faultMesh, faultBoundary, 
                            sieveMesh->getIntSection(data.faultLabel),
                            data.faultId,
-                           firstFaultVertex, firstFaultCell,
+                           firstFaultVertex, firstLagrangeVertex, firstFaultCell,
                            useLagrangeConstraints);
   // Need to copy coordinates from mesh to fault mesh since we are not
   // using create() instead of createParallel().

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestEqKinSrc.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestEqKinSrc.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestEqKinSrc.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -237,8 +237,9 @@
   mesh->coordsys(&cs);
 
   // Create fault mesh
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh->sieveMesh()->getIntSection(faultLabel)->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh->sieveMesh()->getIntSection(faultLabel)->size();
+  int firstFaultCell      = mesh->sieveMesh()->getIntSection(faultLabel)->size();
   const bool useLagrangeConstraints = true;
   if (useLagrangeConstraints) {
     firstFaultCell += mesh->sieveMesh()->getIntSection(faultLabel)->size();
@@ -251,7 +252,7 @@
   CohesiveTopology::create(mesh, *faultMesh, faultBoundary, 
                            sieveMesh->getIntSection(faultLabel),
                            faultId,
-                           firstFaultVertex, firstFaultCell,
+                           firstFaultVertex, firstLagrangeVertex, firstFaultCell,
                            useLagrangeConstraints);
   // Need to copy coordinates from mesh to fault mesh since we are not
   // using create() instead of createParallel().

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestFaultCohesive.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestFaultCohesive.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestFaultCohesive.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -21,6 +21,9 @@
 #include "pylith/utils/array.hh" // USES int_array, double_array
 #include "pylith/meshio/MeshIOAscii.hh" // USES MeshIOAscii
 
+#include "spatialdata/units/Nondimensional.hh" // USES Nondimensional
+#include "spatialdata/geocoords/CSCart.hh" // USES CSCart
+
 #include "data/CohesiveDataLine2.hh" // USES CohesiveDataLine2
 #include "data/CohesiveDataTri3.hh" // USES CohesiveDataTri3
 #include "data/CohesiveDataTri3b.hh" // USES CohesiveDataTri3b
@@ -491,15 +494,23 @@
   iohandler.interpolate(false);
   iohandler.read(&mesh);
 
+  spatialdata::geocoords::CSCart cs;
+  spatialdata::units::Nondimensional normalizer;
+  cs.setSpaceDim(mesh.dimension());
+  cs.initialize();
+  mesh.coordsys(&cs);
+  mesh.nondimensionalize(normalizer);
+  
   CPPUNIT_ASSERT(0 != fault);
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh.sieveMesh()->getIntSection("fault")->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh.sieveMesh()->getIntSection("fault")->size();
+  int firstFaultCell      = mesh.sieveMesh()->getIntSection("fault")->size();
   if (dynamic_cast<FaultCohesive*>(fault)->useLagrangeConstraints()) {
     firstFaultCell += mesh.sieveMesh()->getIntSection("fault")->size();
   }
   fault->id(1);
   fault->label("fault");
-  fault->adjustTopology(&mesh, &firstFaultVertex, &firstFaultCell, flipFault);
+  fault->adjustTopology(&mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, flipFault);
   //mesh->view(data.filename);
 
   CPPUNIT_ASSERT_EQUAL(data.cellDim, mesh.dimension());
@@ -626,22 +637,24 @@
   iohandler.read(&mesh);
 
   CPPUNIT_ASSERT(0 != faultA);
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh.sieveMesh()->getIntSection("faultA")->size() + mesh.sieveMesh()->getIntSection("faultB")->size();
+  CPPUNIT_ASSERT(0 != faultB);
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh.sieveMesh()->getIntSection("faultA")->size() + mesh.sieveMesh()->getIntSection("faultB")->size();
+  int firstFaultCell      = mesh.sieveMesh()->getIntSection("faultA")->size() + mesh.sieveMesh()->getIntSection("faultB")->size();
   if (dynamic_cast<FaultCohesive*>(faultA)->useLagrangeConstraints()) {
     firstFaultCell += mesh.sieveMesh()->getIntSection("faultA")->size();
   }
+  if (dynamic_cast<FaultCohesive*>(faultB)->useLagrangeConstraints()) {
+    firstFaultCell += mesh.sieveMesh()->getIntSection("faultB")->size();
+  }
+
   faultA->id(1);
   faultA->label("faultA");
-  faultA->adjustTopology(&mesh, &firstFaultVertex, &firstFaultCell, flipFaultA);
+  faultA->adjustTopology(&mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, flipFaultA);
 
-  CPPUNIT_ASSERT(0 != faultB);
-  if (dynamic_cast<FaultCohesive*>(faultB)->useLagrangeConstraints()) {
-    firstFaultCell += mesh.sieveMesh()->getIntSection("faultB")->size();
-  }
   faultB->id(2);
   faultB->label("faultB");
-  faultB->adjustTopology(&mesh, &firstFaultVertex, &firstFaultCell, flipFaultB);
+  faultB->adjustTopology(&mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, flipFaultB);
 
   //sieveMesh->view(data.filename);
   CPPUNIT_ASSERT_EQUAL(data.cellDim, mesh.dimension());

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestFaultCohesiveKin.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestFaultCohesiveKin.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestFaultCohesiveKin.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -263,10 +263,7 @@
       for (int i=0; i < fiberDimE; ++i) {
 	const int index = iVertex*spaceDim+i;
 	const double valE = valsE[index];
-	if (fabs(valE) > tolerance)
-	  CPPUNIT_ASSERT_DOUBLES_EQUAL(1.0, vals[i]/valE, tolerance);
-	else
-	  CPPUNIT_ASSERT_DOUBLES_EQUAL(valE, vals[i], tolerance);
+	CPPUNIT_ASSERT_DOUBLES_EQUAL(valE, vals[i], tolerance);
       } // for
     } // for
   } // Integrate residual with disp (as opposed to disp increment).
@@ -294,10 +291,7 @@
       for (int i=0; i < fiberDimE; ++i) {
 	const int index = iVertex*spaceDim+i;
 	const double valE = valsE[index];
-	if (fabs(valE) > tolerance)
-	  CPPUNIT_ASSERT_DOUBLES_EQUAL(1.0, vals[i]/valE, tolerance);
-	else
-	  CPPUNIT_ASSERT_DOUBLES_EQUAL(valE, vals[i], tolerance);
+	CPPUNIT_ASSERT_DOUBLES_EQUAL(valE, vals[i], tolerance);
       } // for
     } // for
   } // Integrate residual with disp increment.
@@ -857,13 +851,6 @@
   
   //mesh->debug(true); // DEBUGGING
   
-  spatialdata::geocoords::CSCart cs;
-  spatialdata::units::Nondimensional normalizer;
-  cs.setSpaceDim(mesh->dimension());
-  cs.initialize();
-  mesh->coordsys(&cs);
-  mesh->nondimensionalize(normalizer);
-  
   _quadrature->initialize(_data->basis, _data->numQuadPts, _data->numBasis,
 			  _data->basisDeriv,
 			  _data->numQuadPts, _data->numBasis, _data->cellDim,
@@ -903,8 +890,9 @@
     names[i][1] = '\0';
   } // for
   
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh->sieveMesh()->getIntSection(_data->label)->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh->sieveMesh()->getIntSection(_data->label)->size();
+  int firstFaultCell      = mesh->sieveMesh()->getIntSection(_data->label)->size();
   if (fault->useLagrangeConstraints()) {
     firstFaultCell += mesh->sieveMesh()->getIntSection(_data->label)->size();
   }
@@ -913,8 +901,15 @@
   fault->quadrature(_quadrature);
   
   fault->eqsrcs(const_cast<const char**>(names), nsrcs, sources, nsrcs);
-  fault->adjustTopology(mesh, &firstFaultVertex, &firstFaultCell, _flipFault);
+  fault->adjustTopology(mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
   
+  spatialdata::geocoords::CSCart cs;
+  spatialdata::units::Nondimensional normalizer;
+  cs.setSpaceDim(mesh->dimension());
+  cs.initialize();
+  mesh->coordsys(&cs);
+  mesh->nondimensionalize(normalizer);
+
   const double upDir[] = { 0.0, 0.0, 1.0 };
   const double normalDir[] = { 1.0, 0.0, 0.0 };
   

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestLiuCosSlipFn.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestLiuCosSlipFn.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestLiuCosSlipFn.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -393,8 +393,9 @@
   mesh->coordsys(&cs);
 
   // Create fault mesh
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh->sieveMesh()->getIntSection(faultLabel)->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh->sieveMesh()->getIntSection(faultLabel)->size();
+  int firstFaultCell      = mesh->sieveMesh()->getIntSection(faultLabel)->size();
   const bool useLagrangeConstraints = true;
   if (useLagrangeConstraints) {
     firstFaultCell += mesh->sieveMesh()->getIntSection(faultLabel)->size();
@@ -406,9 +407,8 @@
                                 *mesh, sieveMesh->getIntSection(faultLabel));
   CohesiveTopology::create(mesh, *faultMesh, faultBoundary, 
                            sieveMesh->getIntSection(faultLabel),
-                           faultId,
-                           firstFaultVertex, firstFaultCell,
-                           useLagrangeConstraints);
+                           faultId, firstFaultVertex, firstLagrangeVertex,
+			   firstFaultCell, useLagrangeConstraints);
   // Need to copy coordinates from mesh to fault mesh since we are not
   // using create() instead of createParallel().
   const ALE::Obj<SieveSubMesh>& faultSieveMesh = faultMesh->sieveMesh();
@@ -466,8 +466,9 @@
 
   // Create fault mesh
   topology::SubMesh faultMesh;
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
+  int firstFaultCell      = mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
   const bool useLagrangeConstraints = true;
   if (useLagrangeConstraints) {
     firstFaultCell += mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
@@ -479,9 +480,8 @@
                                 mesh, sieveMesh->getIntSection(data.faultLabel));
   CohesiveTopology::create(&mesh, faultMesh, faultBoundary, 
                            sieveMesh->getIntSection(data.faultLabel),
-                           data.faultId,
-                           firstFaultVertex, firstFaultCell,
-                           useLagrangeConstraints);
+                           data.faultId, firstFaultVertex, firstLagrangeVertex,
+			   firstFaultCell, useLagrangeConstraints);
   // Need to copy coordinates from mesh to fault mesh since we are not
   // using create() instead of createParallel().
   const ALE::Obj<SieveSubMesh>& faultSieveMesh = faultMesh.sieveMesh();

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestStepSlipFn.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestStepSlipFn.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/faults/TestStepSlipFn.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -309,8 +309,9 @@
   mesh->coordsys(&cs);
 
   // Create fault mesh
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh->sieveMesh()->getIntSection(faultLabel)->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh->sieveMesh()->getIntSection(faultLabel)->size();
+  int firstFaultCell      = mesh->sieveMesh()->getIntSection(faultLabel)->size();
   const bool useLagrangeConstraints = true;
   if (useLagrangeConstraints) {
     firstFaultCell += mesh->sieveMesh()->getIntSection(faultLabel)->size();
@@ -323,7 +324,7 @@
   CohesiveTopology::create(mesh, *faultMesh, faultBoundary, 
                            sieveMesh->getIntSection(faultLabel),
                            faultId,
-                           firstFaultVertex, firstFaultCell,
+                           firstFaultVertex, firstLagrangeVertex, firstFaultCell,
                            useLagrangeConstraints);
   // Need to copy coordinates from mesh to fault mesh since we are not
   // using create() instead of createParallel().
@@ -376,8 +377,9 @@
 
   // Create fault mesh
   topology::SubMesh faultMesh;
-  int firstFaultVertex = 0;
-  int firstFaultCell   = mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
+  int firstFaultCell      = mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
   const bool useLagrangeConstraints = true;
   if (useLagrangeConstraints) {
     firstFaultCell += mesh.sieveMesh()->getIntSection(data.faultLabel)->size();
@@ -390,7 +392,7 @@
   CohesiveTopology::create(&mesh, faultMesh, faultBoundary, 
                            sieveMesh->getIntSection(data.faultLabel),
                            data.faultId,
-                           firstFaultVertex, firstFaultCell,
+                           firstFaultVertex, firstLagrangeVertex, firstFaultCell,
                            useLagrangeConstraints);
   // Need to copy coordinates from mesh to fault mesh since we are not
   // using create() instead of createParallel().

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKBCMesh.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKBCMesh.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKBCMesh.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -208,14 +208,15 @@
 
   if (0 != _data->faultLabel) {
     faults::FaultCohesiveKin fault;
-    int firstFaultVertex = 0;
-    int firstFaultCell   = _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
+    int firstFaultVertex    = 0;
+    int firstLagrangeVertex = _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
+    int firstFaultCell      = _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
     if (fault.useLagrangeConstraints()) {
       firstFaultCell += _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
     }
     fault.label(_data->faultLabel);
     fault.id(_data->faultId);
-    fault.adjustTopology(_mesh, &firstFaultVertex, &firstFaultCell, _flipFault);
+    fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
   } // if
 
   delete _submesh; _submesh = new topology::SubMesh(*_mesh, _data->bcLabel);

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKFaultMesh.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKFaultMesh.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKFaultMesh.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -211,15 +211,17 @@
   iohandler.read(_mesh);
 
   faults::FaultCohesiveKin fault;
-  int firstFaultVertex = 0;
-  int firstFaultCell   = _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
+  int firstFaultVertex    = 0;
+  int firstLagrangeVertex = _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
+  int firstFaultCell      = _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
   const bool constraintCell = true;
   if (constraintCell) {
     firstFaultCell += _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
   }
   fault.label(_data->faultLabel);
   fault.id(_data->faultId);
-  fault.adjustTopology(_mesh, &firstFaultVertex, &firstFaultCell, _flipFault);
+  fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, 
+		       &firstFaultCell, _flipFault);
   std::map<Mesh::point_type, Mesh::point_type> cohesiveToFault;
   faults::CohesiveTopology::createFaultParallel(_faultMesh, &cohesiveToFault,
 						*_mesh, _data->faultId,

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKMesh.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKMesh.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKMesh.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -82,7 +82,7 @@
   const char* format = "%4.1f";
   writer.timeFormat(format);
   CPPUNIT_ASSERT_EQUAL(std::string(format), writer._timeFormat);
-} // testInterpolate
+} // testTimeFormat
 
 // ----------------------------------------------------------------------
 // Test timeConstant()
@@ -94,9 +94,21 @@
   const double value = 4.5;
   writer.timeConstant(value);
   CPPUNIT_ASSERT_EQUAL(value, writer._timeConstant);
-} // testInterpolate
+} // testTimeConstant
 
 // ----------------------------------------------------------------------
+// Test precision()
+void
+pylith::meshio::TestDataWriterVTKMesh::testPrecision(void)
+{ // testPrecision
+  DataWriterVTK<topology::Mesh, MeshField> writer;
+
+  const int value = 4;
+  writer.precision(value);
+  CPPUNIT_ASSERT_EQUAL(value, writer._precision);
+} // testPrecision
+
+// ----------------------------------------------------------------------
 // Test openTimeStep() and closeTimeStep()
 void
 pylith::meshio::TestDataWriterVTKMesh::testTimeStep(void)
@@ -269,14 +281,15 @@
 
   if (0 != _data->faultLabel) {
     faults::FaultCohesiveKin fault;
-    int firstFaultVertex = 0;
-    int firstFaultCell   = _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
+    int firstFaultVertex    = 0;
+    int firstLagrangeVertex = _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
+    int firstFaultCell      = _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
     if (fault.useLagrangeConstraints()) {
       firstFaultCell += _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
     }
     fault.label(_data->faultLabel);
     fault.id(_data->faultId);
-    fault.adjustTopology(_mesh, &firstFaultVertex, &firstFaultCell, _flipFault);
+    fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
   } // if
 } // _initialize
 

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKMesh.hh
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKMesh.hh	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKMesh.hh	2010-04-10 00:35:22 UTC (rev 16524)
@@ -43,6 +43,7 @@
   CPPUNIT_TEST( testFilename );
   CPPUNIT_TEST( testTimeFormat );
   CPPUNIT_TEST( testTimeConstant );
+  CPPUNIT_TEST( testPrecision );
   CPPUNIT_TEST( testVtkFilename );
 
   CPPUNIT_TEST_SUITE_END();
@@ -68,6 +69,9 @@
   /// Test timeConstant()
   void testTimeConstant(void);
 
+  /// Test precision()
+  void testPrecision(void);
+
   /// Test openTimeStep() and closeTimeStep()
   void testTimeStep(void);
 

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKSubMesh.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKSubMesh.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/TestDataWriterVTKSubMesh.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -209,14 +209,15 @@
 
   if (0 != _data->faultLabel) {
     faults::FaultCohesiveKin fault;
-    int firstFaultVertex = 0;
-    int firstFaultCell   = _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
+    int firstFaultVertex    = 0;
+    int firstLagrangeVertex = _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
+    int firstFaultCell      = _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
     if (fault.useLagrangeConstraints()) {
       firstFaultCell += _mesh->sieveMesh()->getIntSection(_data->faultLabel)->size();
     }
     fault.label(_data->faultLabel);
     fault.id(_data->faultId);
-    fault.adjustTopology(_mesh, &firstFaultVertex, &firstFaultCell, _flipFault);
+    fault.adjustTopology(_mesh, &firstFaultVertex, &firstLagrangeVertex, &firstFaultCell, _flipFault);
   } // if
 
   CPPUNIT_ASSERT(0 != _data->bcLabel);

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_bc_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_bc_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_bc_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,14 +3,14 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 8 double
--1 -1 1
--1 1 1
-0 -1 1
-0 1 1
-1 -1 1
-1 1 1
-0 -1 1
-0 1 1
+-1.000000e+00 -1.000000e+00 1.000000e+00
+-1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
+1.000000e+00 -1.000000e+00 1.000000e+00
+1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
 CELLS 2 10
 4  6 7 1 0
 4  4 5 3 2
@@ -19,13 +19,13 @@
 9
 CELL_DATA 2
 VECTORS traction double
-1.1 2.2 3.3
-4.4 5.5 6.6
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
+2.100000e+00
+3.200000e+00
 SCALARS other double 6
 LOOKUP_TABLE default
-1.2 2.3 3.4 4.5 5.6 6.7
-7.8 8.9 9 10.1 11.2 12.3
+1.200000e+00 2.300000e+00 3.400000e+00 4.500000e+00 5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00 9.000000e+00 1.010000e+01 1.120000e+01 1.230000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_bc_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_bc_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_bc_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,14 +3,14 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 8 double
--1 -1 1
--1 1 1
-0 -1 1
-0 1 1
-1 -1 1
-1 1 1
-0 -1 1
-0 1 1
+-1.000000e+00 -1.000000e+00 1.000000e+00
+-1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
+1.000000e+00 -1.000000e+00 1.000000e+00
+1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
 CELLS 2 10
 4  6 7 1 0
 4  4 5 3 2

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_bc_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_bc_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_bc_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,14 +3,14 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 8 double
--1 -1 1
--1 1 1
-0 -1 1
-0 1 1
-1 -1 1
-1 1 1
-0 -1 1
-0 1 1
+-1.000000e+00 -1.000000e+00 1.000000e+00
+-1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
+1.000000e+00 -1.000000e+00 1.000000e+00
+1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
 CELLS 2 10
 4  6 7 1 0
 4  4 5 3 2
@@ -19,31 +19,31 @@
 9
 POINT_DATA 8
 VECTORS displacements double
-1.1 2.2 3.3
-4.4 5.5 6.6
-7.7 8.8 9.9
-10.1 11.2 12.3
-1.2 2.3 3.4
-4.5 5.6 6.7
-7.8 8.9 9
-10.2 11.3 12.4
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
+7.700000e+00 8.800000e+00 9.900000e+00
+1.010000e+01 1.120000e+01 1.230000e+01
+1.200000e+00 2.300000e+00 3.400000e+00
+4.500000e+00 5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00 9.000000e+00
+1.020000e+01 1.130000e+01 1.240000e+01
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
-6.5
-7.6
-8.7
-9.8
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
+6.500000e+00
+7.600000e+00
+8.700000e+00
+9.800000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
-1.3 2.4
-3.5 4.6
-5.7 6.8
-7.9 8
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00
+1.300000e+00 2.400000e+00
+3.500000e+00 4.600000e+00
+5.700000e+00 6.800000e+00
+7.900000e+00 8.000000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,22 +3,22 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 16 double
--1 -1 -1
--1 1 -1
-0 -1 -1
-0 1 -1
-1 -1 -1
-1 1 -1
--1 -1 1
--1 1 1
-0 -1 1
-0 1 1
-1 -1 1
-1 1 1
-0 -1 -1
-0 1 -1
-0 -1 1
-0 1 1
+-1.000000e+00 -1.000000e+00 -1.000000e+00
+-1.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+1.000000e+00 -1.000000e+00 -1.000000e+00
+1.000000e+00 1.000000e+00 -1.000000e+00
+-1.000000e+00 -1.000000e+00 1.000000e+00
+-1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
+1.000000e+00 -1.000000e+00 1.000000e+00
+1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
 CELLS 2 18
 8  0 12 13 1 6 14 15 7
 8  2 4 5 3 8 10 11 9
@@ -27,13 +27,13 @@
 12
 CELL_DATA 2
 VECTORS traction double
-1.1 2.2 3.3
-4.4 5.5 6.6
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
+2.100000e+00
+3.200000e+00
 SCALARS other double 6
 LOOKUP_TABLE default
-1.2 2.3 3.4 4.5 5.6 6.7
-1.1 2.2 3.3 4.4 5.5 6.6
+1.200000e+00 2.300000e+00 3.400000e+00 4.500000e+00 5.600000e+00 6.700000e+00
+1.100000e+00 2.200000e+00 3.300000e+00 4.400000e+00 5.500000e+00 6.600000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_fault_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_fault_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_fault_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,20 +3,20 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
-0 -1 -1
-0 1 -1
-0 -1 1
-0 1 1
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
 CELLS 1 5
 4  2 3 1 0
 CELL_TYPES 1
 9
 CELL_DATA 1
 VECTORS traction double
-1.1 2.2 3.3
+1.100000e+00 2.200000e+00 3.300000e+00
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
+2.100000e+00
 SCALARS other double 6
 LOOKUP_TABLE default
-1.2 2.3 3.4 4.5 5.6 6.7
+1.200000e+00 2.300000e+00 3.400000e+00 4.500000e+00 5.600000e+00 6.700000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_fault_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_fault_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_fault_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
-0 -1 -1
-0 1 -1
-0 -1 1
-0 1 1
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
 CELLS 1 5
 4  2 3 1 0
 CELL_TYPES 1

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_fault_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_fault_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_fault_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,29 +3,29 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
-0 -1 -1
-0 1 -1
-0 -1 1
-0 1 1
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
 CELLS 1 5
 4  2 3 1 0
 CELL_TYPES 1
 9
 POINT_DATA 4
 VECTORS displacements double
-1.1 2.2 3.3
-4.4 5.5 6.6
-7.7 8.8 9.9
-10.1 11.2 12.3
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
+7.700000e+00 8.800000e+00 9.900000e+00
+1.010000e+01 1.120000e+01 1.230000e+01
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_mat_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_mat_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_mat_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,32 +3,32 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 16 double
--1 -1 -1
--1 1 -1
-0 -1 -1
-0 1 -1
-1 -1 -1
-1 1 -1
--1 -1 1
--1 1 1
-0 -1 1
-0 1 1
-1 -1 1
-1 1 1
-0 -1 -1
-0 1 -1
-0 -1 1
-0 1 1
+-1.000000e+00 -1.000000e+00 -1.000000e+00
+-1.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+1.000000e+00 -1.000000e+00 -1.000000e+00
+1.000000e+00 1.000000e+00 -1.000000e+00
+-1.000000e+00 -1.000000e+00 1.000000e+00
+-1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
+1.000000e+00 -1.000000e+00 1.000000e+00
+1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
 CELLS 1 9
 8  2 4 5 3 8 10 11 9
 CELL_TYPES 1
 12
 CELL_DATA 1
 VECTORS traction double
-1.1 2.2 3.3
+1.100000e+00 2.200000e+00 3.300000e+00
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
+2.100000e+00
 SCALARS other double 6
 LOOKUP_TABLE default
-1.2 2.3 3.4 4.5 5.6 6.7
+1.200000e+00 2.300000e+00 3.400000e+00 4.500000e+00 5.600000e+00 6.700000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_mat_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_mat_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_mat_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,22 +3,22 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 16 double
--1 -1 -1
--1 1 -1
-0 -1 -1
-0 1 -1
-1 -1 -1
-1 1 -1
--1 -1 1
--1 1 1
-0 -1 1
-0 1 1
-1 -1 1
-1 1 1
-0 -1 -1
-0 1 -1
-0 -1 1
-0 1 1
+-1.000000e+00 -1.000000e+00 -1.000000e+00
+-1.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+1.000000e+00 -1.000000e+00 -1.000000e+00
+1.000000e+00 1.000000e+00 -1.000000e+00
+-1.000000e+00 -1.000000e+00 1.000000e+00
+-1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
+1.000000e+00 -1.000000e+00 1.000000e+00
+1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
 CELLS 1 9
 8  2 4 5 3 8 10 11 9
 CELL_TYPES 1

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_mat_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_mat_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_mat_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,77 +3,77 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 16 double
--1 -1 -1
--1 1 -1
-0 -1 -1
-0 1 -1
-1 -1 -1
-1 1 -1
--1 -1 1
--1 1 1
-0 -1 1
-0 1 1
-1 -1 1
-1 1 1
-0 -1 -1
-0 1 -1
-0 -1 1
-0 1 1
+-1.000000e+00 -1.000000e+00 -1.000000e+00
+-1.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+1.000000e+00 -1.000000e+00 -1.000000e+00
+1.000000e+00 1.000000e+00 -1.000000e+00
+-1.000000e+00 -1.000000e+00 1.000000e+00
+-1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
+1.000000e+00 -1.000000e+00 1.000000e+00
+1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
 CELLS 1 9
 8  2 4 5 3 8 10 11 9
 CELL_TYPES 1
 12
 POINT_DATA 16
 VECTORS displacements double
-1.1 2.2 3.3
-4.4 5.5 6.6
-7.7 8.8 9.9
-10.1 11.2 12.3
-1.2 2.3 3.4
-4.5 5.6 6.7
-7.8 8.9 9
-10.2 11.3 12.4
-1.3 2.4 3.5
-4.6 5.7 6.8
-7.9 8 9.1
-10.2 11.3 12.4
-13.5 14.6 15.7
-16.8 17.9 18.1
-19.2 20.3 21.4
-22.5 23.6 24.7
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
+7.700000e+00 8.800000e+00 9.900000e+00
+1.010000e+01 1.120000e+01 1.230000e+01
+1.200000e+00 2.300000e+00 3.400000e+00
+4.500000e+00 5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00 9.000000e+00
+1.020000e+01 1.130000e+01 1.240000e+01
+1.300000e+00 2.400000e+00 3.500000e+00
+4.600000e+00 5.700000e+00 6.800000e+00
+7.900000e+00 8.000000e+00 9.100000e+00
+1.020000e+01 1.130000e+01 1.240000e+01
+1.350000e+01 1.460000e+01 1.570000e+01
+1.680000e+01 1.790000e+01 1.810000e+01
+1.920000e+01 2.030000e+01 2.140000e+01
+2.250000e+01 2.360000e+01 2.470000e+01
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
-6.5
-7.6
-8.7
-9.8
-10
-12.1
-11.1
-13.1
-14.1
-15.1
-16.1
-17.1
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
+6.500000e+00
+7.600000e+00
+8.700000e+00
+9.800000e+00
+1.000000e+01
+1.210000e+01
+1.110000e+01
+1.310000e+01
+1.410000e+01
+1.510000e+01
+1.610000e+01
+1.710000e+01
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
-1.3 2.4
-3.5 4.6
-5.7 6.8
-7.9 8
-1.3 2.4
-3.5 4.6
-5.7 6.8
-8 1.4
-2.5 3.6
-4.8 1.5
-2.6 3.7
-4.8 5.9
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00
+1.300000e+00 2.400000e+00
+3.500000e+00 4.600000e+00
+5.700000e+00 6.800000e+00
+7.900000e+00 8.000000e+00
+1.300000e+00 2.400000e+00
+3.500000e+00 4.600000e+00
+5.700000e+00 6.800000e+00
+8.000000e+00 1.400000e+00
+2.500000e+00 3.600000e+00
+4.800000e+00 1.500000e+00
+2.600000e+00 3.700000e+00
+4.800000e+00 5.900000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_surf_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_surf_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_surf_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 -1 1
--1 1 1
-0 -1 1
-0 1 1
-1 -1 1
-1 1 1
+-1.000000e+00 -1.000000e+00 1.000000e+00
+-1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
+1.000000e+00 -1.000000e+00 1.000000e+00
+1.000000e+00 1.000000e+00 1.000000e+00
 CELLS 2 10
 4  2 3 1 0
 4  4 5 3 2
@@ -17,13 +17,13 @@
 9
 CELL_DATA 2
 VECTORS traction double
-1.1 2.2 3.3
-4.4 5.5 6.6
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
+2.100000e+00
+3.200000e+00
 SCALARS other double 6
 LOOKUP_TABLE default
-1.2 2.3 3.4 4.5 5.6 6.7
-7.8 8.9 9 10.1 11.2 12.3
+1.200000e+00 2.300000e+00 3.400000e+00 4.500000e+00 5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00 9.000000e+00 1.010000e+01 1.120000e+01 1.230000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_surf_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_surf_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_surf_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 -1 1
--1 1 1
-0 -1 1
-0 1 1
-1 -1 1
-1 1 1
+-1.000000e+00 -1.000000e+00 1.000000e+00
+-1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
+1.000000e+00 -1.000000e+00 1.000000e+00
+1.000000e+00 1.000000e+00 1.000000e+00
 CELLS 2 10
 4  2 3 1 0
 4  4 5 3 2

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_surf_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_surf_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_surf_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 -1 1
--1 1 1
-0 -1 1
-0 1 1
-1 -1 1
-1 1 1
+-1.000000e+00 -1.000000e+00 1.000000e+00
+-1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
+1.000000e+00 -1.000000e+00 1.000000e+00
+1.000000e+00 1.000000e+00 1.000000e+00
 CELLS 2 10
 4  2 3 1 0
 4  4 5 3 2
@@ -17,25 +17,25 @@
 9
 POINT_DATA 6
 VECTORS displacements double
-7.8 8.9 9
-10.2 11.3 12.4
-1.3 2.4 3.5
-4.6 5.7 6.8
-7.9 8.1 9.2
-10.3 11.4 12.5
+7.800000e+00 8.900000e+00 9.000000e+00
+1.020000e+01 1.130000e+01 1.240000e+01
+1.300000e+00 2.400000e+00 3.500000e+00
+4.600000e+00 5.700000e+00 6.800000e+00
+7.900000e+00 8.100000e+00 9.200000e+00
+1.030000e+01 1.140000e+01 1.250000e+01
 SCALARS pressure double 1
 LOOKUP_TABLE default
-8.7
-9.8
-10.9
-11.8
-12.7
-13.6
+8.700000e+00
+9.800000e+00
+1.090000e+01
+1.180000e+01
+1.270000e+01
+1.360000e+01
 SCALARS other double 2
 LOOKUP_TABLE default
-5.7 6.8
-7.9 8
-8.1 8.2
-9.2 9.3
-10.4 10.5
-11.5 11.6
+5.700000e+00 6.800000e+00
+7.900000e+00 8.000000e+00
+8.100000e+00 8.200000e+00
+9.200000e+00 9.300000e+00
+1.040000e+01 1.050000e+01
+1.150000e+01 1.160000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,22 +3,22 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 16 double
--1 -1 -1
--1 1 -1
-0 -1 -1
-0 1 -1
-1 -1 -1
-1 1 -1
--1 -1 1
--1 1 1
-0 -1 1
-0 1 1
-1 -1 1
-1 1 1
-0 -1 -1
-0 1 -1
-0 -1 1
-0 1 1
+-1.000000e+00 -1.000000e+00 -1.000000e+00
+-1.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+1.000000e+00 -1.000000e+00 -1.000000e+00
+1.000000e+00 1.000000e+00 -1.000000e+00
+-1.000000e+00 -1.000000e+00 1.000000e+00
+-1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
+1.000000e+00 -1.000000e+00 1.000000e+00
+1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
 CELLS 2 18
 8  0 12 13 1 6 14 15 7
 8  2 4 5 3 8 10 11 9

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/hex8_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,22 +3,22 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 16 double
--1 -1 -1
--1 1 -1
-0 -1 -1
-0 1 -1
-1 -1 -1
-1 1 -1
--1 -1 1
--1 1 1
-0 -1 1
-0 1 1
-1 -1 1
-1 1 1
-0 -1 -1
-0 1 -1
-0 -1 1
-0 1 1
+-1.000000e+00 -1.000000e+00 -1.000000e+00
+-1.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+1.000000e+00 -1.000000e+00 -1.000000e+00
+1.000000e+00 1.000000e+00 -1.000000e+00
+-1.000000e+00 -1.000000e+00 1.000000e+00
+-1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
+1.000000e+00 -1.000000e+00 1.000000e+00
+1.000000e+00 1.000000e+00 1.000000e+00
+0.000000e+00 -1.000000e+00 -1.000000e+00
+0.000000e+00 1.000000e+00 -1.000000e+00
+0.000000e+00 -1.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 1.000000e+00
 CELLS 2 18
 8  0 12 13 1 6 14 15 7
 8  2 4 5 3 8 10 11 9
@@ -27,55 +27,55 @@
 12
 POINT_DATA 16
 VECTORS displacements double
-1.1 2.2 3.3
-4.4 5.5 6.6
-7.7 8.8 9.9
-10.1 11.2 12.3
-1.2 2.3 3.4
-4.5 5.6 6.7
-7.8 8.9 9
-10.2 11.3 12.4
-1.3 2.4 3.5
-4.6 5.7 6.8
-7.9 8 9.1
-10.2 11.3 12.4
-13.5 14.6 15.7
-16.8 17.9 18.1
-19.2 20.3 21.4
-22.5 23.6 24.7
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
+7.700000e+00 8.800000e+00 9.900000e+00
+1.010000e+01 1.120000e+01 1.230000e+01
+1.200000e+00 2.300000e+00 3.400000e+00
+4.500000e+00 5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00 9.000000e+00
+1.020000e+01 1.130000e+01 1.240000e+01
+1.300000e+00 2.400000e+00 3.500000e+00
+4.600000e+00 5.700000e+00 6.800000e+00
+7.900000e+00 8.000000e+00 9.100000e+00
+1.020000e+01 1.130000e+01 1.240000e+01
+1.350000e+01 1.460000e+01 1.570000e+01
+1.680000e+01 1.790000e+01 1.810000e+01
+1.920000e+01 2.030000e+01 2.140000e+01
+2.250000e+01 2.360000e+01 2.470000e+01
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
-6.5
-7.6
-8.7
-9.8
-10
-12.1
-11.1
-13.1
-14.1
-15.1
-16.1
-17.1
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
+6.500000e+00
+7.600000e+00
+8.700000e+00
+9.800000e+00
+1.000000e+01
+1.210000e+01
+1.110000e+01
+1.310000e+01
+1.410000e+01
+1.510000e+01
+1.610000e+01
+1.710000e+01
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
-1.3 2.4
-3.5 4.6
-5.7 6.8
-7.9 8
-1.3 2.4
-3.5 4.6
-5.7 6.8
-8 1.4
-2.5 3.6
-4.8 1.5
-2.6 3.7
-4.8 5.9
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00
+1.300000e+00 2.400000e+00
+3.500000e+00 4.600000e+00
+5.700000e+00 6.800000e+00
+7.900000e+00 8.000000e+00
+1.300000e+00 2.400000e+00
+3.500000e+00 4.600000e+00
+5.700000e+00 6.800000e+00
+8.000000e+00 1.400000e+00
+2.500000e+00 3.600000e+00
+4.800000e+00 1.500000e+00
+2.600000e+00 3.700000e+00
+4.800000e+00 5.900000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0.0 0.0
-0 0.0 0.0
-1 0.0 0.0
-0 0.0 0.0
+-1.000000e+00 0.0 0.0
+0.000000e+00 0.0 0.0
+1.000000e+00 0.0 0.0
+0.000000e+00 0.0 0.0
 CELLS 2 6
 2  0 3
 2  1 2
@@ -15,13 +15,13 @@
 3
 CELL_DATA 2
 VECTORS traction double
-1.1 0.0 0.0
-2.2 0.0 0.0
+1.100000e+00 0.0 0.0
+2.200000e+00 0.0 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-2.2
+2.100000e+00
+2.200000e+00
 SCALARS other double 1
 LOOKUP_TABLE default
-1.2
-2.3
+1.200000e+00
+2.300000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_mat_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_mat_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_mat_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,20 +3,20 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0.0 0.0
-0 0.0 0.0
-1 0.0 0.0
-0 0.0 0.0
+-1.000000e+00 0.0 0.0
+0.000000e+00 0.0 0.0
+1.000000e+00 0.0 0.0
+0.000000e+00 0.0 0.0
 CELLS 1 3
 2  0 3
 CELL_TYPES 1
 3
 CELL_DATA 1
 VECTORS traction double
-1.1 0.0 0.0
+1.100000e+00 0.0 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
+2.100000e+00
 SCALARS other double 1
 LOOKUP_TABLE default
-1.2
+1.200000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_mat_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_mat_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_mat_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0.0 0.0
-0 0.0 0.0
-1 0.0 0.0
-0 0.0 0.0
+-1.000000e+00 0.0 0.0
+0.000000e+00 0.0 0.0
+1.000000e+00 0.0 0.0
+0.000000e+00 0.0 0.0
 CELLS 1 3
 2  0 3
 CELL_TYPES 1

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_mat_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_mat_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_mat_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,29 +3,29 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0.0 0.0
-0 0.0 0.0
-1 0.0 0.0
-0 0.0 0.0
+-1.000000e+00 0.0 0.0
+0.000000e+00 0.0 0.0
+1.000000e+00 0.0 0.0
+0.000000e+00 0.0 0.0
 CELLS 1 3
 2  0 3
 CELL_TYPES 1
 3
 POINT_DATA 4
 VECTORS displacements double
-1.1 0.0 0.0
-2.2 0.0 0.0
-3.3 0.0 0.0
-4.4 0.0 0.0
+1.100000e+00 0.0 0.0
+2.200000e+00 0.0 0.0
+3.300000e+00 0.0 0.0
+4.400000e+00 0.0 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_surf_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_surf_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_surf_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,16 +3,16 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 1 double
--1 0.0 0.0
+-1.000000e+00 0.0 0.0
 CELLS 1 2
 1  0
 CELL_TYPES 1
 -1
 CELL_DATA 1
 VECTORS traction double
-1.1 2.2 0.0
+1.100000e+00 2.200000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
+2.100000e+00
 VECTORS other double
-1.2 2.3 3.4
+1.200000e+00 2.300000e+00 3.400000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_surf_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_surf_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_surf_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,7 +3,7 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 1 double
--1 0.0 0.0
+-1.000000e+00 0.0 0.0
 CELLS 1 2
 1  0
 CELL_TYPES 1

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_surf_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_surf_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_surf_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,17 +3,17 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 1 double
--1 0.0 0.0
+-1.000000e+00 0.0 0.0
 CELLS 1 2
 1  0
 CELL_TYPES 1
 -1
 POINT_DATA 1
 VECTORS displacements double
-1.1 2.2 0.0
+1.100000e+00 2.200000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
+2.100000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
+1.200000e+00 2.300000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0.0 0.0
-0 0.0 0.0
-1 0.0 0.0
-0 0.0 0.0
+-1.000000e+00 0.0 0.0
+0.000000e+00 0.0 0.0
+1.000000e+00 0.0 0.0
+0.000000e+00 0.0 0.0
 CELLS 2 6
 2  0 3
 2  1 2

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/line2_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0.0 0.0
-0 0.0 0.0
-1 0.0 0.0
-0 0.0 0.0
+-1.000000e+00 0.0 0.0
+0.000000e+00 0.0 0.0
+1.000000e+00 0.0 0.0
+0.000000e+00 0.0 0.0
 CELLS 2 6
 2  0 3
 2  1 2
@@ -15,19 +15,19 @@
 3
 POINT_DATA 4
 VECTORS displacements double
-1.1 0.0 0.0
-2.2 0.0 0.0
-3.3 0.0 0.0
-4.4 0.0 0.0
+1.100000e+00 0.0 0.0
+2.200000e+00 0.0 0.0
+3.300000e+00 0.0 0.0
+4.400000e+00 0.0 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_cell_filter_t12.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_cell_filter_t12.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_cell_filter_t12.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0 0.0
-0 -1 0.0
-0 1 0.0
-1 0 0.0
+-1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
 CELLS 2 8
 3  0 1 2
 3  1 3 2
@@ -16,5 +16,5 @@
 CELL_DATA 2
 SCALARS field data double 1
 LOOKUP_TABLE default
-1.125
-2.125
+1.125000e+00
+2.125000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_cell_t12.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_cell_t12.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_cell_t12.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0 0.0
-0 -1 0.0
-0 1 0.0
-1 0 0.0
+-1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
 CELLS 2 8
 3  0 1 2
 3  1 3 2
@@ -16,5 +16,5 @@
 CELL_DATA 2
 SCALARS field data double 2
 LOOKUP_TABLE default
-1.1 1.2
-2.1 2.2
+1.100000e+00 1.200000e+00
+2.100000e+00 2.200000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_t12.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_t12.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_t12.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0 0.0
-0 -1 0.0
-0 1 0.0
-1 0 0.0
+-1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
 CELLS 2 8
 3  0 1 2
 3  1 3 2

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_vertex_filter_t12.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_vertex_filter_t12.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_vertex_filter_t12.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0 0.0
-0 -1 0.0
-0 1 0.0
-1 0 0.0
+-1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
 CELLS 2 8
 3  0 1 2
 3  1 3 2
@@ -16,7 +16,7 @@
 POINT_DATA 4
 SCALARS field data double 1
 LOOKUP_TABLE default
-1.62788
-3.04138
-4.45533
-5.86941
+1.627882e+00
+3.041381e+00
+4.455334e+00
+5.869412e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_vertex_t12.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_vertex_t12.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/output_vertex_t12.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0 0.0
-0 -1 0.0
-0 1 0.0
-1 0 0.0
+-1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
 CELLS 2 8
 3  0 1 2
 3  1 3 2
@@ -15,7 +15,7 @@
 5
 POINT_DATA 4
 VECTORS field data double
-1.1 1.2 0.0
-2.1 2.2 0.0
-3.1 3.2 0.0
-4.1 4.2 0.0
+1.100000e+00 1.200000e+00 0.0
+2.100000e+00 2.200000e+00 0.0
+3.100000e+00 3.200000e+00 0.0
+4.100000e+00 4.200000e+00 0.0

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_bc_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_bc_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_bc_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,9 +3,9 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 3 double
--1 -1 0.0
-0 -1 0.0
-1 -1 0.0
+-1.000000e+00 -1.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+1.000000e+00 -1.000000e+00 0.0
 CELLS 2 6
 2  0 1
 2  1 2
@@ -14,12 +14,12 @@
 3
 CELL_DATA 2
 VECTORS traction double
-1.1 2.2 0.0
-3.3 4.4 0.0
+1.100000e+00 2.200000e+00 0.0
+3.300000e+00 4.400000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
+2.100000e+00
+3.200000e+00
 VECTORS other double
-1.2 2.3 3.4
-4.5 5.6 6.7
+1.200000e+00 2.300000e+00 3.400000e+00
+4.500000e+00 5.600000e+00 6.700000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_bc_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_bc_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_bc_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,9 +3,9 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 3 double
--1 -1 0.0
-0 -1 0.0
-1 -1 0.0
+-1.000000e+00 -1.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+1.000000e+00 -1.000000e+00 0.0
 CELLS 2 6
 2  0 1
 2  1 2

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_bc_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_bc_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_bc_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,9 +3,9 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 3 double
--1 -1 0.0
-0 -1 0.0
-1 -1 0.0
+-1.000000e+00 -1.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+1.000000e+00 -1.000000e+00 0.0
 CELLS 2 6
 2  0 1
 2  1 2
@@ -14,16 +14,16 @@
 3
 POINT_DATA 3
 VECTORS displacements double
-1.1 2.2 0.0
-3.3 4.4 0.0
-5.5 6.6 0.0
+1.100000e+00 2.200000e+00 0.0
+3.300000e+00 4.400000e+00 0.0
+5.500000e+00 6.600000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
+2.100000e+00
+3.200000e+00
+4.300000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 -1 0.0
--1 1 0.0
-0 -1 0.0
-0 1 0.0
-1 -1 0.0
-1 1 0.0
+-1.000000e+00 -1.000000e+00 0.0
+-1.000000e+00 1.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 -1.000000e+00 0.0
+1.000000e+00 1.000000e+00 0.0
 CELLS 2 10
 4  0 2 3 1
 4  2 4 5 3
@@ -17,12 +17,12 @@
 9
 CELL_DATA 2
 VECTORS traction double
-1.1 2.2 0.0
-3.3 4.4 0.0
+1.100000e+00 2.200000e+00 0.0
+3.300000e+00 4.400000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-2.2
+2.100000e+00
+2.200000e+00
 VECTORS other double
-1.2 2.3 3.4
-4.5 5.6 6.7
+1.200000e+00 2.300000e+00 3.400000e+00
+4.500000e+00 5.600000e+00 6.700000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_fault_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_fault_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_fault_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,17 +3,17 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 2 double
-0 -1 0.0
-0 1 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
 CELLS 1 3
 2  1 0
 CELL_TYPES 1
 3
 CELL_DATA 1
 VECTORS traction double
-1.1 2.2 0.0
+1.100000e+00 2.200000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
+2.100000e+00
 VECTORS other double
-1.2 2.3 3.4
+1.200000e+00 2.300000e+00 3.400000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_fault_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_fault_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_fault_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,8 +3,8 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 2 double
-0 -1 0.0
-0 1 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
 CELLS 1 3
 2  1 0
 CELL_TYPES 1

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_fault_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_fault_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_fault_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,21 +3,21 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 2 double
-0 -1 0.0
-0 1 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
 CELLS 1 3
 2  1 0
 CELL_TYPES 1
 3
 POINT_DATA 2
 VECTORS displacements double
-1.1 2.2 0.0
-3.3 4.4 0.0
+1.100000e+00 2.200000e+00 0.0
+3.300000e+00 4.400000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
+2.100000e+00
+3.200000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_mat_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_mat_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_mat_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,21 +3,21 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 -1 0.0
--1 1 0.0
-0 -1 0.0
-0 1 0.0
-1 -1 0.0
-1 1 0.0
+-1.000000e+00 -1.000000e+00 0.0
+-1.000000e+00 1.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 -1.000000e+00 0.0
+1.000000e+00 1.000000e+00 0.0
 CELLS 1 5
 4  2 4 5 3
 CELL_TYPES 1
 9
 CELL_DATA 1
 VECTORS traction double
-1.1 2.2 0.0
+1.100000e+00 2.200000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
+2.100000e+00
 VECTORS other double
-1.2 2.3 3.4
+1.200000e+00 2.300000e+00 3.400000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_mat_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_mat_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_mat_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 -1 0.0
--1 1 0.0
-0 -1 0.0
-0 1 0.0
-1 -1 0.0
-1 1 0.0
+-1.000000e+00 -1.000000e+00 0.0
+-1.000000e+00 1.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 -1.000000e+00 0.0
+1.000000e+00 1.000000e+00 0.0
 CELLS 1 5
 4  2 4 5 3
 CELL_TYPES 1

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_mat_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_mat_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_mat_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,37 +3,37 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 -1 0.0
--1 1 0.0
-0 -1 0.0
-0 1 0.0
-1 -1 0.0
-1 1 0.0
+-1.000000e+00 -1.000000e+00 0.0
+-1.000000e+00 1.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 -1.000000e+00 0.0
+1.000000e+00 1.000000e+00 0.0
 CELLS 1 5
 4  2 4 5 3
 CELL_TYPES 1
 9
 POINT_DATA 6
 VECTORS displacements double
-1.1 2.2 0.0
-3.3 4.4 0.0
-5.5 6.6 0.0
-7.7 8.8 0.0
-9.9 10.1 0.0
-11.2 12.3 0.0
+1.100000e+00 2.200000e+00 0.0
+3.300000e+00 4.400000e+00 0.0
+5.500000e+00 6.600000e+00 0.0
+7.700000e+00 8.800000e+00 0.0
+9.900000e+00 1.010000e+01 0.0
+1.120000e+01 1.230000e+01 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
-6.5
-7.6
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
+6.500000e+00
+7.600000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
-9.8 7.6
-6.5 5.4
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00
+9.800000e+00 7.600000e+00
+6.500000e+00 5.400000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_surf_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_surf_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_surf_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,9 +3,9 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 3 double
--1 -1 0.0
-0 -1 0.0
-1 -1 0.0
+-1.000000e+00 -1.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+1.000000e+00 -1.000000e+00 0.0
 CELLS 2 6
 2  0 1
 2  1 2
@@ -14,12 +14,12 @@
 3
 CELL_DATA 2
 VECTORS traction double
-1.1 2.2 0.0
-3.3 4.4 0.0
+1.100000e+00 2.200000e+00 0.0
+3.300000e+00 4.400000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
+2.100000e+00
+3.200000e+00
 VECTORS other double
-1.2 2.3 3.4
-4.5 5.6 6.7
+1.200000e+00 2.300000e+00 3.400000e+00
+4.500000e+00 5.600000e+00 6.700000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_surf_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_surf_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_surf_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,9 +3,9 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 3 double
--1 -1 0.0
-0 -1 0.0
-1 -1 0.0
+-1.000000e+00 -1.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+1.000000e+00 -1.000000e+00 0.0
 CELLS 2 6
 2  0 1
 2  1 2

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_surf_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_surf_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_surf_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,9 +3,9 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 3 double
--1 -1 0.0
-0 -1 0.0
-1 -1 0.0
+-1.000000e+00 -1.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+1.000000e+00 -1.000000e+00 0.0
 CELLS 2 6
 2  0 1
 2  1 2
@@ -14,16 +14,16 @@
 3
 POINT_DATA 3
 VECTORS displacements double
-1.1 2.2 0.0
-5.5 6.6 0.0
-9.9 10 0.0
+1.100000e+00 2.200000e+00 0.0
+5.500000e+00 6.600000e+00 0.0
+9.900000e+00 1.000000e+01 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-4.3
-6.5
+2.100000e+00
+4.300000e+00
+6.500000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-5.6 6.7
-9.1 10.2
+1.200000e+00 2.300000e+00
+5.600000e+00 6.700000e+00
+9.100000e+00 1.020000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 -1 0.0
--1 1 0.0
-0 -1 0.0
-0 1 0.0
-1 -1 0.0
-1 1 0.0
+-1.000000e+00 -1.000000e+00 0.0
+-1.000000e+00 1.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 -1.000000e+00 0.0
+1.000000e+00 1.000000e+00 0.0
 CELLS 2 10
 4  0 2 3 1
 4  2 4 5 3

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/quad4_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 -1 0.0
--1 1 0.0
-0 -1 0.0
-0 1 0.0
-1 -1 0.0
-1 1 0.0
+-1.000000e+00 -1.000000e+00 0.0
+-1.000000e+00 1.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 -1.000000e+00 0.0
+1.000000e+00 1.000000e+00 0.0
 CELLS 2 10
 4  0 2 3 1
 4  2 4 5 3
@@ -17,25 +17,25 @@
 9
 POINT_DATA 6
 VECTORS displacements double
-1.1 2.2 0.0
-3.3 4.4 0.0
-5.5 6.6 0.0
-7.7 8.8 0.0
-9.9 10.1 0.0
-11.2 12.3 0.0
+1.100000e+00 2.200000e+00 0.0
+3.300000e+00 4.400000e+00 0.0
+5.500000e+00 6.600000e+00 0.0
+7.700000e+00 8.800000e+00 0.0
+9.900000e+00 1.010000e+01 0.0
+1.120000e+01 1.230000e+01 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
-6.5
-7.6
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
+6.500000e+00
+7.600000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
-9.8 7.6
-6.5 5.4
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00
+9.800000e+00 7.600000e+00
+6.500000e+00 5.400000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_bc_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_bc_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_bc_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 0 0
-0 -1 0
-0 1 0
-1 0 0
-0 -1 0
-0 1 0
+-1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
+1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
 CELLS 2 8
 3  0 2 1
 3  4 5 3
@@ -17,13 +17,13 @@
 5
 CELL_DATA 2
 VECTORS traction double
-1.1 2.2 3.3
-4.4 5.5 6.6
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
+2.100000e+00
+3.200000e+00
 SCALARS other double 6
 LOOKUP_TABLE default
-1.2 2.3 3.4 4.5 5.6 6.7
-7.8 8.9 9 10.1 11.2 12.3
+1.200000e+00 2.300000e+00 3.400000e+00 4.500000e+00 5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00 9.000000e+00 1.010000e+01 1.120000e+01 1.230000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_bc_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_bc_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_bc_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 0 0
-0 -1 0
-0 1 0
-1 0 0
-0 -1 0
-0 1 0
+-1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
+1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
 CELLS 2 8
 3  0 2 1
 3  4 5 3

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_bc_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_bc_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_bc_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 0 0
-0 -1 0
-0 1 0
-1 0 0
-0 -1 0
-0 1 0
+-1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
+1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
 CELLS 2 8
 3  0 2 1
 3  4 5 3
@@ -17,25 +17,25 @@
 5
 POINT_DATA 6
 VECTORS displacements double
-1.1 2.2 3.3
-4.4 5.5 6.6
-7.7 8.8 9.9
-10 11.1 12.2
-13.3 14.4 15.5
-16.6 17.7 18.8
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
+7.700000e+00 8.800000e+00 9.900000e+00
+1.000000e+01 1.110000e+01 1.220000e+01
+1.330000e+01 1.440000e+01 1.550000e+01
+1.660000e+01 1.770000e+01 1.880000e+01
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
-6.5
-7.6
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
+6.500000e+00
+7.600000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
-9 10.1
-11.2 12.3
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00
+9.000000e+00 1.010000e+01
+1.120000e+01 1.230000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,14 +3,14 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 8 double
--1 0 0
-0 -1 0
-0 0 1
-0 1 0
-1 0 0
-0 -1 0
-0 0 1
-0 1 0
+-1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
+1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
 CELLS 2 10
 4  1 2 3 0
 4  5 7 6 4
@@ -19,13 +19,13 @@
 10
 CELL_DATA 2
 VECTORS traction double
-1.1 2.2 3.3
-4.4 5.5 6.6
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
+2.100000e+00
+3.200000e+00
 SCALARS other double 6
 LOOKUP_TABLE default
-1.2 2.3 3.4 4.5 5.6 6.7
-7.8 8.9 9 10.1 11.2 12.3
+1.200000e+00 2.300000e+00 3.400000e+00 4.500000e+00 5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00 9.000000e+00 1.010000e+01 1.120000e+01 1.230000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_fault_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_fault_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_fault_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,19 +3,19 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 3 double
-0 -1 0
-0 0 1
-0 1 0
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
 CELLS 1 4
 3  1 0 2
 CELL_TYPES 1
 5
 CELL_DATA 1
 VECTORS traction double
-1.1 2.2 3.3
+1.100000e+00 2.200000e+00 3.300000e+00
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
+2.100000e+00
 SCALARS other double 6
 LOOKUP_TABLE default
-1.2 2.3 3.4 4.5 5.6 6.7
+1.200000e+00 2.300000e+00 3.400000e+00 4.500000e+00 5.600000e+00 6.700000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_fault_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_fault_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_fault_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,9 +3,9 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 3 double
-0 -1 0
-0 0 1
-0 1 0
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
 CELLS 1 4
 3  1 0 2
 CELL_TYPES 1

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_fault_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_fault_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_fault_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,25 +3,25 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 3 double
-0 -1 0
-0 0 1
-0 1 0
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
 CELLS 1 4
 3  1 0 2
 CELL_TYPES 1
 5
 POINT_DATA 3
 VECTORS displacements double
-1.1 2.2 3.3
-4.4 5.5 6.6
-7.7 8.8 9.9
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
+7.700000e+00 8.800000e+00 9.900000e+00
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
+2.100000e+00
+3.200000e+00
+4.300000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_mat_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_mat_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_mat_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,14 +3,14 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 8 double
--1 0 0
-0 -1 0
-0 0 1
-0 1 0
-1 0 0
-0 -1 0
-0 0 1
-0 1 0
+-1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
+1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
 CELLS 2 10
 4  1 2 3 0
 4  5 7 6 4
@@ -19,13 +19,13 @@
 10
 CELL_DATA 2
 VECTORS traction double
-1.1 2.2 3.3
-4.4 5.5 6.6
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
+2.100000e+00
+3.200000e+00
 SCALARS other double 6
 LOOKUP_TABLE default
-1.2 2.3 3.4 4.5 5.6 6.7
-7.8 8.9 9 10.1 11.2 12.3
+1.200000e+00 2.300000e+00 3.400000e+00 4.500000e+00 5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00 9.000000e+00 1.010000e+01 1.120000e+01 1.230000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_mat_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_mat_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_mat_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,14 +3,14 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 8 double
--1 0 0
-0 -1 0
-0 0 1
-0 1 0
-1 0 0
-0 -1 0
-0 0 1
-0 1 0
+-1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
+1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
 CELLS 2 10
 4  1 2 3 0
 4  5 7 6 4

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_mat_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_mat_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_mat_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,14 +3,14 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 8 double
--1 0 0
-0 -1 0
-0 0 1
-0 1 0
-1 0 0
-0 -1 0
-0 0 1
-0 1 0
+-1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
+1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
 CELLS 2 10
 4  1 2 3 0
 4  5 7 6 4
@@ -19,31 +19,31 @@
 10
 POINT_DATA 8
 VECTORS displacements double
-1.1 2.2 3.3
-4.4 5.5 6.6
-7.7 8.8 9.9
-10 11.1 12.2
-13.3 14.4 15.5
-16.6 17.7 18.8
-19.9 20 21.1
-22.2 23.3 24.4
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
+7.700000e+00 8.800000e+00 9.900000e+00
+1.000000e+01 1.110000e+01 1.220000e+01
+1.330000e+01 1.440000e+01 1.550000e+01
+1.660000e+01 1.770000e+01 1.880000e+01
+1.990000e+01 2.000000e+01 2.110000e+01
+2.220000e+01 2.330000e+01 2.440000e+01
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
-6.5
-7.6
-8.7
-9.8
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
+6.500000e+00
+7.600000e+00
+8.700000e+00
+9.800000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
-9 10.1
-11.2 12.3
-13.4 14.5
-15.6 16.7
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00
+9.000000e+00 1.010000e+01
+1.120000e+01 1.230000e+01
+1.340000e+01 1.450000e+01
+1.560000e+01 1.670000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_surf_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_surf_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_surf_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0 0
-0 -1 0
-0 1 0
-1 0 0
+-1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
+1.000000e+00 0.000000e+00 0.000000e+00
 CELLS 2 8
 3  0 2 1
 3  1 2 3
@@ -15,13 +15,13 @@
 5
 CELL_DATA 2
 VECTORS traction double
-1.1 2.2 3.3
-4.4 5.5 6.6
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
+2.100000e+00
+3.200000e+00
 SCALARS other double 6
 LOOKUP_TABLE default
-1.2 2.3 3.4 4.5 5.6 6.7
-7.8 8.9 9 10.1 11.2 12.3
+1.200000e+00 2.300000e+00 3.400000e+00 4.500000e+00 5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00 9.000000e+00 1.010000e+01 1.120000e+01 1.230000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_surf_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_surf_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_surf_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0 0
-0 -1 0
-0 1 0
-1 0 0
+-1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
+1.000000e+00 0.000000e+00 0.000000e+00
 CELLS 2 8
 3  0 2 1
 3  1 2 3

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_surf_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_surf_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_surf_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,10 +3,10 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
--1 0 0
-0 -1 0
-0 1 0
-1 0 0
+-1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
+1.000000e+00 0.000000e+00 0.000000e+00
 CELLS 2 8
 3  0 2 1
 3  1 2 3
@@ -15,19 +15,19 @@
 5
 POINT_DATA 4
 VECTORS displacements double
-1.1 2.2 3.3
-4.4 5.5 6.6
-10 11.1 12.2
-13.3 14.4 15.5
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
+1.000000e+01 1.110000e+01 1.220000e+01
+1.330000e+01 1.440000e+01 1.550000e+01
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-5.4
-6.5
+2.100000e+00
+3.200000e+00
+5.400000e+00
+6.500000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-7.8 8.9
-9 10.1
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+7.800000e+00 8.900000e+00
+9.000000e+00 1.010000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,14 +3,14 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 8 double
--1 0 0
-0 -1 0
-0 0 1
-0 1 0
-1 0 0
-0 -1 0
-0 0 1
-0 1 0
+-1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
+1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
 CELLS 2 10
 4  1 2 3 0
 4  5 7 6 4

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tet4_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,14 +3,14 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 8 double
--1 0 0
-0 -1 0
-0 0 1
-0 1 0
-1 0 0
-0 -1 0
-0 0 1
-0 1 0
+-1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
+1.000000e+00 0.000000e+00 0.000000e+00
+0.000000e+00 -1.000000e+00 0.000000e+00
+0.000000e+00 0.000000e+00 1.000000e+00
+0.000000e+00 1.000000e+00 0.000000e+00
 CELLS 2 10
 4  1 2 3 0
 4  5 7 6 4
@@ -19,31 +19,31 @@
 10
 POINT_DATA 8
 VECTORS displacements double
-1.1 2.2 3.3
-4.4 5.5 6.6
-7.7 8.8 9.9
-10 11.1 12.2
-13.3 14.4 15.5
-16.6 17.7 18.8
-19.9 20 21.1
-22.2 23.3 24.4
+1.100000e+00 2.200000e+00 3.300000e+00
+4.400000e+00 5.500000e+00 6.600000e+00
+7.700000e+00 8.800000e+00 9.900000e+00
+1.000000e+01 1.110000e+01 1.220000e+01
+1.330000e+01 1.440000e+01 1.550000e+01
+1.660000e+01 1.770000e+01 1.880000e+01
+1.990000e+01 2.000000e+01 2.110000e+01
+2.220000e+01 2.330000e+01 2.440000e+01
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
-6.5
-7.6
-8.7
-9.8
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
+6.500000e+00
+7.600000e+00
+8.700000e+00
+9.800000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
-9 10.1
-11.2 12.3
-13.4 14.5
-15.6 16.7
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00
+9.000000e+00 1.010000e+01
+1.120000e+01 1.230000e+01
+1.340000e+01 1.450000e+01
+1.560000e+01 1.670000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_bc_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_bc_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_bc_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,17 +3,17 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 2 double
-0 -1 0.0
-1 0 0.0
+0.000000e+00 -1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
 CELLS 1 3
 2  0 1
 CELL_TYPES 1
 3
 CELL_DATA 1
 VECTORS traction double
-1.1 2.2 0.0
+1.100000e+00 2.200000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
+2.100000e+00
 VECTORS other double
-1.2 2.3 3.4
+1.200000e+00 2.300000e+00 3.400000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_bc_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_bc_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_bc_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,8 +3,8 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 2 double
-0 -1 0.0
-1 0 0.0
+0.000000e+00 -1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
 CELLS 1 3
 2  0 1
 CELL_TYPES 1

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_bc_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_bc_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_bc_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,21 +3,21 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 2 double
-0 -1 0.0
-1 0 0.0
+0.000000e+00 -1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
 CELLS 1 3
 2  0 1
 CELL_TYPES 1
 3
 POINT_DATA 2
 VECTORS displacements double
-1.1 2.2 0.0
-3.3 4.4 0.0
+1.100000e+00 2.200000e+00 0.0
+3.300000e+00 4.400000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
+2.100000e+00
+3.200000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 0 0.0
-0 -1 0.0
-0 1 0.0
-1 0 0.0
-0 -1 0.0
-0 1 0.0
+-1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
 CELLS 2 8
 3  0 4 5
 3  1 3 2
@@ -17,12 +17,12 @@
 5
 CELL_DATA 2
 VECTORS traction double
-1.1 2.2 0.0
-3.3 4.4 0.0
+1.100000e+00 2.200000e+00 0.0
+3.300000e+00 4.400000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-2.2
+2.100000e+00
+2.200000e+00
 VECTORS other double
-1.2 2.3 3.4
-4.5 5.6 6.7
+1.200000e+00 2.300000e+00 3.400000e+00
+4.500000e+00 5.600000e+00 6.700000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_fault_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_fault_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_fault_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,17 +3,17 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 2 double
-0 -1 0.0
-0 1 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
 CELLS 1 3
 2  1 0
 CELL_TYPES 1
 3
 CELL_DATA 1
 VECTORS traction double
-1.1 2.2 0.0
+1.100000e+00 2.200000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
+2.100000e+00
 VECTORS other double
-1.2 2.3 3.4
+1.200000e+00 2.300000e+00 3.400000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_fault_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_fault_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_fault_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,8 +3,8 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 2 double
-0 -1 0.0
-0 1 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
 CELLS 1 3
 2  1 0
 CELL_TYPES 1

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_fault_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_fault_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_fault_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,21 +3,21 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 2 double
-0 -1 0.0
-0 1 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
 CELLS 1 3
 2  1 0
 CELL_TYPES 1
 3
 POINT_DATA 2
 VECTORS displacements double
-1.1 2.2 0.0
-3.3 4.4 0.0
+1.100000e+00 2.200000e+00 0.0
+3.300000e+00 4.400000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
+2.100000e+00
+3.200000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_mat_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_mat_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_mat_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,21 +3,21 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 0 0.0
-0 -1 0.0
-0 1 0.0
-1 0 0.0
-0 -1 0.0
-0 1 0.0
+-1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
 CELLS 1 4
 3  1 3 2
 CELL_TYPES 1
 5
 CELL_DATA 1
 VECTORS traction double
-1.1 2.2 0.0
+1.100000e+00 2.200000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
+2.100000e+00
 VECTORS other double
-1.2 2.3 3.4
+1.200000e+00 2.300000e+00 3.400000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_mat_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_mat_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_mat_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 0 0.0
-0 -1 0.0
-0 1 0.0
-1 0 0.0
-0 -1 0.0
-0 1 0.0
+-1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
 CELLS 1 4
 3  1 3 2
 CELL_TYPES 1

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_mat_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_mat_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_mat_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,37 +3,37 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 0 0.0
-0 -1 0.0
-0 1 0.0
-1 0 0.0
-0 -1 0.0
-0 1 0.0
+-1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
 CELLS 1 4
 3  1 3 2
 CELL_TYPES 1
 5
 POINT_DATA 6
 VECTORS displacements double
-1.1 2.2 0.0
-3.3 4.4 0.0
-5.5 6.6 0.0
-7.7 8.8 0.0
-9.9 10 0.0
-11.1 12.2 0.0
+1.100000e+00 2.200000e+00 0.0
+3.300000e+00 4.400000e+00 0.0
+5.500000e+00 6.600000e+00 0.0
+7.700000e+00 8.800000e+00 0.0
+9.900000e+00 1.000000e+01 0.0
+1.110000e+01 1.220000e+01 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
-6.5
-7.6
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
+6.500000e+00
+7.600000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
-9 10.1
-11.2 12.3
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00
+9.000000e+00 1.010000e+01
+1.120000e+01 1.230000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_surf_cell_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_surf_cell_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_surf_cell_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,17 +3,17 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 2 double
-0 -1 0.0
-1 0 0.0
+0.000000e+00 -1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
 CELLS 1 3
 2  0 1
 CELL_TYPES 1
 3
 CELL_DATA 1
 VECTORS traction double
-1.1 2.2 0.0
+1.100000e+00 2.200000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
+2.100000e+00
 VECTORS other double
-1.2 2.3 3.4
+1.200000e+00 2.300000e+00 3.400000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_surf_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_surf_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_surf_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,8 +3,8 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 2 double
-0 -1 0.0
-1 0 0.0
+0.000000e+00 -1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
 CELLS 1 3
 2  0 1
 CELL_TYPES 1

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_surf_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_surf_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_surf_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,21 +3,21 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 2 double
-0 -1 0.0
-1 0 0.0
+0.000000e+00 -1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
 CELLS 1 3
 2  0 1
 CELL_TYPES 1
 3
 POINT_DATA 2
 VECTORS displacements double
-3.3 4.4 0.0
-7.7 8.8 0.0
+3.300000e+00 4.400000e+00 0.0
+7.700000e+00 8.800000e+00 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-3.2
-5.4
+3.200000e+00
+5.400000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-3.4 4.5
-7.8 8.9
+3.400000e+00 4.500000e+00
+7.800000e+00 8.900000e+00

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 0 0.0
-0 -1 0.0
-0 1 0.0
-1 0 0.0
-0 -1 0.0
-0 1 0.0
+-1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
 CELLS 2 8
 3  0 4 5
 3  1 3 2

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_vertex_t10.vtk
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_vertex_t10.vtk	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/meshio/data/tri3_vertex_t10.vtk	2010-04-10 00:35:22 UTC (rev 16524)
@@ -3,12 +3,12 @@
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 6 double
--1 0 0.0
-0 -1 0.0
-0 1 0.0
-1 0 0.0
-0 -1 0.0
-0 1 0.0
+-1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
+1.000000e+00 0.000000e+00 0.0
+0.000000e+00 -1.000000e+00 0.0
+0.000000e+00 1.000000e+00 0.0
 CELLS 2 8
 3  0 4 5
 3  1 3 2
@@ -17,25 +17,25 @@
 5
 POINT_DATA 6
 VECTORS displacements double
-1.1 2.2 0.0
-3.3 4.4 0.0
-5.5 6.6 0.0
-7.7 8.8 0.0
-9.9 10 0.0
-11.1 12.2 0.0
+1.100000e+00 2.200000e+00 0.0
+3.300000e+00 4.400000e+00 0.0
+5.500000e+00 6.600000e+00 0.0
+7.700000e+00 8.800000e+00 0.0
+9.900000e+00 1.000000e+01 0.0
+1.110000e+01 1.220000e+01 0.0
 SCALARS pressure double 1
 LOOKUP_TABLE default
-2.1
-3.2
-4.3
-5.4
-6.5
-7.6
+2.100000e+00
+3.200000e+00
+4.300000e+00
+5.400000e+00
+6.500000e+00
+7.600000e+00
 SCALARS other double 2
 LOOKUP_TABLE default
-1.2 2.3
-3.4 4.5
-5.6 6.7
-7.8 8.9
-9 10.1
-11.2 12.3
+1.200000e+00 2.300000e+00
+3.400000e+00 4.500000e+00
+5.600000e+00 6.700000e+00
+7.800000e+00 8.900000e+00
+9.000000e+00 1.010000e+01
+1.120000e+01 1.230000e+01

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/topology/TestMeshOps.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/topology/TestMeshOps.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/topology/TestMeshOps.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -18,6 +18,7 @@
 
 #include "pylith/topology/Mesh.hh" // USES Mesh
 #include "pylith/meshio/MeshIOAscii.hh" // USES MeshIOAscii
+#include <stdexcept> // USES std::runtime_error
 
 // ----------------------------------------------------------------------
 CPPUNIT_TEST_SUITE_REGISTRATION( pylith::topology::TestMeshOps );

Modified: short/3D/PyLith/tags/v1.4.3/unittests/libtests/utils/TestEventLogger.cc
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/libtests/utils/TestEventLogger.cc	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/libtests/utils/TestEventLogger.cc	2010-04-10 00:35:22 UTC (rev 16524)
@@ -146,7 +146,7 @@
   logger.className("my class");
   logger.initialize();
 
-  const char* stages[] = { "stage A", "stage B", "stage C" };
+  const char* stages[] = { "stage A1", "stage B1", "stage C1" };
   const int numStages = 3;
   int ids[numStages];
 
@@ -171,7 +171,7 @@
   logger.className("my class");
   logger.initialize();
 
-  const char* stages[] = { "stage A", "stage B", "stage C" };
+  const char* stages[] = { "stage A2", "stage B2", "stage C2" };
   const int numStages = 3;
 
   for (int i=0; i < numStages; ++i)
@@ -189,8 +189,8 @@
        ++s_iter, ++i)
     CPPUNIT_ASSERT_EQUAL(s_iter->second, ids[i]);
   
-  const int idNew = logger.stageId("stage D");
-  CPPUNIT_ASSERT_EQUAL(idNew, logger.stageId("stage D"));
+  const int idNew = logger.stageId("stage D2");
+  CPPUNIT_ASSERT_EQUAL(idNew, logger.stageId("stage D2"));
 } // testStageId
 
 // ----------------------------------------------------------------------
@@ -202,7 +202,7 @@
   logger.className("my class");
   logger.initialize();
 
-  const char* stages[] = { "stage A", "stage B", "stage C" };
+  const char* stages[] = { "stage A3", "stage B3", "stage C3" };
   const int numStages = 3;
   int ids[numStages];
 

Modified: short/3D/PyLith/tags/v1.4.3/unittests/pytests/faults/TestFaultCohesiveKin.py
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/pytests/faults/TestFaultCohesiveKin.py	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/pytests/faults/TestFaultCohesiveKin.py	2010-04-10 00:35:22 UTC (rev 16524)
@@ -108,9 +108,12 @@
     fault.inventory.faultLabel = "fault"
     fault._configure()
 
+    nvertices = fault.numVertices(mesh)
     firstFaultVertex = 0
-    firstFaultCell = 0
-    fault.adjustTopology(mesh, firstFaultVertex, firstFaultCell)
+    firstLagrangeVertex = nvertices
+    firstFaultCell      = 2*nvertices
+    fault.adjustTopology(mesh, firstFaultVertex, firstLagrangeVertex,
+                         firstFaultCell)
 
     # We should really add something here to check to make sure things
     # actually initialized correctly    
@@ -336,9 +339,12 @@
     eqsrc.inventory.slipfn = slipfn
     eqsrc._configure()
 
+    nvertices = fault.numVertices(mesh)
     firstFaultVertex = 0
-    firstFaultCell = 0
-    fault.adjustTopology(mesh, firstFaultVertex, firstFaultCell)
+    firstLagrangeVertex = nvertices
+    firstFaultCell      = 2*nvertices
+    fault.adjustTopology(mesh, firstFaultVertex, firstLagrangeVertex,
+                         firstFaultCell)
     fault.preinitialize(mesh)
     fault.timeStep(dt)
     fault.verifyConfiguration()

Modified: short/3D/PyLith/tags/v1.4.3/unittests/pytests/utils/TestEventLogger.py
===================================================================
--- short/3D/PyLith/tags/v1.4.3/unittests/pytests/utils/TestEventLogger.py	2010-04-09 21:12:26 UTC (rev 16523)
+++ short/3D/PyLith/tags/v1.4.3/unittests/pytests/utils/TestEventLogger.py	2010-04-10 00:35:22 UTC (rev 16524)
@@ -108,7 +108,7 @@
     logger.className("logging A")
     logger.initialize()
 
-    stages = ["stage 1" , "stage 2" , "stage 3"]
+    stages = ["stage 1a" , "stage 2a" , "stage 3a"]
     id = {}
     for stage in stages:
       id[stage] = logger.registerStage(stage)
@@ -125,18 +125,18 @@
     logger = EventLogger()
     logger.className("logging A")
     logger.initialize()
-    stages = ["stage 1" , "stage 2" , "stage 3"]
+    stages = ["stage 1b" , "stage 2b" , "stage 3b"]
     for stage in stages:
       logger.registerStage(stage)
 
-    logger.stagePush("stage 2")
+    logger.stagePush("stage 2b")
     logger.stagePop()
 
-    logger.stagePush("stage 1")
+    logger.stagePush("stage 1b")
     logger.stagePop()
 
-    logger.stagePush("stage 3")
-    logger.stagePush("stage 1")
+    logger.stagePush("stage 3b")
+    logger.stagePush("stage 1b")
     logger.stagePop()
     logger.stagePop()
     return



More information about the CIG-COMMITS mailing list