[cig-commits] r18400 - short/3D/PyLith/trunk/libsrc/meshio

[knepley at geodynamics.org]

Author: knepley
Date: 2011-05-20 12:01:26 -0700 (Fri, 20 May 2011)
New Revision: 18400

Modified:
   short/3D/PyLith/trunk/libsrc/meshio/MeshBuilder.cc
Log:
Add orientation section so we can build interpolated meshes


Modified: short/3D/PyLith/trunk/libsrc/meshio/MeshBuilder.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/meshio/MeshBuilder.cc	2011-05-20 16:20:15 UTC (rev 18399)
+++ short/3D/PyLith/trunk/libsrc/meshio/MeshBuilder.cc	2011-05-20 19:01:26 UTC (rev 18400)
@@ -120,13 +120,17 @@
       // Same old thing
       ALE::Obj<SieveFlexMesh::sieve_type> s =
 	new SieveFlexMesh::sieve_type(sieve->comm(), sieve->debug());
+      ALE::Obj<SieveFlexMesh::arrow_section_type> orientation = new SieveFlexMesh::arrow_section_type(sieve->comm(), sieve->debug());
 
+      s->setDebug(2);
       ALE::SieveBuilder<SieveFlexMesh>::buildTopology(s, meshDim, 
-                                                  numCells, 
-                                                  const_cast<int*>(&cells[0]), 
-                                                  numVertices, 
-                                                  interpolate,
-                                                  numCorners);
+                                                      numCells, 
+                                                      const_cast<int*>(&cells[0]), 
+                                                      numVertices, 
+                                                      interpolate,
+                                                      numCorners,
+                                                      -1,
+                                                      orientation);
       std::map<SieveFlexMesh::point_type,SieveFlexMesh::point_type> renumbering;
       ALE::ISieveConverter::convertSieve(*s, *sieve, renumbering);
     } // if/else