1 SUBROUTINE SGBTRF
( M
, KL
, KU
, AB
, LDAB
)
3 * Slightly modified version of
20000601-1.f that still ICES with
4 * CVS
20010118 g77 on mips
-sgi
-irix6
.5
/-mabi
=64.
6 * Originally derived from LAPACK
3.0 test suite failure
.
8 * David Billinghurst
, (David
.Billinghurst@riotinto
.com
.au
)
11 INTEGER KL
, KU
, LDAB
, M
14 INTEGER J
, JB
, JJ
, JP
, KV
, KM
, F
15 REAL WORK13
(65,64), WORK31
(65,64)
21 JP
= F
( KM
+1, AB
( KV
+1, JJ
) )
22 CALL SSWAP
( JB
, AB
( KV
+1+JJ
-J
, J
), LDAB
-1,
23 $ AB
( KV
+JP
+JJ
-J
, J
), LDAB
-1 )