diff options
Diffstat (limited to 'src/lapack.F90')
| -rw-r--r-- | src/lapack.F90 | 51 |
1 files changed, 51 insertions, 0 deletions
diff --git a/src/lapack.F90 b/src/lapack.F90 new file mode 100644 index 0000000..5bcd70b --- /dev/null +++ b/src/lapack.F90 @@ -0,0 +1,51 @@ +! $Header$ + +#include "cctk.h" + +module lapack + implicit none + public + + interface geev + subroutine sgeev (jobvl, jobvr, n, a, lda, wr, wi, vl, ldvl, vr, & + & ldvr, work, lwork, info) + character jobvl, jobvr + integer info, lda, ldvl, ldvr, lwork, n + real a(lda,n), vl(ldvl,n), vr(ldvr,n), & + & wi(n), work(lwork), wr(n) + end subroutine sgeev + subroutine dgeev (jobvl, jobvr, n, a, lda, wr, wi, vl, ldvl, vr, & + & ldvr, work, lwork, info) + character jobvl, jobvr + integer info, lda, ldvl, ldvr, lwork, n + double precision a(lda,n), vl(ldvl,n), vr(ldvr,n), & + & wi(n), work(lwork), wr(n) + end subroutine dgeev + end interface + + interface gesv + subroutine sgesv (n, nrhs, a, lda, ipiv, b, ldb, info) + implicit none + integer n + integer nrhs + integer lda + real a(lda,n) + integer ipiv(n) + integer ldb + real b(ldb,nrhs) + integer info + end subroutine sgesv + subroutine dgesv (n, nrhs, a, lda, ipiv, b, ldb, info) + implicit none + integer n + integer nrhs + integer lda + double precision a(lda,n) + integer ipiv(n) + integer ldb + double precision b(ldb,nrhs) + integer info + end subroutine dgesv + end interface + +end module lapack |