OSDN Git Service

* gcc.c-torture/execute/ieee/rbug.c: Force FP to extended-precision
[pf3gnuchains/gcc-fork.git] / gcc / testsuite / g77.f-torture / compile / 20000601-2.f
1       SUBROUTINE SGBTRF( M, KL, KU, AB, LDAB )
2
3 *  Slightly modified version of 20000601-1.f that still ICES with
4 *  CVS 20010118 g77 on mips-sgi-irix6.5/-mabi=64.
5 *
6 *  Originally derived from LAPACK 3.0 test suite failure.
7 *
8 *  David Billinghurst, (David.Billinghurst@riotinto.com.au)
9 *  18 January 2001
10
11       INTEGER   KL, KU, LDAB, M
12       REAL      AB( LDAB, * )
13
14       INTEGER   J, JB, JJ, JP, KV, KM
15       REAL      WORK13(65,64), WORK31(65,64)
16       KV = KU + KL
17       DO J = 1, M
18          JB = MIN( 1, M-J+1 )
19          DO JJ = J, J + JB - 1
20             KM = MIN( KL, M-JJ )
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 )
24          END DO
25       END DO
26       RETURN
27       END