[cig-commits] r16463 - short/3D/PyLith/trunk/libsrc/topology

[brad at geodynamics.org]

Author: brad
Date: 2010-03-30 14:39:56 -0700 (Tue, 30 Mar 2010)
New Revision: 16463

Modified:
   short/3D/PyLith/trunk/libsrc/topology/Mesh.cc
Log:
Small cleanup of nondimensionalization of coordinates.

Modified: short/3D/PyLith/trunk/libsrc/topology/Mesh.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/topology/Mesh.cc	2010-03-30 21:39:19 UTC (rev 16462)
+++ short/3D/PyLith/trunk/libsrc/topology/Mesh.cc	2010-03-30 21:39:56 UTC (rev 16463)
@@ -120,15 +120,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