[cig-commits] r17282 - short/3D/PyLith/trunk/libsrc/meshio
brad at geodynamics.org
brad at geodynamics.org
Fri Oct 15 11:59:54 PDT 2010
Author: brad
Date: 2010-10-15 11:59:54 -0700 (Fri, 15 Oct 2010)
New Revision: 17282
Modified:
short/3D/PyLith/trunk/libsrc/meshio/DataWriterHDF5.cc
Log:
Changes to debugging output. Now turned off.
Modified: short/3D/PyLith/trunk/libsrc/meshio/DataWriterHDF5.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/meshio/DataWriterHDF5.cc 2010-10-15 18:56:49 UTC (rev 17281)
+++ short/3D/PyLith/trunk/libsrc/meshio/DataWriterHDF5.cc 2010-10-15 18:59:54 UTC (rev 17282)
@@ -169,6 +169,16 @@
field.createVector();
vector = field.vector();
}
+
+#if 0 // TEMPORARY DEBUGGING
+ const char* vecname = 0;
+ PetscObjectGetName((PetscObject) vector, &vecname);
+ std::cout << "NAME field: " << field.label()
+ << ", section: " << field.section()->getName()
+ << ", vec: " << vecname
+ << std::endl;
+#endif
+
// TODO: Create scatter if necessary
field.createScatter();
field.scatterSectionToVector();
@@ -209,15 +219,13 @@
vector = field.vector();
}
-#if 1 // TEMPORARY DEBUGGING
- field.view("CELL FIELD");
- const ALE::Obj<typename mesh_type::SieveMesh>& sieveMesh = field.mesh().sieveMesh();
- assert(!sieveMesh.isNull());
- const ALE::Obj<typename mesh_type::RealSection>& section = field.section();
- assert(!section.isNull());
- const ALE::Obj<typename mesh_type::SieveMesh::order_type>& globalOrder = sieveMesh->getFactory()->getGlobalOrder(sieveMesh, section->getName(), section);
- std::cout << "GLOBAL ORDER LOCAL SIZE: " << globalOrder->getLocalSize() << std::endl;
-
+#if 0 // TEMPORARY DEBUGGING
+ const char* vecname = 0;
+ PetscObjectGetName((PetscObject) vector, &vecname);
+ std::cout << "NAME field: " << field.label()
+ << ", section: " << field.section()->getName()
+ << ", vec: " << vecname
+ << std::endl;
#endif
// TODO: Create scatter only if necessary
field.createScatter();
More information about the CIG-COMMITS
mailing list