1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
|
! $Header$
#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(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=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
|