[MOAB-dev] r5500 - MOAB/trunk/tools/mbcoupler
iulian at mcs.anl.gov
iulian at mcs.anl.gov
Thu Apr 26 06:51:01 CDT 2012
Author: iulian
Date: 2012-04-26 06:51:01 -0500 (Thu, 26 Apr 2012)
New Revision: 5500
Modified:
MOAB/trunk/tools/mbcoupler/Coupler.cpp
Log:
spectral target flag and some formatting
Modified: MOAB/trunk/tools/mbcoupler/Coupler.cpp
===================================================================
--- MOAB/trunk/tools/mbcoupler/Coupler.cpp 2012-04-25 23:37:02 UTC (rev 5499)
+++ MOAB/trunk/tools/mbcoupler/Coupler.cpp 2012-04-26 11:51:01 UTC (rev 5500)
@@ -191,9 +191,9 @@
}
// repeat for target sets
spectral_sets.empty();
- // now initialize a source spectral element !
+ // now initialize a target spectral element !
_spectralTarget = new moab::Element::SpectralHex(sem_dims[0]);
- specSou = true;
+ specTar = true;
}
_ntot = sem_dims[0]*sem_dims[1]*sem_dims[2];
@@ -585,7 +585,7 @@
result = mbImpl->get_entities_by_dimension(treeiter.handle(), 3, range_leaf, false);
if(result != MB_SUCCESS) std::cout << "Problem getting leaf in a range" << std::endl;
- // loop over the range_leaf
+ // loop over the range_leaf
for(Range::iterator iter = range_leaf.begin(); iter != range_leaf.end(); iter++)
{
//test to find out in which entity the point is
@@ -620,7 +620,6 @@
spcHex->set_gl_points((double*)xval, (double*)yval, (double*)zval);
try{
-
tmp_nat_coords =spcHex->ievaluate(CartVect(xyz));
}
catch (Element::Map::EvaluationError) {
@@ -629,8 +628,6 @@
continue; // it is possible that the point is outside, so it will not converge
}
-
-
}
else
{
@@ -648,7 +645,6 @@
return result;
}
-
if (etype == MBHEX) {
Element::LinearHex hexmap(coords_vert);
try {
@@ -689,69 +685,70 @@
More information about the moab-dev
mailing list