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

[knepley at geodynamics.org]

Author: knepley
Date: 2006-08-12 13:45:52 -0700 (Sat, 12 Aug 2006)
New Revision: 4277

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc
Log:
Now split nodes seem to work


Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc	2006-08-12 15:42:37 UTC (rev 4276)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc	2006-08-12 20:45:52 UTC (rev 4277)
@@ -177,7 +177,7 @@
 
 // Process mesh
 
-PetscErrorCode MeshView_Sieve(ALE::Obj<ALE::Mesh> mesh, PetscViewer viewer);
+PetscErrorCode MeshView_Sieve(const ALE::Obj<ALE::Mesh>& mesh, PetscViewer viewer);
 
 char pypylith3d_processMesh__doc__[] = "";
 char pypylith3d_processMesh__name__[] = "processMesh";
@@ -206,11 +206,12 @@
   PetscInt         *boundaryVertices;
   PetscScalar      *boundaryValues;
   PetscInt          numBoundaryVertices, numBoundaryComponents;
+  int               debugFlag = 0;
   PetscErrorCode    ierr;
 
   ierr = MPI_Comm_rank(comm, &rank);
   sprintf(meshOutputFile, "%s.%d", meshInputFile, rank);
-  mesh = ALE::PyLith::Builder::readMesh(comm, 3, meshInputFile, false, (bool) interpolateMesh, 1);
+  mesh = ALE::PyLith::Builder::readMesh(comm, 3, meshInputFile, false, (bool) interpolateMesh, debugFlag);
   int numElements = mesh->getTopologyNew()->heightStratum(0, 0)->size();
   ierr = MPI_Bcast(&numElements, 1, MPI_INT, 0, comm);
   debug << journal::at(__HERE__) << "[" << rank << "]Created new PETSc Mesh for " << meshInputFile << journal::endl;