1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
|
/* config.h -- compile-time configuration for AHFinderDirect */
/* $Header$ */
/*
* This header file is #included by C++ and C files, and also by
* Fortran 77 files which are run through a C-like preprocessor. Alas,
* due to a Cactus bug on some systems, this means that all preprocessor
* commands must start in column 1. :( :( :(
*/
/*
* prerequisites:
* "cctk.h" or "fake_cctk.h"
*/
/******************************************************************************/
/*
* definitions for C/C++ only
*/
#ifdef CCODE
typedef CCTK_REAL fp;
/* Fortran 'integer' type */
typedef CCTK_INT integer;
#endif /* CCODE */
/******************************************************************************/
/*
* definitions for C++, C, and Fortran
*/
/*
* FIXME: this assumes fp == C 'double'
* CCTK_REAL_PRECISION_{4,8,16} are helpful, but not quite enough)
*/
#undef FP_IS_FLOAT
#define FP_IS_DOUBLE
#if defined(FP_IS_FLOAT)
#define FP_SCANF_FORMAT "%f"
#elif defined(FP_IS_DOUBLE)
#define FP_SCANF_FORMAT "%lf"
#else
#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.
*/
#ifdef FINITE_DIFF_ORDER
/* someone (probably a Makefile) has already defined this */
/* ==> leave it alone */
#else
#define FINITE_DIFF_ORDER 4
#endif
/**************************************/
/*
* 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.
*/
/* 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
/**************************************/
/*
* The following #ifdefs and #defines shouldn't need changing; they
* decode the previous set of Jacobian storage format and linear solver
* #defines to determine which storage formats we are using.
*/
#ifdef HAVE_DENSE_JACOBIAN__LAPACK
#define HAVE_DENSE_JACOBIAN
#endif
#ifdef HAVE_ROW_SPARSE_JACOBIAN__ILUCG
#define HAVE_ROW_SPARSE_JACOBIAN
#endif
|