1 files changed, 52 insertions, 0 deletions
diff --git a/src/constants.F90 b/src/constants.F90
new file mode 100644
index 0000000..d9b948d
--- /dev/null
+++ b/src/constants.F90
@@ -0,0 +1,52 @@
+! $Header$
+
+#include "cctk.h"
+#include "cctk_Parameters.h"
+
+module constants
+  implicit none
+  DECLARE_CCTK_PARAMETERS
+  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