2 c { dg-options "-std=legacy" }
5 * Derived from LAPACK
3.0 routine CHGEQZ
6 * Fails on i686
-pc
-cygwin with gcc
-2.97 snapshots at
-O2 and higher
9 * David Billinghurst
, (David
.Billinghurst@riotinto
.com
)
11 * Rewritten by Toon Moene
(toon@moene
.indiv
.nluug
.nl
)
16 A
(4,3) = (0.05,0.2)/3.0E-7
18 A
(5,4) = (-2.0E-07,2.0E-07)
21 SUBROUTINE CHGEQZ
( N
, A
)
23 ABS1
( X
) = ABS
( REAL( X
) ) + ABS
( AIMAG
( X
) )
27 TEMP2
= ABS1
( A
( J
+1, J
) )
28 TEMPR
= MAX
( TEMP
, TEMP2
)
29 IF( TEMPR
.LT
. 1.0 .AND
. TEMPR
.NE
. 0.0 ) THEN
33 IF ( ABS1
(A
(J
,J
-1))*TEMP2
.LE
. TEMP
) GO TO 90
35 c Should not reach here, but need a statement
37 90 IF ( I
.NE
. 4 ) THEN
38 PRINT*
,'I =', I
, ' but should be 4'