Fix bug in three timelevel initialisation:

Fix bug in three timelevel initialisation:
Change order of time level flipping and restricting.

ggf.cc: Remove commented out code.

darcs-hash:20030810105209-07bb3-40c12c5fb838cd56a99f26c8ca71529bfe11a900.gz

1 files changed, 4 insertions, 6 deletions

diff --git a/Carpet/Carpet/src/Restrict.cc b/Carpet/Carpet/src/Restrict.cc
index e392af3ca..1467cdd58 100644
--- a/Carpet/Carpet/src/Restrict.cc
+++ b/Carpet/Carpet/src/Restrict.cc
@@ -10,7 +10,7 @@
 #include "carpet.hh"
 
 extern "C" {
-  static const char* rcsid = "$Header: /home/eschnett/C/carpet/Carpet/Carpet/Carpet/src/Restrict.cc,v 1.19 2003/07/20 21:03:43 schnetter Exp $";
+  static const char* rcsid = "$Header: /home/eschnett/C/carpet/Carpet/Carpet/Carpet/src/Restrict.cc,v 1.20 2003/08/10 12:52:09 schnetter Exp $";
   CCTK_FILEVERSION(Carpet_Carpet_Restrict_cc);
 }
 
@@ -40,11 +40,9 @@ namespace Carpet {
             // use background time here (which may not be modified by
             // the user)
             const CCTK_REAL time = tt->time (tl, reflevel, mglevel);
-            if (tl==0) {
-              const CCTK_REAL time1 = tt->time (tl, reflevel, mglevel);
-              const CCTK_REAL time2 = cgh->cctk_time / cgh->cctk_delta_time;
-              assert (fabs((time1 - time2) / (fabs(time1) + fabs(time2) + fabs(cgh->cctk_delta_time))) < 1e-12);
-            }
+            const CCTK_REAL time1 = tt->time (0, reflevel, mglevel) / tt->get_delta (0, 0);
+            const CCTK_REAL time2 = cgh->cctk_time / cgh->cctk_delta_time;
+            assert (fabs((time1 - time2) / (fabs(time1) + fabs(time2) + fabs(cgh->cctk_delta_time))) < 1e-12);
             
             if (reflevel < hh->reflevels()-1) {