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
|
! $Header$
#include "cctk.h"
module constants
implicit none
private
public delta2, delta3, delta4
public eta4
public epsilon2, epsilon3
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 :: pi = 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068_rk
end module constants
|
