diff options
author | schnetter <> | 2004-05-29 09:36:00 +0000 |
---|---|---|
committer | schnetter <> | 2004-05-29 09:36:00 +0000 |
commit | bd136c754f465d09935dc58ce4624dd671441d9c (patch) | |
tree | 86ff04d836bdaaf2236fa3dff134c43ef1ab4e5c /Carpet/CarpetLib/src | |
parent | a400ba67d672bcf97c9524094b3f651261dd7437 (diff) |
Copy the whole exterior during regridding instead of relying on
Copy the whole exterior during regridding instead of relying on
boundary prolongation.
darcs-hash:20040529093622-07bb3-743ef08b661877e18b6772d31ecb6e18134773da.gz
Diffstat (limited to 'Carpet/CarpetLib/src')
-rw-r--r-- | Carpet/CarpetLib/src/ggf.cc | 9 |
1 files changed, 6 insertions, 3 deletions
diff --git a/Carpet/CarpetLib/src/ggf.cc b/Carpet/CarpetLib/src/ggf.cc index 282c9d524..466bb0689 100644 --- a/Carpet/CarpetLib/src/ggf.cc +++ b/Carpet/CarpetLib/src/ggf.cc @@ -1,4 +1,4 @@ -// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/ggf.cc,v 1.43 2004/05/21 18:13:41 schnetter Exp $ +// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/ggf.cc,v 1.44 2004/05/29 11:36:22 schnetter Exp $ #include <assert.h> #include <math.h> @@ -111,8 +111,9 @@ void ggf<D>::recompose_fill (comm_state<D>& state, const int rl) for (int tl=tmin; tl<=tmax; ++tl) { // Find out which regions need to be prolongated + // (Copy the exterior because some variables are not prolongated) // TODO: do this once in the dh instead of for each variable here - ibset work (d.boxes.at(rl).at(c).at(ml).interior); + ibset work (d.boxes.at(rl).at(c).at(ml).exterior); // Copy from old storage, if possible // TODO: copy only from interior regions? @@ -171,7 +172,7 @@ void ggf<D>::recompose_fill (comm_state<D>& state, const int rl) } // if rl // Note that work need not be empty here; in this case, not - // everything could be initialised. This is okay on + // everything could be initialised. This is okay on outer // boundaries. // TODO: check this. @@ -228,6 +229,8 @@ void ggf<D>::recompose_sync (comm_state<D>& state, const int rl) } // for c } + + // Cycle the time levels by rotating the data sets template<int D> void ggf<D>::cycle (int rl, int c, int ml) { |