[cig-commits] r5094 - short/3D/PyLith/branches/pylith-0.8/pylith3d/module

[knepley at geodynamics.org]

Author: knepley
Date: 2006-10-26 10:03:00 -0700 (Thu, 26 Oct 2006)
New Revision: 5094

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc
Log:
Fixed problem with BC
 - Was figuring numElements incorrectly


Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc	2006-10-26 00:28:07 UTC (rev 5093)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc	2006-10-26 17:03:00 UTC (rev 5094)
@@ -215,6 +215,7 @@
   ierr = MeshCreatePyLith(comm, 3, meshInputFile, PETSC_FALSE, (PetscTruth) interpolateMesh, &mesh);
   ierr = MeshGetMesh(mesh, m);
   m->setDebug(debugFlag);
+  int numElements = m->getTopology()->heightStratum(0, 0)->size();
   debug << journal::at(__HERE__) << "[" << rank << "]Created new PETSc Mesh for " << meshInputFile << journal::endl;
   m = ALE::New::Distribution<ALE::Mesh::topology_type>::distributeMesh(m, partitioner);
   ierr = MeshSetMesh(mesh, m);
@@ -225,7 +226,6 @@
   ALE::Mesh::foliated_section_type::patch_type patch      = 0;
   std::set<int> seen;
 
-  int numElements = m->getTopology()->heightStratum(0, 0)->size();
   ierr = MPI_Bcast(&numElements, 1, MPI_INT, 0, comm);
   boundaries->setTopology(m->getTopology());
   // Reverse order allows newer conditions to override older, as required by PyLith