Dont use initial guess origins from parameters when tracking from gridscalar

Problem: My horizon will appear sometimes during the simulation at some position. I don't know this position in advance, so I set origin_* and AHFinderDirect::initial_guesscoord_sphere*_center to zero initially in my par-file. However, I have a grid scalar which tracks the coordinate location where the horizon will eventually appear.
When the horizon finder starts to search for a horizon, it will first setup the
coordinate ellipsoid. This routine is executed *before* the new origin is set
from the grid scalar.
The tracking occurs in the routine Newton(...). The ellipsoid is set before
Newton(...) gets executed. Hence, the ellipsoid uses the value that got set via
parameters (which would be zero in my case).

Patch by Christian Reisswig



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

1 files changed, 9 insertions, 0 deletions

diff --git a/src/driver/driver.hh b/src/driver/driver.hh
index 6569ba8..c998ce9 100644
--- a/src/driver/driver.hh
+++ b/src/driver/driver.hh
@@ -431,6 +431,10 @@ void setup_initial_guess(patch_system& ps,
 enum initial_guess_method
   decode_initial_guess_method(const char initial_guess_method_string[]);
 
+
+void set_initial_guess_parameters(struct AH_data& AH_data, const int hn,
+                                  const fp origin_x, const fp origin_y, const fp origin_z);
+
 // Newton.cc
 // returns true for success, false for failure to converge
 void Newton(const cGH* GH,
@@ -447,6 +451,11 @@ void Newton(const cGH* GH,
 	    const struct verbose_info& verbose_info,
 	    struct iteration_status_buffers& isb);
 
+// Tracks coordinate origin
+void track_origin(const cGH* const cctkGH, patch_system& ps, 
+                  struct AH_data* const AH_data_ptr, 
+                  const int hn, const bool print_algorithm_highlights);
+
 // io.cc
 void input_gridfn(patch_system& ps, int unknown_gfn,
 		  const struct IO_info& IO_info, const char base_file_name[],