From 676fce297a14688119c36fc23d832e8725a33cd6 Mon Sep 17 00:00:00 2001 From: schnetter Date: Sun, 23 Oct 2005 17:26:46 +0000 Subject: Correct sign of epsilon3. Add constants delta1 and epsilon1. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinUtils/TGRtensor/trunk@36 b716e942-a2de-43ad-8f52-f3dfe468e4e7 --- src/constants.F90 | 17 ++++++++++++----- 1 file 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, & -- cgit v1.2.3