* add support for sparse-matrix Jacobians ==> works!

* change default in param.ccl to use this
* change default in src/include/config.h to default to no longer
  link in LAPACK routines
* update documentation


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

1 files changed, 16 insertions, 7 deletions

diff --git a/src/include/config.h b/src/include/config.h
index 492da5a..cf35c7a 100644
--- a/src/include/config.h
+++ b/src/include/config.h
@@ -47,6 +47,9 @@ typedef CCTK_INT integer;
 #error "don't know fp datatype!"
 #endif
 
+
+/**************************************/
+
 /*
  * The angular finite differencing in our multipatch system can be
  * either 2nd order or 4th order, configurable here.
@@ -58,12 +61,18 @@ typedef CCTK_INT integer;
 #define FINITE_DIFF_ORDER	4
 #endif
 
+/**************************************/
+
 /*
- * What types of Jacobian matrix storage do we want to compile in
- * support for?  N.b. each of these requires linking with the corresponding
- * linear-solver library; see "../elliptic/Jacobian.hh" for details
- * on the storage formats, or "../make.configuration.defn" (if it exists)
- * for details on the libraries.
+ * What types of Jacobian matrix storage and linear solvers do we want to
+ * compile in support for?  Note that each of these requires linking with
+ * the corresponding linear-solver library; see "../elliptic/Jacobian.hh"
+ * for details on the storage formats, or "../make.configuration.defn"
+ * (if it exists) for details on the libraries.
  */
-#define HAVE_DENSE_JACOBIAN
-#define HAVE_ROW_SPARSE_JACOBIAN
+
+/* store as (Fortran) dense matrix, solve with LAPACK */
+#undef HAVE_DENSE_JACOBIAN__LAPACK
+
+/* store as row-oriented sparse matrix, solve with ILUCG */
+#define HAVE_ROW_SPARSE_JACOBIAN__ILUCG