31
32
33
34#include "implicit_f.inc"
35
36
37
38 INTEGER N
40 . mt(n,*), st(*)
41
42
43
44 INTEGER I,II,I1,M,K,INFO,IPIV(N)
45#ifndef MYREAL8
46 DOUBLE PRECISION :: MTDB(N,N)
47 DOUBLE PRECISION :: STDB(N)
48#endif
49
50#ifdef WITHOUT_LINALG
51 WRITE(6,*) "Error: Lapack needed for Flexible bodies"
52 CALL flush(6)
54#else
55#ifdef MYREAL8
56 CALL dgetrf(n,n,mt,n,ipiv,info)
57 CALL dgetrs(
'N',n,1,mt,n,ipiv,st,n,info)
58#else
59 mtdb(1:n,1:n) = mt(1:n,1:n)
60 CALL dgetrf(n,n,mtdb,n,ipiv,info)
61 stdb(1:n) = st(1:n)
62 CALL dgetrs(
'N',n,1,mtdb,n,ipiv,stdb,n,info)
63 st(1:n) = stdb(1:n)
64#endif
65#endif
66
subroutine dgetrs(trans, n, nrhs, a, lda, ipiv, b, ldb, info)
DGETRS
subroutine dgetrf(m, n, a, lda, ipiv, info)
DGETRF