aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorschnetter <schnetter@b716e942-a2de-43ad-8f52-f3dfe468e4e7>2010-04-07 22:59:31 +0000
committerschnetter <schnetter@b716e942-a2de-43ad-8f52-f3dfe468e4e7>2010-04-07 22:59:31 +0000
commit1839a7bfae273274a5738b7debdb5d2b80207b70 (patch)
tree04d6c670b14f7abbec1226d56b2df2198ffcfc95
parent192e5a703e89bf974eaffc42d44f68658ac4a7ba (diff)
Remove simple routines
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinUtils/TGRtensor/trunk@50 b716e942-a2de-43ad-8f52-f3dfe468e4e7
-rw-r--r--src/adm_metric.F90100
1 files changed, 0 insertions, 100 deletions
diff --git a/src/adm_metric.F90 b/src/adm_metric.F90
index a9ef4a4..1286569 100644
--- a/src/adm_metric.F90
+++ b/src/adm_metric.F90
@@ -19,11 +19,6 @@ module adm_metric
!!$ public calc_extcurvdot
!!$ public calc_3metricdot2
- public calc_3metricdot_simple
- public calc_3metricderivdot_simple
- public calc_extcurvdot_simple
- public calc_3metricdot2_simple
-
!!$ public simplify_t4
contains
@@ -317,64 +312,9 @@ contains
end do
end subroutine calc_3metricderivdot
-!!$ subroutine calc_extcurvdot (gg,gu,ri, kk, t44, kk_dot)
-!!$ CCTK_REAL, intent(in) :: gg(3,3), gu(3,3)
-!!$ CCTK_REAL, intent(in) :: ri(3,3)
-!!$ CCTK_REAL, intent(in) :: kk(3,3)
-!!$ CCTK_REAL, intent(out) :: kk_dot(3,3)
-!!$ integer :: i,j,l,m
-!!$
-!!$ ! d/dt K_ij = R_ij - 2 K_il g^lm K_mj + g^lm K_lm K_ij
-!!$ ! - 8 pi (T_ij - 1/2 g_ij T)
-!!$
-!!$ do i=1,3
-!!$ do j=1,3
-!!$ kk_dot(i,j) = ri(i,j) - 8*pi * tt(i,j)
-!!$ do l=1,3
-!!$ do m=1,3
-!!$ kk_dot(i,j) = kk_dot(i,j) &
-!!$ + gu(l,m) * (-2 * kk(i,l) * kk(m,j) + kk(l,m) * kk(i,j)) &
-!!$ + 4*pi * gg(i,j) * gu(l,m) * tt(l,m)
-!!$ end do
-!!$ end do
-!!$ end do
-!!$ end do
-!!$ end subroutine calc_extcurvdot_simple
-!!$
-!!$ subroutine calc_3metricdot2 (kk_dot, gg_dot2)
-!!$ CCTK_REAL, intent(in) :: kk_dot(3,3)
-!!$ CCTK_REAL, intent(out) :: gg_dot2(3,3)
-!!$ integer :: i,j
-!!$
-!!$ ! d/dt g_ij = -2 alpha K_ij + g_kj beta^k,i + g_ik beta^k,j + beta^k g_ij,k
-!!$
-!!$ ! d^2/dt^2 g_ij = -2 alpha K_ij + g_kj beta^k,i + g_ik beta^k,j + beta^k g_ij,k
-!!$
-!!$ do i=1,3
-!!$ do j=1,3
-!!$ gg_dot2(i,j) = -2 * kk_dot(i,j)
-!!$ end do
-!!$ end do
-!!$ end subroutine calc_3metricdot2
- subroutine calc_extcurv_simple (gg_dot, kk)
- CCTK_REAL, intent(in) :: gg_dot(3,3)
- CCTK_REAL, intent(out) :: kk(3,3)
- integer :: i,j
-
- ! d/dt g_ij = -2 K_ij
-
- do i=1,3
- do j=1,3
- kk(i,j) = - gg_dot(i,j)
- kk(i,j) = kk(i,j) / 2
- end do
- end do
-
- end subroutine calc_extcurv_simple
-
subroutine calc_3metricdot_simple (kk, gg_dot)
CCTK_REAL, intent(in) :: kk(3,3)
CCTK_REAL, intent(out) :: gg_dot(3,3)
@@ -405,46 +345,6 @@ contains
end do
end subroutine calc_3metricderivdot_simple
- subroutine calc_extcurvdot_simple (gg,gu,ri, kk, tt, kk_dot)
- CCTK_REAL, intent(in) :: gg(3,3), gu(3,3)
- CCTK_REAL, intent(in) :: ri(3,3)
- CCTK_REAL, intent(in) :: kk(3,3)
- CCTK_REAL, intent(in) :: tt(3,3)
- CCTK_REAL, intent(out) :: kk_dot(3,3)
- integer :: i,j,l,m
-
- ! d/dt K_ij = R_ij - 2 K_il g^lm K_mj + g^lm K_lm K_ij
- ! - 8 pi (T_ij - 1/2 g_ij T)
-
- do i=1,3
- do j=1,3
- kk_dot(i,j) = ri(i,j) - 8*pi * tt(i,j)
- do l=1,3
- do m=1,3
- kk_dot(i,j) = kk_dot(i,j) &
- + gu(l,m) * (-2 * kk(i,l) * kk(m,j) + kk(l,m) * kk(i,j)) &
- + 4*pi * gg(i,j) * gu(l,m) * tt(l,m)
- end do
- end do
- end do
- end do
- end subroutine calc_extcurvdot_simple
-
- subroutine calc_3metricdot2_simple (kk_dot, gg_dot2)
- CCTK_REAL, intent(in) :: kk_dot(3,3)
- CCTK_REAL, intent(out) :: gg_dot2(3,3)
- integer :: i,j
-
- ! d/dt g_ij = -2 K_ij
- ! d^2/dt^2 g_ij = -2 d/dt K_ij
-
- do i=1,3
- do j=1,3
- gg_dot2(i,j) = -2 * kk_dot(i,j)
- end do
- end do
- end subroutine calc_3metricdot2_simple
-
!!$ subroutine simplify_t4 (t4, alfa, beta, t4s)