global: Add varying refinement factors

Add support for varying refinement factors.  The spatial refinement
factors can be different in different directions, can be different
from the time refinement factor, and can be different on each level.
(However, the underlying spatial transport operators do currently not
handle any factors except two.)

darcs-hash:20050501205010-891bb-8d3a74abaad55ee6c77ef18d51fca2a2b69740de.gz

1 files changed, 6 insertions, 4 deletions

diff --git a/Carpet/Carpet/src/modes.cc b/Carpet/Carpet/src/modes.cc
index cfa3728b9..7c2a922c3 100644
--- a/Carpet/Carpet/src/modes.cc
+++ b/Carpet/Carpet/src/modes.cc
@@ -236,9 +236,10 @@ namespace Carpet {
     Checkpoint ("Entering level mode");
     
     reflevel = rl;
-    reflevelfact = ipow(reffact, reflevel);
-    ivect::ref(cgh->cctk_levfac) = reflevelfact;
-    cgh->cctk_timefac = reflevelfact;
+    timereflevelfact = timereffacts.at (reflevel);
+    spacereflevelfact = spacereffacts.at (reflevel);
+    ivect::ref(cgh->cctk_levfac) = spacereflevelfact;
+    cgh->cctk_timefac = timereflevelfact;
     
     // Set current time
     assert (mglevel>=0 and mglevel<(int)leveltimes.size());
@@ -272,7 +273,8 @@ namespace Carpet {
     }
     
     reflevel = -1;
-    reflevelfact = -deadbeef;
+    timereflevelfact = -deadbeef;
+    spacereflevelfact = ivect(-deadbeef);
     ivect::ref(cgh->cctk_levfac) = -deadbeef;
     cgh->cctk_timefac = -deadbeef;