diff options
| author | shawley <> | 2002-06-07 15:25:00 +0000 |
|---|---|---|
| committer | shawley <> | 2002-06-07 15:25:00 +0000 |
| commit | f33e05cc2427b48c4717ba55bb867fc0c11b0389 (patch) | |
| tree | 58aaa2bdc0d0c592f24ba2c17abad61af7635a9b /Carpet | |
| parent | 7f2e2ec109fc496f781c39ce5633c9ca6d32c8cb (diff) | |
added calls to tt-> set_delta..., also setting current time explicitly to
added calls to tt-> set_delta..., also setting current time explicitly to
zero after "last flip"
darcs-hash:20020607152509-e415b-f0e4e568c0325ebc849292f6637c021cb2e7cacf.gz
Diffstat (limited to 'Carpet')
| -rw-r--r-- | Carpet/Carpet/src/Initialise.cc | 10 |
1 files changed, 6 insertions, 4 deletions
diff --git a/Carpet/Carpet/src/Initialise.cc b/Carpet/Carpet/src/Initialise.cc index 62e271938..9d1b5d129 100644 --- a/Carpet/Carpet/src/Initialise.cc +++ b/Carpet/Carpet/src/Initialise.cc @@ -9,7 +9,7 @@ #include "carpet.hh" -static const char* rcsid = "$Header: /home/eschnett/C/carpet/Carpet/Carpet/Carpet/src/Initialise.cc,v 1.15 2002/06/07 16:44:42 schnetter Exp $"; +static const char* rcsid = "$Header: /home/eschnett/C/carpet/Carpet/Carpet/Carpet/src/Initialise.cc,v 1.16 2002/06/07 17:25:09 shawley Exp $"; CCTK_FILEVERSION(Carpet_Initialise_cc) @@ -117,7 +117,8 @@ namespace Carpet { // erik: what about arrays? FlipTimeLevels(cgh); for (int rl=0; rl<hh->reflevels(); ++rl) { - tt->set_time (rl, mglevel, 2*tt->get_delta (rl, mglevel) - tt->get_time (rl, mglevel)); + tt->set_time (rl, mglevel, -1*tt->get_time (rl, mglevel) ); + tt->set_delta (rl, mglevel,-1*tt->get_delta (rl, mglevel)); } // Keep track of which direction (in time) we're integrating time_dir *= -1; @@ -149,7 +150,8 @@ namespace Carpet { BEGIN_MGLEVEL_LOOP(cgh) { FlipTimeLevels(cgh); for (int rl=0; rl<hh->reflevels(); ++rl) { - tt->set_time (rl, mglevel, 2*tt->get_delta (rl, mglevel) - tt->get_time (rl, mglevel)); + tt->set_time (rl, mglevel, -1*tt->get_time (rl, mglevel) ); + tt->set_delta (rl, mglevel,-1*tt->get_delta (rl, mglevel)); } } END_MGLEVEL_LOOP(cgh); time_dir *= -1; @@ -184,7 +186,7 @@ namespace Carpet { BEGIN_MGLEVEL_LOOP(cgh) { FlipTimeLevels(cgh); for (int rl=0; rl<hh->reflevels(); ++rl) { - tt->set_time (rl, mglevel, 2*tt->get_delta (rl, mglevel) - tt->get_time (rl, mglevel)); + tt->set_time (rl, mglevel, 0); } } END_MGLEVEL_LOOP(cgh); time_dir *= -1; |