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
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
|
/* config.h -- compile-time configuration for AHFinderDirect */
/* $Header$ */
#ifndef AHFINDERDIRECT__CONFIG_H
#define AHFINDERDIRECT__CONFIG_H
/*
* 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
*/
/* CCTK_VWarn() "warning level" for fatal errors (terminate execution) */
/* note low-level software uses error_exit() instead of CCTK_VWarn(), */
/* with exit codes defined in "stdc.h" */
#define FATAL_ERROR (-1)
/******************************************************************************/
/*
* ***** definitions for integer and floating point datatypes *****
*/
/*
* FIXME: this assumes fp == C 'double'
* CCTK_REAL_PRECISION_{4,8,16} are helpful, but not quite enough)
*/
#include <cctk.h>
#undef FP_IS_FLOAT
#undef FP_IS_DOUBLE
#ifdef CCTK_REAL_PRECISION_4
#define FP_IS_FLOAT
#endif
#ifdef CCTK_REAL_PRECISION_8
#define FP_IS_DOUBLE
#endif
#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
* 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
/* store as row-oriented sparse matrix, solve with UMFPACK */
#define HAVE_ROW_SPARSE_JACOBIAN__UMFPACK
/******************************************************************************/
/******************************************************************************/
/******************************************************************************/
/*
* 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
#ifdef HAVE_ROW_SPARSE_JACOBIAN__UMFPACK
#define HAVE_ROW_SPARSE_JACOBIAN
#endif
#endif /* AHFINDERDIRECT__CONFIG_H */
|
