2 * Derived from LAPACK
3.0 routine CHGEQZ
3 * Fails on i686
-pc
-cygwin with gcc
-2.97 snapshots at
-O2 and higher
6 * David Billinghurst
, (David
.Billinghurst@riotinto
.com
)
8 * Rewritten by Toon Moene
(toon@moene
.indiv
.nluug
.nl
)
13 A
(4,3) = (0.05,0.2)/3.0E-7
15 A
(5,4) = (-2.0E-07,2.0E-07)
18 SUBROUTINE CHGEQZ
( N
, A
)
20 ABS1
( X
) = ABS
( REAL( X
) ) + ABS
( AIMAG
( X
) )
24 TEMP2
= ABS1
( A
( J
+1, J
) )
25 TEMPR
= MAX
( TEMP
, TEMP2
)
26 IF( TEMPR
.LT
. 1.0 .AND
. TEMPR
.NE
. 0.0 ) THEN
30 IF ( ABS1
(A
(J
,J
-1))*TEMP2
.LE
. TEMP
) GO TO 90
32 c Should not reach here, but need a statement
34 90 IF ( I
.NE
. 4 ) THEN
35 PRINT*
,'I =', I
, ' but should be 4'