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
|
! $Header$
#include "cctk.h"
module constants
implicit none
private
public delta2, delta3, delta4
public eta4
public epsilon2, epsilon3, epsilon4
public pi
CCTK_REAL, parameter :: zero = 0
integer, parameter :: rk = kind(zero)
CCTK_REAL, parameter :: delta2(2,2) &
= reshape((/ 1,0, 0,1 /), (/2,2/))
CCTK_REAL, parameter :: delta3(3,3) &
= reshape((/ 1,0,0, 0,1,0, 0,0,1 /), (/3,3/))
CCTK_REAL, parameter :: delta4(0:3,0:3) &
= reshape((/ 1,0,0,0, 0,1,0,0, 0,0,1,0, 0,0,0,1 /), (/4,4/))
CCTK_REAL, parameter :: eta4(0:3,0:3) &
= reshape((/ -1,0,0,0, 0,1,0,0, 0,0,1,0, 0,0,0,1 /), (/4,4/))
CCTK_REAL, parameter :: epsilon2(2,2) &
= reshape ((/ 0,+1, -1,0 /), (/2,2/))
CCTK_REAL, parameter :: epsilon3(3,3,3) &
= reshape ((/ 0,0,0, 0,0,+1, 0,-1,0, &
& 0,0,-1, 0,0,0, +1,0,0, &
& 0,+1,0, -1,0,0, 0,0,0 /), (/3,3,3/))
CCTK_REAL, parameter :: epsilon4(0:3,0:3,0:3,0:3) &
= reshape ((/ 0,0,0,0, 0,0,0,0, 0,0,0,0, 0,0,0,0, &
& 0,0,0,0, 0,0,0,0, 0,0,0,+1, 0,0,-1,0, &
& 0,0,0,0, 0,0,0,-1, 0,0,0,0, 0,+1,0,0, &
& 0,0,0,0, 0,0,+1,0, 0,-1,0,0, 0,0,0,0, &
& 0,0,0,0, 0,0,0,0, 0,0,0,-1, 0,0,+1,0, &
& 0,0,0,0, 0,0,0,0, 0,0,0,0, 0,0,0,0, &
& 0,0,0,+1, 0,0,0,0, 0,0,0,0, -1,0,0,0, &
& 0,0,-1,0, 0,0,0,0, +1,0,0,0, 0,0,0,0, &
& 0,0,0,0, 0,0,0,+1, 0,0,0,0, 0,-1,0,0, &
& 0,0,0,-1, 0,0,0,0, 0,0,0,0, +1,0,0,0, &
& 0,0,0,0, 0,0,0,0, 0,0,0,0, 0,0,0,0, &
& 0,+1,0,0, -1,0,0,0, 0,0,0,0, 0,0,0,0, &
& 0,0,0,0, 0,0,-1,0, 0,+1,0,0, 0,0,0,0, &
& 0,0,+1,0, 0,0,0,0, -1,0,0,0, 0,0,0,0, &
& 0,-1,0,0, +1,0,0,0, 0,0,0,0, 0,0,0,0, &
& 0,0,0,0, 0,0,0,0, 0,0,0,0, 0,0,0,0 /), &
& (/4,4,4,4/))
CCTK_REAL, parameter :: pi = 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068_rk
end module constants
|