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

[brad at geodynamics.org]

Author: brad
Date: 2010-09-13 11:47:41 -0700 (Mon, 13 Sep 2010)
New Revision: 17188

Modified:
   short/3D/PyLith/trunk/libsrc/topology/RefineUniform.cc
Log:
Fixed small bug in refinement.

Modified: short/3D/PyLith/trunk/libsrc/topology/RefineUniform.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/topology/RefineUniform.cc	2010-09-13 14:07:10 UTC (rev 17187)
+++ short/3D/PyLith/trunk/libsrc/topology/RefineUniform.cc	2010-09-13 18:47:41 UTC (rev 17188)
@@ -108,19 +108,18 @@
   assert(!sieveMesh.isNull());
   const ALE::Obj<SieveMesh>& newSieveMesh = newMesh->sieveMesh();
   assert(!newSieveMesh.isNull());
+  ALE::Obj<SieveMesh::sieve_type> newSieve =
+    new SieveMesh::sieve_type(mesh.comm(), mesh.debug());
+  newSieveMesh->setSieve(newSieve);
 
   std::map<edge_type, point_type> edge2vertex;
    
 #if 0 // ORIGINAL VERSION
-  ALE::Obj<SieveMesh::sieve_type> newSieve =
-    new SieveMesh::sieve_type(mesh.comm(), mesh.debug());
-
-  newSieveMesh->setSieve(newSieve);
   ALE::MeshBuilder<Mesh>::refineTetrahedra(*mesh.sieveMesh(), * newSieveMesh,
 					   edge2vertex);
 #else
   // Is arg to CellRefiner constructor the new mesh or the old mesh?
-  ALE::MeshBuilder<SieveMesh>::CellRefiner<SieveMesh,edge_type> refiner(*newSieveMesh);
+  ALE::MeshBuilder<SieveMesh>::CellRefiner<SieveMesh,edge_type> refiner(*sieveMesh);
 
   ALE::MeshBuilder<SieveMesh>::refineGeneral< SieveMesh,
     ALE::MeshBuilder<SieveMesh>::CellRefiner<SieveMesh,edge_type> >(*sieveMesh, *newSieveMesh, refiner);