TMP2 = MIN( TMP2, -PIVMIN )
END IF
*
-* A series of compiler directives to defeat vectorization
-* for the next loop
-*
-*$PL$ CMCHAR=' '
-CDIR$ NEXTSCALAR
-C$DIR SCALAR
-CDIR$ NEXT SCALAR
-CVD$L NOVECTOR
-CDEC$ NOVECTOR
-CVD$ NOVECTOR
-*VDIR NOVECTOR
-*VOCL LOOP,SCALAR
-CIBM PREFER SCALAR
-*$PL$ CMCHAR='*'
-*
DO 90 J = 2, N
TMP2 = D( J ) - E2( J-1 ) / TMP2 - TMP1
IF( TMP2.LE.PIVMIN ) THEN
100 CONTINUE
KL = KLNEW
*
-* End of Serial Version of the loop
-*
END IF
*
* Check for convergence
WORK( N ) = ZERO
PIVMIN = ONE
*
-*DIR$ NOVECTOR
DO 10 J = 2, N
TMP1 = E( J-1 )**2
IF( ABS( D( J )*D( J-1 ) )*ULP**2+SAFEMN.GT.TMP1 ) THEN
* to underflow. (E.g., less than SMALL.)
*
*
-* A series of compiler directives to defeat vectorization
-* for the next loop
-*
-*$PL$ CMCHAR=' '
-CDIR$ NEXTSCALAR
-C$DIR SCALAR
-CDIR$ NEXT SCALAR
-CVD$L NOVECTOR
-CDEC$ NOVECTOR
-CVD$ NOVECTOR
-*VDIR NOVECTOR
-*VOCL LOOP,SCALAR
-CIBM PREFER SCALAR
-*$PL$ CMCHAR='*'
-*
DO 120 JW = 1, NW
-*
-*$PL$ CMCHAR=' '
-CDIR$ NEXTSCALAR
-C$DIR SCALAR
-CDIR$ NEXT SCALAR
-CVD$L NOVECTOR
-CDEC$ NOVECTOR
-CVD$ NOVECTOR
-*VDIR NOVECTOR
-*VOCL LOOP,SCALAR
-CIBM PREFER SCALAR
-*$PL$ CMCHAR='*'
-*
DO 110 JA = 1, NA
SUMS( JA, JW ) = ZERO
SUMP( JA, JW ) = ZERO
110 CONTINUE
120 CONTINUE
*
-*$PL$ CMCHAR=' '
-CDIR$ NEXTSCALAR
-C$DIR SCALAR
-CDIR$ NEXT SCALAR
-CVD$L NOVECTOR
-CDEC$ NOVECTOR
-CVD$ NOVECTOR
-*VDIR NOVECTOR
-*VOCL LOOP,SCALAR
-CIBM PREFER SCALAR
-*$PL$ CMCHAR='*'
-*
DO 130 JA = 1, NA
IF( ILCPLX ) THEN
SUM( JA, 1 ) = -ACOEF*SUMS( JA, 1 ) +
* Compute the number of eigenvalues in the initial intervals.
*
MOUT = 0
-CDIR$ NOVECTOR
DO 30 JI = 1, MINP
DO 20 JP = 1, 2
TMP1 = D( 1 ) - AB( JI, JP )
TMP2 = MIN( TMP2, -PIVMIN )
END IF
*
-* A series of compiler directives to defeat vectorization
-* for the next loop
-*
-*$PL$ CMCHAR=' '
-CDIR$ NEXTSCALAR
-C$DIR SCALAR
-CDIR$ NEXT SCALAR
-CVD$L NOVECTOR
-CDEC$ NOVECTOR
-CVD$ NOVECTOR
-*VDIR NOVECTOR
-*VOCL LOOP,SCALAR
-CIBM PREFER SCALAR
-*$PL$ CMCHAR='*'
-*
DO 90 J = 2, N
TMP2 = D( J ) - E2( J-1 ) / TMP2 - TMP1
IF( TMP2.LE.PIVMIN ) THEN
100 CONTINUE
KL = KLNEW
*
-* End of Serial Version of the loop
-*
END IF
*
* Check for convergence
WORK( N ) = ZERO
PIVMIN = ONE
*
-CDIR$ NOVECTOR
DO 10 J = 2, N
TMP1 = E( J-1 )**2
IF( ABS( D( J )*D( J-1 ) )*ULP**2+SAFEMN.GT.TMP1 ) THEN
* to underflow. (E.g., less than SMALL.)
*
*
-* A series of compiler directives to defeat vectorization
-* for the next loop
-*
-*$PL$ CMCHAR=' '
-CDIR$ NEXTSCALAR
-C$DIR SCALAR
-CDIR$ NEXT SCALAR
-CVD$L NOVECTOR
-CDEC$ NOVECTOR
-CVD$ NOVECTOR
-*VDIR NOVECTOR
-*VOCL LOOP,SCALAR
-CIBM PREFER SCALAR
-*$PL$ CMCHAR='*'
-*
DO 120 JW = 1, NW
-*
-*$PL$ CMCHAR=' '
-CDIR$ NEXTSCALAR
-C$DIR SCALAR
-CDIR$ NEXT SCALAR
-CVD$L NOVECTOR
-CDEC$ NOVECTOR
-CVD$ NOVECTOR
-*VDIR NOVECTOR
-*VOCL LOOP,SCALAR
-CIBM PREFER SCALAR
-*$PL$ CMCHAR='*'
-*
DO 110 JA = 1, NA
SUMS( JA, JW ) = ZERO
SUMP( JA, JW ) = ZERO
110 CONTINUE
120 CONTINUE
*
-*$PL$ CMCHAR=' '
-CDIR$ NEXTSCALAR
-C$DIR SCALAR
-CDIR$ NEXT SCALAR
-CVD$L NOVECTOR
-CDEC$ NOVECTOR
-CVD$ NOVECTOR
-*VDIR NOVECTOR
-*VOCL LOOP,SCALAR
-CIBM PREFER SCALAR
-*$PL$ CMCHAR='*'
-*
DO 130 JA = 1, NA
IF( ILCPLX ) THEN
SUM( JA, 1 ) = -ACOEF*SUMS( JA, 1 ) +