start adding support for sparse matrix Jacobians (doesn't work yet)

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

1 files changed, 7 insertions, 6 deletions

diff --git a/src/driver/driver.hh b/src/driver/driver.hh
index 119bafe..9747002 100644
--- a/src/driver/driver.hh
+++ b/src/driver/driver.hh
@@ -181,7 +181,7 @@ struct	BH_diagnostics
 struct	AH_info
 	{
 	patch_system* ps_ptr;
-	Jacobian* Jac_ptr;
+	Jacobian_matrix* Jac_ptr;
 
 	struct initial_guess_info initial_guess_info;
 
@@ -245,7 +245,7 @@ enum initial_guess_method
 // Newton.cc
 // returns true for success, false for failure to converge
 bool Newton_solve(patch_system& ps,
-		  Jacobian& Jac,
+		  Jacobian_matrix& Jac,
 		  const struct cactus_grid_info& cgi,
 		  const struct geometry_info& gi,
 		  const struct Jacobian_info& Jacobian_info,
@@ -262,10 +262,11 @@ void input_gridfn(patch_system& ps, int unknown_gfn,
 void output_gridfn(patch_system& ps, int unknown_gfn,
 		   const struct IO_info& IO_info, const char base_file_name[],
 		   int hn, bool print_msg_flag, int AHF_iteration = 0);
-void print_Jacobians(const patch_system& ps,
-		     const Jacobian* Jac_NP, const Jacobian* Jac_SD_FDdr,
-		     const struct IO_info& IO_info, const char base_file_name[],
-		     int hn, bool print_msg_flag, int AHF_iteration = 0);
+void output_Jacobians(const patch_system& ps,
+		      const Jacobian_matrix* Jac_NP_ptr,
+		      const Jacobian_matrix* Jac_SD_FDdr_ptr,
+		      const struct IO_info& IO_info, const char base_file_name[],
+		      int hn, bool print_msg_flag, int AHF_iteration = 0);
 FILE* setup_BH_diagnostics_output_file(const struct IO_info& IO_info,
 				       int hn, int N_horizons);
 void output_BH_diagnostics_fn(const struct BH_diagnostics& BH_diagnostics,