diff options
author | Erik Schnetter <schnetter@aei.mpg.de> | 2004-09-18 13:30:00 +0000 |
---|---|---|
committer | Erik Schnetter <schnetter@aei.mpg.de> | 2004-09-18 13:30:00 +0000 |
commit | 837059042828992bd744de5a574e5c5237fdd7cb (patch) | |
tree | f5a071cf315dfc05c5e50fcaba42169bc5ecc4c1 /Carpet/CarpetRegrid/src | |
parent | 6bd1e3636aad48e92b6ed7fb438443beaaa54608 (diff) |
Correct spacing setup on nonzero convergence levels
Do not take basemglevel into account when checking the spacing while
translating the refined regions coordinates into indices.
darcs-hash:20040918133001-891bb-14797eb044358682cd349b54c308c20f405e7b6e.gz
Diffstat (limited to 'Carpet/CarpetRegrid/src')
-rw-r--r-- | Carpet/CarpetRegrid/src/manualcoordinatelist.cc | 4 |
1 files changed, 3 insertions, 1 deletions
diff --git a/Carpet/CarpetRegrid/src/manualcoordinatelist.cc b/Carpet/CarpetRegrid/src/manualcoordinatelist.cc index fc9961602..361c979a3 100644 --- a/Carpet/CarpetRegrid/src/manualcoordinatelist.cc +++ b/Carpet/CarpetRegrid/src/manualcoordinatelist.cc @@ -158,7 +158,9 @@ namespace CarpetRegrid { bbvect const & ob = newobss.at(rl-1).at(c); // TODO: // assert (domain_from_coordbase); - rvect const spacing = base_spacing * pow(CCTK_REAL(mgfact), basemglevel) / ipow(reffact, rl); + // TODO: why can basemglevel not be used here? + // rvect const spacing = base_spacing * pow(CCTK_REAL(mgfact), basemglevel) / ipow(reffact, rl); + rvect const spacing = base_spacing / ipow(reffact, rl); assert (all(abs(ext.stride() - spacing) < spacing * 1.0e-10)); ManualCoordinates_OneLevel (cctkGH, hh, rl, refinement_levels, |