From 7872a7bec468f8a9c82e3a90f60951994864b33c Mon Sep 17 00:00:00 2001 From: Erik Schnetter Date: Thu, 19 Apr 2007 02:25:00 +0000 Subject: CarpetSlab: Update to recent changes darcs-hash:20070419022546-dae7b-73710f7f14788fe66b3c7e44efb01637ad780343.gz --- Carpet/CarpetSlab/src/GetHyperslab.cc | 11 +++++------ Carpet/CarpetSlab/src/slab.cc | 10 ++++------ 2 files changed, 9 insertions(+), 12 deletions(-) (limited to 'Carpet/CarpetSlab') diff --git a/Carpet/CarpetSlab/src/GetHyperslab.cc b/Carpet/CarpetSlab/src/GetHyperslab.cc index 1c9e687cc..9a4e07117 100644 --- a/Carpet/CarpetSlab/src/GetHyperslab.cc +++ b/Carpet/CarpetSlab/src/GetHyperslab.cc @@ -171,7 +171,7 @@ namespace CarpetSlab { // Create collector data object void* myhdata = rank==collect_proc ? hdata : 0; - gdata* const alldata = mydata->make_typed(-1); + gdata* const alldata = mydata->make_typed (-1, error_centered, op_sync); alldata->allocate (hextent, collect_proc, myhdata); // Done with the temporary stuff @@ -186,10 +186,9 @@ namespace CarpetSlab { mydata = (*myff)(tl, rl, component, mglevel); // Calculate overlapping extents - const bboxset myextents - = ((mydd->boxes.at(mglevel).at(rl).at(component).sync_not - | mydd->boxes.at(mglevel).at(rl).at(component).interior) - & hextent); + bboxset myextents = + mydd->boxes.at(mglevel).at(rl).at(component).interior & hextent; + myextents.normalize(); // Loop over overlapping extents for (bboxset::const_iterator ext_iter = myextents.begin(); @@ -213,7 +212,7 @@ namespace CarpetSlab { for (int proc=0; procmake_typed(-1); + tmpdata.at(proc) = mydata->make_typed (-1, error_centered, op_sync); tmpdata.at(proc)->allocate (alldata->extent(), proc, myhdata); tmpdata.at(proc)->copy_from (state.at(proc), alldata, alldata->extent()); } diff --git a/Carpet/CarpetSlab/src/slab.cc b/Carpet/CarpetSlab/src/slab.cc index 7bb364033..ba88e16a6 100644 --- a/Carpet/CarpetSlab/src/slab.cc +++ b/Carpet/CarpetSlab/src/slab.cc @@ -204,7 +204,7 @@ namespace CarpetSlab { // Create collector data object void* myhdata = rank==collect_proc ? hdata : 0; - gdata* const alldata = mydata->make_typed(-1); + gdata* const alldata = mydata->make_typed (-1, error_centered, op_sync); alldata->allocate (hextent, collect_proc, myhdata); // Done with the temporary stuff @@ -219,10 +219,8 @@ namespace CarpetSlab { mydata = (*myff)(tl, rl, component, mglevel); // Calculate overlapping extents - const bboxset myextents - = ((mydd->boxes.at(mglevel).at(rl).at(component).sync_not - | mydd->boxes.at(mglevel).at(rl).at(component).interior) - & hextent); + const bboxset myextents = + mydd->boxes.at(mglevel).at(rl).at(component).interior & hextent; // Loop over overlapping extents for (bboxset::const_iterator ext_iter = myextents.begin(); @@ -246,7 +244,7 @@ namespace CarpetSlab { for (int proc=0; procmake_typed(-1); + tmpdata.at(proc) = mydata->make_typed (-1, error_centered, op_sync); tmpdata.at(proc)->allocate (alldata->extent(), proc, myhdata); tmpdata.at(proc)->copy_from (state.at(proc), alldata, alldata->extent()); } -- cgit v1.2.3