diff options
author | Erik Schnetter <schnetter@cct.lsu.edu> | 2005-11-19 21:21:00 +0000 |
---|---|---|
committer | Erik Schnetter <schnetter@cct.lsu.edu> | 2005-11-19 21:21:00 +0000 |
commit | a8d081b00bc4c07166f6827b48026a792acb3dee (patch) | |
tree | 2094122167ad2ce619f8abb9e7429a128276300c /Carpet/Carpet/src/modes.cc | |
parent | 0e0d0a99b6966b6ff900d754eea6a43af36ffe24 (diff) |
UNDO: Carpet: Store grid spacing per map instead of globally
Store the coarse grid spacing per individual map instead of globally.
Set and unset the grid spacing when entering or leaving singlemap
mode, not global mode. This breaks code that tries to access the
coarse grid spacing in the cGH in global mode.
darcs-hash:20051119212129-dae7b-ca0ccee9a08c988e44b9024b69e17da2d78d8cba.gz
Diffstat (limited to 'Carpet/Carpet/src/modes.cc')
-rw-r--r-- | Carpet/Carpet/src/modes.cc | 40 |
1 files changed, 12 insertions, 28 deletions
diff --git a/Carpet/Carpet/src/modes.cc b/Carpet/Carpet/src/modes.cc index 7e63816b3..9669097f5 100644 --- a/Carpet/Carpet/src/modes.cc +++ b/Carpet/Carpet/src/modes.cc @@ -75,13 +75,12 @@ namespace Carpet { // TODO: this could also just be "mglevel" instead cgh->cctk_convlevel = basemglevel + mglevel; - // Set time delta + // Set time and space delta cgh->cctk_delta_time = delta_time * mglevelfact; -// // Set space delta -// for (int d=0; d<dim; ++d) { -// cgh->cctk_origin_space[d] = origin_space.at(mglevel)[d]; -// cgh->cctk_delta_space[d] = delta_space[d] * mglevelfact; -// } + for (int d=0; d<dim; ++d) { + cgh->cctk_origin_space[d] = origin_space.at(mglevel)[d]; + cgh->cctk_delta_space[d] = delta_space[d] * mglevelfact; + } // Set array information for (int group=0; group<CCTK_NumGroups(); ++group) { @@ -162,16 +161,15 @@ namespace Carpet { if (mglevel == -1) return; // early return - // Save and unset time delta + // Save and unset time and space delta delta_time = cgh->cctk_delta_time / mglevelfact; cgh->cctk_delta_time = 0.0; -// // Save and unset space delta -// for (int d=0; d<dim; ++d) { -// origin_space.at(mglevel)[d] = cgh->cctk_origin_space[d]; -// delta_space[d] = cgh->cctk_delta_space[d] / mglevelfact; -// cgh->cctk_origin_space[d] = -424242.0; -// cgh->cctk_delta_space[d] = 0.0; -// } + for (int d=0; d<dim; ++d) { + origin_space.at(mglevel)[d] = cgh->cctk_origin_space[d]; + delta_space[d] = cgh->cctk_delta_space[d] / mglevelfact; + cgh->cctk_origin_space[d] = -424242.0; + cgh->cctk_delta_space[d] = 0.0; + } // Set array information for (int group=0; group<CCTK_NumGroups(); ++group) { @@ -293,12 +291,6 @@ namespace Carpet { carpetGH.map = map = m; - // Set space delta - for (int d=0; d<dim; ++d) { - cgh->cctk_origin_space[d] = origin_space.at(map).at(mglevel)[d]; - cgh->cctk_delta_space[d] = delta_space.at(map)[d] * mglevelfact; - } - // Set grid shape const ibbox& coarseext = vdd.at(map)->bases.at(mglevel).at(0 ).exterior; const ibbox& baseext = vdd.at(map)->bases.at(mglevel).at(reflevel).exterior; @@ -331,14 +323,6 @@ namespace Carpet { if (map == -1) return; // early return - // Save and unset space delta - for (int d=0; d<dim; ++d) { - origin_space.at(map).at(mglevel)[d] = cgh->cctk_origin_space[d]; - delta_space.at(map)[d] = cgh->cctk_delta_space[d] / mglevelfact; - cgh->cctk_origin_space[d] = -424242.0; - cgh->cctk_delta_space[d] = 0.0; - } - // Unset grid shape ivect::ref(cgh->cctk_levoff) = deadbeef; ivect::ref(cgh->cctk_levoffdenom) = 0; |