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
|
#define Power(x, y) (pow(x,y))
#define Sqrt(x) (sqrt(x))
#ifdef KRANC_C
# define Abs(x) (fabs(x))
# define Min(x, y) (fmin(x,y))
# define Min3(x, y, z) (fmin(fmin((x), (y)), (z)))
# define Max(x, y) (fmax(x,y))
# define IfThen(x,y,z) ((x) ? (y) : (z))
#else
# define Abs(x) (abs(x))
# define IntAbs(x) (abs(x))
# define Min(x, y) (min(x,y))
# define Max(x, y) (max(x,y))
# define IfThen(x,y,z) ((x)*(y) + (1-(x))*(z))
#endif
#define Exp(x) (exp(x))
#define Log(x) (log(x))
#define Sin(x) (sin(x))
#define Cos(x) (cos(x))
#define Tan(x) (tan(x))
#define Sec(x) (1.0/cos(x))
#define Csc(x) (1.0/sin(x))
#define Cot(x) (1.0/tan(x))
#define ArcSin(x) (asin(x))
#define ArcCos(x) (acos(x))
#define ArcTan(x) (atan(x))
#define ArcSec(x) (cos(1.0/(x)))
#define ArcCsc(x) (sin(1.0/(x)))
#define ArcCot(x) (tan(1.0/(x)))
#define Sinh(x) (sinh(x))
#define Cosh(x) (cosh(x))
#define Tanh(x) (tanh(x))
#define Sech(x) (1.0/cosh(x))
#define Csch(x) (1.0/sinh(x))
#define Coth(x) (1.0/tanh(x))
#define Sign(x) (sgn(x))
#ifdef KRANC_C
# define ToReal(x) ((CCTK_REAL)(x))
#else
# define ToReal(x) (real((x),kind(khalf)))
#endif
#define MinMod(x, y) ((x) * (y) < 0 ? 0 : (fabs((x)) < fabs((y)) ? (x) : (y)))
#define VanLeer(x, y) ((x) * (y) < 0 ? 0 : (Min3(2*fabs(x),2*fabs(y),0.5*(fabs(x)+fabs(y)))*Sign((x)+(y))))
#ifdef KRANC_C
# define E M_E
# define Pi M_PI
#else
# define E 2.71828182845904523536029d0
# define Pi 3.14159265358979323846264d0
#endif
#define StepFunction(x) ((x)>0)
|