1 files changed, 26 insertions, 0 deletions
diff --git a/src/tensor4.F90 b/src/tensor4.F90
new file mode 100644
index 0000000..ad03df7
--- /dev/null
+++ b/src/tensor4.F90
@@ -0,0 +1,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