[cig-commits] r14963 - short/3D/PyLith/trunk/libsrc/materials
brad at geodynamics.org
brad at geodynamics.org
Sun May 10 13:07:47 PDT 2009
Author: brad
Date: 2009-05-10 13:07:47 -0700 (Sun, 10 May 2009)
New Revision: 14963
Modified:
short/3D/PyLith/trunk/libsrc/materials/Material.cc
Log:
Removed debugging info.
Modified: short/3D/PyLith/trunk/libsrc/materials/Material.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/materials/Material.cc 2009-05-10 19:47:51 UTC (rev 14962)
+++ short/3D/PyLith/trunk/libsrc/materials/Material.cc 2009-05-10 20:07:47 UTC (rev 14963)
@@ -254,7 +254,6 @@
int stateVarIndex = -1;
_findField(&propertyIndex, &stateVarIndex, name);
- std::cout << "AA" << std::endl;
// Get cell information
const ALE::Obj<SieveMesh>& sieveMesh = field->mesh().sieveMesh();
assert(!sieveMesh.isNull());
@@ -264,11 +263,9 @@
const SieveMesh::label_sequence::iterator cellsBegin = cells->begin();
const SieveMesh::label_sequence::iterator cellsEnd = cells->end();
- std::cout << "BB" << std::endl;
topology::FieldBase::VectorFieldEnum fieldType = topology::FieldBase::OTHER;
if (propertyIndex >= 0) { // If field is a property
- std::cout << "CC" << std::endl;
int propOffset = 0;
const string_vector& properties = _metadata.properties();
assert(propertyIndex < properties.size());
@@ -277,37 +274,30 @@
_metadata.fiberDim(properties[i].c_str(), Metadata::PROPERTY);
const int fiberDim = _metadata.fiberDim(name, Metadata::PROPERTY);
- std::cout << "DD" << std::endl;
// :TODO: Get scale information
// Get properties section
- std::cout << "EE" << std::endl;
const ALE::Obj<RealSection>& propertiesSection = _properties->section();
assert(!propertiesSection.isNull());
const int totalPropsFiberDimLocal = (cells->size() > 0) ?
propertiesSection->getFiberDimension(*cells->begin()) : 0;
- std::cout << "FF" << std::endl;
int totalPropsFiberDim = 0;
MPI_Allreduce((void *) &totalPropsFiberDimLocal,
(void *) &totalPropsFiberDim, 1,
MPI_INT, MPI_MAX, field->mesh().comm());
assert(totalPropsFiberDim > 0);
- std::cout << "GG" << std::endl;
const int numPropsQuadPt = _numPropsQuadPt;
const int numQuadPts = totalPropsFiberDim / numPropsQuadPt;
assert(totalPropsFiberDim == numQuadPts * numPropsQuadPt);
const int totalFiberDim = numQuadPts * fiberDim;
- std::cout << "HH" << std::endl;
// Allocate buffer for property field if necessary.
- std::cout << "II" << std::endl;
const ALE::Obj<RealSection>& fieldSection = field->section();
bool useCurrentField = !fieldSection.isNull();
if (!fieldSection.isNull()) {
// check fiber dimension
const int totalFiberDimCurrentLocal = (cells->size() > 0) ?
fieldSection->getFiberDimension(*cells->begin()) : 0;
- std::cout << "JJ" << std::endl;
int totalFiberDimCurrent = 0;
MPI_Allreduce((void *) &totalFiberDimCurrentLocal,
(void *) &totalFiberDimCurrent, 1,
@@ -315,16 +305,13 @@
assert(totalFiberDimCurrent > 0);
useCurrentField = totalFiberDim == totalFiberDimCurrent;
} // if
- std::cout << "KK" << std::endl;
if (!useCurrentField) {
field->newSection(cells, totalFiberDim);
field->allocate();
} // if
- std::cout << "NN" << std::endl;
assert(!fieldSection.isNull());
field->label(name);
fieldType = _metadata.fieldType(name, Metadata::PROPERTY);
- std::cout << "OO" << std::endl;
// Buffer for property at cell's quadrature points
double_array fieldCell(numQuadPts*fiberDim);
@@ -382,7 +369,6 @@
// check fiber dimension
const int totalFiberDimCurrentLocal = (cells->size() > 0) ?
fieldSection->getFiberDimension(*cells->begin()) : 0;
- std::cout << "JJ" << std::endl;
int totalFiberDimCurrent = 0;
MPI_Allreduce((void *) &totalFiberDimCurrentLocal,
(void *) &totalFiberDimCurrent, 1,
More information about the CIG-COMMITS
mailing list