add some code for the d/dr term in the Jacobian -- not finished yet

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

1 files changed, 26 insertions, 8 deletions

diff --git a/src/gr/AHFinderDirect.hh b/src/gr/AHFinderDirect.hh
index 42b77c9..a8305d2 100644
--- a/src/gr/AHFinderDirect.hh
+++ b/src/gr/AHFinderDirect.hh
@@ -55,6 +55,26 @@ struct	cactus_grid_info
 	const fp* K_dd_33_data;
 	};
 
+//
+// This struct holds the various parameters used in computing the Jacobian
+// matrix.
+//
+struct	Jacobian_info
+	{
+	fp perturbation_amplitude;
+	};
+
+//
+// This struct holds the various parameters used in solving the H(h) = 0
+// equations.
+//
+struct	solver_info
+	{
+	// stuff for Newton_solve()
+	int max_Newton_iterations;
+	fp H_norm_for_convergence;
+	};
+
 //******************************************************************************
 
 //
@@ -70,19 +90,17 @@ void horizon_function(patch_system& ps,
 		      const struct cactus_grid_info& cgi,
 		      const struct geometry_interpolator_info& gii,
 		      bool Jacobian_flag,
-		      jtutil::norm<fp>& H_norms,
+		      jtutil::norm<fp>* H_norms_ptr,
 		      bool msg_flag = true);
 
 // horizon_Jacobian.cc
 Jacobian& create_Jacobian(patch_system& ps,
 			  const char Jacobian_type[]);
-void horizon_Jacobian_SD(patch_system& ps,
-			 Jacobian& Jac);
-void horizon_Jacobian_NP(patch_system& ps,
-			 const struct cactus_grid_info& cgi,
-			 const struct geometry_interpolator_info& ii,
-			 Jacobian& Jac,
-			 fp perturbation_amplitude);
+void horizon_Jacobian(patch_system& ps,
+		      const struct cactus_grid_info& cgi,
+		      const struct geometry_interpolator_info& gii,
+		      const struct Jacobian_info& Jac_info,
+		      Jacobian& Jac);
 
 // Newton.cc
 void Newton_solve(patch_system& ps,