diff options
| author | miguel <miguel@b7a48df3-cbbf-4440-997f-b4b717c9f7fc> | 2001-03-20 10:35:30 +0000 |
|---|---|---|
| committer | miguel <miguel@b7a48df3-cbbf-4440-997f-b4b717c9f7fc> | 2001-03-20 10:35:30 +0000 |
| commit | 892f6795a04d899c52b67b2829e41fc69fb2d9f9 (patch) | |
| tree | cc1c77ddbcea723abdc9970cd13c5655cda5bfc0 /src | |
| parent | 2f5872adf80f6f6eb56d34ab55c16c02659d1162 (diff) | |
Cleaning up.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMConstraints/trunk@45 b7a48df3-cbbf-4440-997f-b4b717c9f7fc
Diffstat (limited to 'src')
| -rw-r--r-- | src/ADMConstraints.F | 13 |
1 files changed, 5 insertions, 8 deletions
diff --git a/src/ADMConstraints.F b/src/ADMConstraints.F index c97828d..453a43c 100644 --- a/src/ADMConstraints.F +++ b/src/ADMConstraints.F @@ -113,12 +113,12 @@ c Calculate constraints. c Calculate the stress energy tensor at this point c ------------------------------------------------ -c Fill with zeroes. +c Initialize stress-energy tensor components. Ttt = 0.0D0; Ttx = 0.0D0; Tty = 0.0D0; Ttz = 0.0D0 - Txx = 0.0D0; Txy = 0.0D0; Txz = 0.0D0 - Tyy = 0.0D0; Tyz = 0.0D0; Tzz = 0.0D0 + Txx = 0.0D0; Tyy = 0.0D0; Tzz = 0.0D0 + Txy = 0.0D0; Txz = 0.0D0; Tyz = 0.0D0 c Inverse metric may be needed for CalcTmunu. @@ -127,8 +127,8 @@ c Inverse metric may be needed for CalcTmunu. det = DETG_DETCG - uxx = UPPERMET_UXX; uxy = UPPERMET_UXY; uxz = UPPERMET_UXZ - uyy = UPPERMET_UYY; uyz = UPPERMET_UYZ; uzz = UPPERMET_UZZ + uxx = UPPERMET_UXX; uyy = UPPERMET_UYY; uzz = UPPERMET_UZZ + uxy = UPPERMET_UXY; uxz = UPPERMET_UXZ; uyz = UPPERMET_UYZ c Call macro for stress energy tensor. @@ -150,7 +150,6 @@ c = (T_00 - 2 beta^i T_{i0} + beta^i beta^j T_{ij})/alpha^2 m_rho = ialp2*Ttt - if (shift_state == SHIFT_ACTIVE) then m_rho = m_rho + ialp2 @@ -164,7 +163,6 @@ c = (T_00 - 2 beta^i T_{i0} + beta^i beta^j T_{ij})/alpha^2 & + betay(i,j,k)**2*Tty & + betaz(i,j,k)**2*Ttz)*2.0D0) - end if ham(i,j,k) = HAMADM_HAMADM - 16.0D0*pi*m_rho @@ -197,7 +195,6 @@ c S_i = - g_{ia} n_b T^{ab} = - g_i^a n^b T_{ab} c = - (T_{i0} - beta^j T_{ij})/alpha m_sx = - ialp*Ttx - m_sy = - ialp*Tty m_sz = - ialp*Ttz |