An additional coord_delta appeared in the local grid geometry setup. With this removed it finds horizons and so on again.

An additional coord_delta appeared in the local grid geometry setup. With this removed it finds horizons and so on again.
However, for me it seems to be an order of magnitude slower than before?!

darcs-hash:20040129113915-58737-8c9c38eabdb0523434e413305dfcb6378e4f4c2b.gz

1 files changed, 4 insertions, 3 deletions

diff --git a/Carpet/CarpetInterp/src/interp.cc b/Carpet/CarpetInterp/src/interp.cc
index 35a006aae..f0a95d3dc 100644
--- a/Carpet/CarpetInterp/src/interp.cc
+++ b/Carpet/CarpetInterp/src/interp.cc
@@ -1,4 +1,4 @@
-// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetInterp/src/interp.cc,v 1.18 2004/01/25 14:57:29 schnetter Exp $
+// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetInterp/src/interp.cc,v 1.19 2004/01/29 12:39:15 hawke Exp $
 
 #include <assert.h>
 #include <math.h>
@@ -21,7 +21,7 @@
 #include "interp.hh"
 
 extern "C" {
-  static char const * const rcsid = "$Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetInterp/src/interp.cc,v 1.18 2004/01/25 14:57:29 schnetter Exp $";
+  static char const * const rcsid = "$Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetInterp/src/interp.cc,v 1.19 2004/01/29 12:39:15 hawke Exp $";
   CCTK_FILEVERSION(Carpet_CarpetInterp_interp_cc);
 }
 
@@ -325,7 +325,8 @@ namespace CarpetInterp {
               for (int d=0; d<dim; ++d) {
                 lsh[d] = cgh->cctk_lsh[d];
                 coord_delta[d] = cgh->cctk_delta_space[d] / cgh->cctk_levfac[d];
-                coord_origin[d] = cgh->cctk_origin_space[d] + coord_delta[d] * (cgh->cctk_lbnd[d] + 1.0 * cgh->cctk_levoff[d] / cgh->cctk_levoffdenom[d]) * coord_delta[d];
+                coord_origin[d] = cgh->cctk_origin_space[d] + (cgh->cctk_lbnd[d] + 1.0 * cgh->cctk_levoff[d] / cgh->cctk_levoffdenom[d]) * coord_delta[d];
+                
               }