aboutsummaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
authormiguel <miguel@b7a48df3-cbbf-4440-997f-b4b717c9f7fc>2001-10-19 13:00:20 +0000
committermiguel <miguel@b7a48df3-cbbf-4440-997f-b4b717c9f7fc>2001-10-19 13:00:20 +0000
commit4700a41c206721908b5eeca981861d7558bbcb3d (patch)
treea0ac42c3311dc621db016aceace393f11c559203 /src
parentf0fc3d176035ae304ce73c9a53c375f4c4842109 (diff)
Inverse metric is needed for some matter terms (BosonStar for example).
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMConstraints/trunk@50 b7a48df3-cbbf-4440-997f-b4b717c9f7fc
Diffstat (limited to 'src')
-rw-r--r--src/ADMConstraints.F12
1 files changed, 12 insertions, 0 deletions
diff --git a/src/ADMConstraints.F b/src/ADMConstraints.F
index edcb597..2fe6c23 100644
--- a/src/ADMConstraints.F
+++ b/src/ADMConstraints.F
@@ -47,6 +47,7 @@ c Various real variables.
CCTK_REAL :: dx,dy,dz
CCTK_REAL :: m_rho,m_sx,m_sy,m_sz
CCTK_REAL :: pi,ialp,ialp2
+ CCTK_REAL :: det,uxx,uyy,uzz,uxy,uxz,uyz
c Temporaries for the Stress-Energy tensor.
@@ -98,6 +99,17 @@ c Calculate constraints.
c Calculate the stress energy tensor at this point
c ------------------------------------------------
+c This may be needed for CalcTmunu
+
+#include "CactusEinstein/Einstein/src/macro/DETG_guts.h"
+
+ det = DETG_DETCG
+
+#include "CactusEinstein/Einstein/src/macro/UPPERMET_guts.h"
+
+ uxx = UPPERMET_UXX; uxy = UPPERMET_UXY; uxz = UPPERMET_UXZ
+ uyy = UPPERMET_UYY; uyz = UPPERMET_UYZ; uzz = UPPERMET_UZZ
+
c Initialize stress-energy tensor components.
Ttt = 0.0D0