diff options
| -rw-r--r-- | src/tensor4.F90 | 45 |
1 files changed, 41 insertions, 4 deletions
diff --git a/src/tensor4.F90 b/src/tensor4.F90 index ad03df7..23c2328 100644 --- a/src/tensor4.F90 +++ b/src/tensor4.F90 @@ -3,24 +3,61 @@ #include "cctk.h" module tensor4 + use matinv implicit none private public calc_4trace + public calc_4inv + public calc_4invderiv + contains subroutine calc_4trace (kk, gu, trk) - CCTK_REAL, intent(in) :: kk(4,4) - CCTK_REAL, intent(in) :: gu(4,4) + CCTK_REAL, intent(in) :: kk(0:3,0:3) + CCTK_REAL, intent(in) :: gu(0:3,0:3) CCTK_REAL, intent(out) :: trk integer :: i,j trk = 0 - do i=1,4 - do j=1,4 + do i=0,3 + do j=0,3 trk = trk + gu(i,j) * kk(i,j) end do end do end subroutine calc_4trace + + + subroutine calc_4inv (gg, gu) + CCTK_REAL, intent(in) :: gg(0:3,0:3) + CCTK_REAL, intent(out) :: gu(0:3,0:3) + call calc_inv4 (gg, gu) + end subroutine calc_4inv + + subroutine calc_4pinv (gu, pgg, pgu) + CCTK_REAL, intent(in) :: gu(0:3,0:3), pgg(0:3,0:3) + CCTK_REAL, intent(out) :: pgu(0:3,0:3) + integer :: i,j,k,l + do i=0,3 + do j=0,3 + pgu(i,j) = 0 + do k=0,3 + do l=0,3 + pgu(i,j) = pgu(i,j) - gu(i,k) * gu(j,l) * pgg(k,l) + end do + end do + end do + end do + end subroutine calc_4pinv + + subroutine calc_4invderiv (gu, dgg, dgu) + CCTK_REAL, intent(in) :: gu(0:3,0:3), dgg(0:3,0:3,0:3) + CCTK_REAL, intent(out) :: dgu(0:3,0:3,0:3) + integer :: i + do i=0,3 + call calc_4pinv (gu, dgg(:,:,i), dgu(:,:,i)) + end do + end subroutine calc_4invderiv + end module tensor4 |