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

[brad at geodynamics.org]

Author: brad
Date: 2010-08-23 16:49:58 -0700 (Mon, 23 Aug 2010)
New Revision: 17111

Modified:
   short/3D/PyLith/trunk/libsrc/topology/RefineUniform.cc
Log:
Started filling in code to call revised Sieve uniform global refinement stuff.

Modified: short/3D/PyLith/trunk/libsrc/topology/RefineUniform.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/topology/RefineUniform.cc	2010-08-23 23:49:12 UTC (rev 17110)
+++ short/3D/PyLith/trunk/libsrc/topology/RefineUniform.cc	2010-08-23 23:49:58 UTC (rev 17111)
@@ -109,16 +109,24 @@
   const ALE::Obj<SieveMesh>& newSieveMesh = newMesh->sieveMesh();
   assert(!newSieveMesh.isNull());
 
+  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());
 
-  std::map<edge_type, point_type> edge2vertex;
-   
-#if 0 // DON'T KNOW NEW INTERFACE 
   newSieveMesh->setSieve(newSieve);
   ALE::MeshBuilder<Mesh>::refineTetrahedra(*mesh.sieveMesh(), * newSieveMesh,
 					   edge2vertex);
+#else
+#if 0 // DON'T KNOW NEW INTERFACE 
+  // Is arg to CellRefiner constructor the new mesh or the old mesh?
+  ALE::MeshBuilder<SieveMesh>::CellRefiner<SieveMesh,edge_type> refiner(*newSieveMesh);
+
+  ALE::MeshBuilder<SieveMesh>::refineGeneral< SieveMesh,
+    ALE::MeshBuilder<SieveMesh>::CellRefiner<SieveMesh,edge_type> >(*sieveMesh, *newSieveMesh, refiner);
 #endif
+#endif
 
   // Fix material ids
   const int numCells = sieveMesh->heightStratum(0)->size();