From 1839a7bfae273274a5738b7debdb5d2b80207b70 Mon Sep 17 00:00:00 2001 From: schnetter Date: Wed, 7 Apr 2010 22:59:31 +0000 Subject: Remove simple routines git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinUtils/TGRtensor/trunk@50 b716e942-a2de-43ad-8f52-f3dfe468e4e7 --- src/adm_metric.F90 | 100 ----------------------------------------------------- 1 file changed, 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) -- cgit v1.2.3