diff options
| author | Erik Schnetter <schnetter@aei.mpg.de> | 2005-08-09 08:49:00 +0000 |
|---|---|---|
| committer | Erik Schnetter <schnetter@aei.mpg.de> | 2005-08-09 08:49:00 +0000 |
| commit | c7dc46fa15685b2dbeaf9fe92fb207bb5911d313 (patch) | |
| tree | 35e07770a6e6e972605bd7f6326b28d79a610b2f /Carpet/Carpet/src/Comm.cc | |
| parent | 20491415a0f46d48ba0471bf9aeefe59668f59a9 (diff) | |
Carpet: Correct problems with negative time steps
When converting from physicsl time steps (delta_time) to Carpet's time
steps (th::delta_time), do not divide by abs(delta_time), but by
delta_time instead.
When evolving backwards in time while initialising three time levels,
change the sign of the physical time steps only, not the sign of
Carpet's time steps.
darcs-hash:20050809084940-891bb-2517b1568696c71278a98db6261515817a90247a.gz
Diffstat (limited to 'Carpet/Carpet/src/Comm.cc')
| -rw-r--r-- | Carpet/Carpet/src/Comm.cc | 3 |
1 files changed, 1 insertions, 2 deletions
diff --git a/Carpet/Carpet/src/Comm.cc b/Carpet/Carpet/src/Comm.cc index 12d7e1394..0e0965ae4 100644 --- a/Carpet/Carpet/src/Comm.cc +++ b/Carpet/Carpet/src/Comm.cc @@ -1,5 +1,4 @@ #include <cassert> -#include <cmath> #include <cstdlib> #include "cctk.h" @@ -140,7 +139,7 @@ namespace Carpet { // use the current time here (which may be modified by the user) const CCTK_REAL time - = (cctkGH->cctk_time - initial_time) / abs(delta_time); + = (cctkGH->cctk_time - initial_time) / delta_time; for (comm_state state(groups.vartype); ! state.done(); state.step()) { for (int group = 0; group < groups.members.size(); ++group) { |