fix a bug which made us use the -1 ("past") time level every 2nd time

step of an evolution:

Previously the code found the data pointers for the Cactus g_ij and K_ij
grid arrays once in its setup phase, then used those all through the
evolution.  The problem is... the "current time level" data pointers
change at each time step, because Cactus rotates the time-level arrays
rather than copying all the (big) grid array data.

This commit changes the coee to just lookup the Cactus variable indices
in setup, then get a fresh set of data pointers each time we look for
an apparent horizon (i.e. at each time step).


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@850 f88db872-0e4f-0410-b76b-b9085cfa78c5

1 files changed, 9 insertions, 1 deletions

diff --git a/src/gr/gr.hh b/src/gr/gr.hh
index 3f86254..0c933f4 100644
--- a/src/gr/gr.hh
+++ b/src/gr/gr.hh
@@ -58,6 +58,15 @@ struct	cactus_grid_info
 	// false ==> the Cactus gij are the physical metric
 	bool Cactus_conformal_metric;
 
+	// Cactus variable indices
+	int g_dd_11_varindex, g_dd_12_varindex, g_dd_13_varindex,
+			      g_dd_22_varindex, g_dd_23_varindex,
+						g_dd_33_varindex;
+	int K_dd_11_varindex, K_dd_12_varindex, K_dd_13_varindex,
+			      K_dd_22_varindex, K_dd_23_varindex,
+						K_dd_33_varindex;
+	int psi_varindex;		// unused if !Cactus_conformal_metric
+
 	// --> Cactus gridfn data for grid posn (0,0,0)
 	const fp* g_dd_11_data;
 	const fp* g_dd_12_data;
@@ -71,7 +80,6 @@ struct	cactus_grid_info
 	const fp* K_dd_22_data;
 	const fp* K_dd_23_data;
 	const fp* K_dd_33_data;
-
 	const fp* psi_data;		// NULL if !Cactus_conformal_metric
 	};