2 subroutine sgbcon
(norm
,n
,kl
,ku
,ab
,ldab
,ipiv
,anorm
,rcond
,work
,iwork
5 C -- LAPACK routine (version 3.0) --
6 C Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
7 C Courant Institute, Argonne National Lab, and Rice University
10 C .. Scalar Arguments ..
12 integer info
,kl
,ku
,ldab
,n
15 C .. Array Arguments ..
16 integer ipiv
(n
),iwork
(n
)
17 real ab
(ldab
,n
),work
(n
)
22 C demonstrate g77 bug at -O -funroll-loops
23 C =====================================================================
27 parameter(one
= 1.0e+0,zero
= 0.0e+0)
32 integer ix
,j
,jp
,kase
,kase1
,kd
,lm
33 real ainvnm
,scale
,smlnum
,t
35 C .. External Functions ..
39 externallsame
,isamax
,sdot
,slamch
41 C .. External Subroutines ..
42 externalsaxpy
,slacon
,slatbs
,srscl
,xerbla
44 C .. Executable Statements ..
49 C the following min() intrinsic provokes this bug
54 C but only when combined with this if block
58 C and this subroutine call
59 call saxpy
(lm
,-t
,ab
(kd
+1,j
),1,work
(j
+1),1)