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
15 REAL WORK13(65,64), WORK31(65,64)
21 JP = MAX( 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 )