blob: ad03df782441c84881659c8814e048cc10a6793f (
plain)
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
|
! $Header$
#include "cctk.h"
module tensor4
implicit none
private
public calc_4trace
contains
subroutine calc_4trace (kk, gu, trk)
CCTK_REAL, intent(in) :: kk(4,4)
CCTK_REAL, intent(in) :: gu(4,4)
CCTK_REAL, intent(out) :: trk
integer :: i,j
trk = 0
do i=1,4
do j=1,4
trk = trk + gu(i,j) * kk(i,j)
end do
end do
end subroutine calc_4trace
end module tensor4
|
