diff options
Diffstat (limited to 'src/constants.F90')
| -rw-r--r-- | src/constants.F90 | 17 |
1 files changed, 12 insertions, 5 deletions
diff --git a/src/constants.F90 b/src/constants.F90 index 5f162dc..54e47ea 100644 --- a/src/constants.F90 +++ b/src/constants.F90 @@ -6,9 +6,9 @@ module constants implicit none private - public delta2, delta3, delta4 + public delta1, delta2, delta3, delta4 public eta4 - public epsilon2, epsilon3, epsilon4 + public epsilon1, epsilon2, epsilon3, epsilon4 public pi @@ -19,6 +19,10 @@ module constants integer, parameter :: rk = kind(zero) + + CCTK_REAL, parameter :: delta1(1) & + = reshape((/ 1 /), (/1/)) + CCTK_REAL, parameter :: delta2(2,2) & = reshape((/ 1,0, 0,1 /), (/2,2/)) @@ -35,13 +39,16 @@ module constants + CCTK_REAL, parameter :: epsilon1(1) & + = reshape ((/ +1 /), (/1/)) + 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/)) + = 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, & |