Add routines for calculating the determinant of a 4d matrix.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinUtils/TGRtensor/trunk@13 b716e942-a2de-43ad-8f52-f3dfe468e4e7

1 files changed, 29 insertions, 0 deletions

diff --git a/src/lapack.F90 b/src/lapack.F90
index cbe41b6..58e64e1 100644
--- a/src/lapack.F90
+++ b/src/lapack.F90
@@ -129,4 +129,33 @@ module lapack
      end subroutine dsysv
   end interface
   
+  
+  
+  ! trf: LU factorisation
+  ! gb (general banded), ge (general), gt (general tridiagonal),
+  ! pb (symmetric positive banded), po (symmetric positive),
+  ! pp (symmetric positive packed), pt (symmetric positive tridiagonal),
+  ! sp (symmetric packed), sy (symmetric)
+  
+  interface sytrf
+     subroutine ssytrf (m, n, a, lda, ipiv, info)
+       implicit none
+       integer m
+       integer n
+       integer lda
+       real    a(lda,n)
+       integer ipiv(n)
+       integer info
+     end subroutine ssytrf
+     subroutine dsytrf (m, n, a, lda, ipiv, info)
+       implicit none
+       integer          m
+       integer          n
+       integer          lda
+       double precision a(lda,n)
+       integer          ipiv(n)
+       integer          info
+     end subroutine dsytrf
+  end interface
+  
 end module lapack