1 SUBROUTINE sgefa1 (A, LDA, N, IPVT, INFO)
2 USE liprec,
ONLY: wp => sp
8 INTEGER,
DIMENSION(N) :: IPVT
9 REAL(WP),
DIMENSION(LDA,N) :: A
11 CALL sgefa (a, lda, n, ipvt, info)
15 SUBROUTINE dgefa1 (A, LDA, N, IPVT, INFO)
16 USE liprec,
ONLY: wp => dp
22 INTEGER,
DIMENSION(N) :: IPVT
23 REAL(WP),
DIMENSION(LDA,N) :: A
25 CALL dgefa (a, lda, n, ipvt, info)
30 SUBROUTINE sgefa (A, LDA, N, IPVT, INFO)
31 USE liprec,
ONLY: wp => sp
37 INTEGER,
DIMENSION(N) :: IPVT
38 REAL(WP),
DIMENSION(LDA,N) :: A
42 REAL(WP),
PARAMETER :: ZERO = 0, one = 1
43 INTEGER :: J, K, KP1, L, NM1
44 INTEGER,
DIMENSION(1) :: ISAMAX
99 isamax = maxloc(abs(a(k:n,k)))
100 l = isamax(1) + k - 1
107 IF (a(l,k) .ne. zero)
THEN
119 element = -one/a(k,k)
121 a(k+1:n,k) = element*a(k+1:n,k)
132 a(k+1:n,j) = a(k+1:n,j) + element*a(k+1:n,k)
141 IF (a(n,n) .eq. zero) info = n
145 SUBROUTINE dgefa (A, LDA, N, IPVT, INFO)
146 USE liprec,
ONLY: wp => dp
152 INTEGER,
DIMENSION(N) :: IPVT
153 REAL(WP),
DIMENSION(LDA,N) :: A
157 REAL(WP),
PARAMETER :: ZERO = 0, one = 1
158 INTEGER :: J, K, KP1, L, NM1
159 INTEGER,
DIMENSION(1) :: ISAMAX
214 isamax = maxloc(abs(a(k:n,k)))
215 l = isamax(1) + k - 1
222 IF (a(l,k) .ne. zero)
THEN
234 element = -one/a(k,k)
236 a(k+1:n,k) = element*a(k+1:n,k)
247 a(k+1:n,j) = a(k+1:n,j) + element*a(k+1:n,k)
256 IF (a(n,n) .eq. zero) info = n