various changes including d/dr terms in Jacobian by numerical perturbation,

tweak I/O parameters,
move printing Jacobian out of Jacobian class into test driver,
drop unused array BLAS routines in jtutil::


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

1 files changed, 11 insertions, 5 deletions

diff --git a/src/gr/AHFinderDirect.hh b/src/gr/AHFinderDirect.hh
index a8305d2..ac79f14 100644
--- a/src/gr/AHFinderDirect.hh
+++ b/src/gr/AHFinderDirect.hh
@@ -73,6 +73,7 @@ struct	solver_info
 	// stuff for Newton_solve()
 	int max_Newton_iterations;
 	fp H_norm_for_convergence;
+	fp Delta_h_norm_for_convergence;
 	};
 
 //******************************************************************************
@@ -101,12 +102,17 @@ void horizon_Jacobian(patch_system& ps,
 		      const struct geometry_interpolator_info& gii,
 		      const struct Jacobian_info& Jac_info,
 		      Jacobian& Jac);
+void horizon_Jacobian_NP(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,
+// return true for success, false for failure to converge
+bool Newton_solve(patch_system& ps,
 		  const struct cactus_grid_info& cgi,
 		  const struct geometry_interpolator_info& gii,
-		  const char Jacobian_type[],
-		  fp perturbation_amplitude,
-		  int max_Newton_iterations,
-		  fp H_norm_for_convergence);
+		  const struct Jacobian_info& Jac_info,
+		  const struct solver_info& solver_info,
+		  Jacobian& Jac);