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