diff options
author | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2003-06-02 18:12:58 +0000 |
---|---|---|
committer | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2003-06-02 18:12:58 +0000 |
commit | 75633824bcba3aa9f9b380f504143636ec9ab3e2 (patch) | |
tree | 9a5109d1c19bef859363996e0d7fbaf273c84338 | |
parent | f2d4c13de1d3d35ab6c2073648949d98fecbec72 (diff) |
define fp datatype roundoff threshold,
add definition for UMFPACK sparse solver (currently *not* compiled by default)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1085 f88db872-0e4f-0410-b76b-b9085cfa78c5
-rw-r--r-- | src/include/config.h | 31 |
1 files changed, 28 insertions, 3 deletions
diff --git a/src/include/config.h b/src/include/config.h index f5312fd..2655296 100644 --- a/src/include/config.h +++ b/src/include/config.h @@ -17,6 +17,8 @@ */ /******************************************************************************/ +/******************************************************************************/ +/******************************************************************************/ /* * definitions for C/C++ only @@ -30,6 +32,8 @@ typedef CCTK_INT integer; #endif /* CCODE */ /******************************************************************************/ +/******************************************************************************/ +/******************************************************************************/ /* * definitions for C++, C, and Fortran @@ -40,7 +44,11 @@ typedef CCTK_INT integer; /* with exit codes defined in "config.h" */ #define FATAL_ERROR (-1) -/**************************************/ +/******************************************************************************/ + +/* + * ***** definitions for integer and floating point datatypes ***** + */ /* * FIXME: this assumes fp == C 'double' @@ -51,13 +59,21 @@ typedef CCTK_INT integer; #if defined(FP_IS_FLOAT) #define FP_SCANF_FORMAT "%f" +#define FP_EPSILON FLT_EPSILON /* from <float.h> */ #elif defined(FP_IS_DOUBLE) #define FP_SCANF_FORMAT "%lf" +#define FP_EPSILON DBL_EPSILON /* from <float.h> */ #else #error "don't know fp datatype!" #endif -/**************************************/ +/* + * FIXME: this should be set somehow at configure time + */ +#define FORTRAN_INTEGER_IS_INT +#undef FORTRAN_INTEGER_IS_LONG + +/******************************************************************************/ /* * The angular finite differencing in our multipatch system can be @@ -70,7 +86,7 @@ typedef CCTK_INT integer; #define FINITE_DIFF_ORDER 4 #endif -/**************************************/ +/******************************************************************************/ /* * What types of Jacobian matrix storage and linear solvers do we want to @@ -86,6 +102,11 @@ typedef CCTK_INT integer; /* store as row-oriented sparse matrix, solve with ILUCG */ #define HAVE_ROW_SPARSE_JACOBIAN__ILUCG +/* store as row-oriented sparse matrix, solve with UMFPACK */ +#undef HAVE_ROW_SPARSE_JACOBIAN__UMFPACK + +/******************************************************************************/ +/******************************************************************************/ /******************************************************************************/ /* @@ -101,4 +122,8 @@ typedef CCTK_INT integer; #define HAVE_ROW_SPARSE_JACOBIAN #endif +#ifdef HAVE_ROW_SPARSE_JACOBIAN__UMFPACK +#define HAVE_ROW_SPARSE_JACOBIAN +#endif + #endif /* AHFINDERDIRECT__CONFIG_H */ |