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

[knepley at geodynamics.org]

Author: knepley
Date: 2010-01-27 12:35:42 -0800 (Wed, 27 Jan 2010)
New Revision: 16187

Modified:
   short/3D/PyLith/trunk/libsrc/meshio/MeshIOSieve.cc
Log:
Put in Mesh serialization routines


Modified: short/3D/PyLith/trunk/libsrc/meshio/MeshIOSieve.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/meshio/MeshIOSieve.cc	2010-01-27 01:44:00 UTC (rev 16186)
+++ short/3D/PyLith/trunk/libsrc/meshio/MeshIOSieve.cc	2010-01-27 20:35:42 UTC (rev 16187)
@@ -48,17 +48,14 @@
 void
 pylith::meshio::MeshIOSieve::_read(void)
 { // _read
-  MPI_Comm comm = _mesh->comm();
-  int rank = 0;
-  int meshDim = 0;
-  int spaceDim = 0;
-  int numVertices = 0;
-  int numCells = 0;
-  int numCorners = 0;
+  _mesh->createSieveMesh(3);
+  const ALE::Obj<topology::Mesh::SieveMesh>& sieveMesh = _mesh->sieveMesh();
+  assert(!sieveMesh.isNull());
 
-
-  // :TODO: STUFF GOES HERE
-  assert(false);
+  ALE::Obj<topology::Mesh::SieveMesh::sieve_type> sieve = 
+    new topology::Mesh::SieveMesh::sieve_type(_mesh->comm());
+  sieveMesh->setSieve(sieve);
+  ALE::MeshSerializer::loadMesh(_filename, *sieveMesh);
 } // read
 
 // ----------------------------------------------------------------------
@@ -66,9 +63,7 @@
 void
 pylith::meshio::MeshIOSieve::_write(void) const
 { // write
-
   ALE::MeshSerializer::writeMesh(_filename, *_mesh->sieveMesh());
-
 } // write