version 1.7, 2010/12/21 13:53:28
|
version 1.8, 2011/07/22 07:38:06
|
Line 2
|
Line 2
|
$ RELTOL, PIVMIN, D, E, E2, NVAL, AB, C, MOUT, |
$ RELTOL, PIVMIN, D, E, E2, NVAL, AB, C, MOUT, |
$ NAB, WORK, IWORK, INFO ) |
$ NAB, WORK, IWORK, INFO ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.2) -- |
* -- LAPACK auxiliary routine (version 3.3.1) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* November 2006 |
* -- April 2011 -- |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
INTEGER IJOB, INFO, MINP, MMAX, MOUT, N, NBMIN, NITMAX |
INTEGER IJOB, INFO, MINP, MMAX, MOUT, N, NBMIN, NITMAX |
Line 251
|
Line 251
|
* Compute the number of eigenvalues in the initial intervals. |
* Compute the number of eigenvalues in the initial intervals. |
* |
* |
MOUT = 0 |
MOUT = 0 |
*DIR$ NOVECTOR |
|
DO 30 JI = 1, MINP |
DO 30 JI = 1, MINP |
DO 20 JP = 1, 2 |
DO 20 JP = 1, 2 |
TMP1 = D( 1 ) - AB( JI, JP ) |
TMP1 = D( 1 ) - AB( JI, JP ) |
Line 407
|
Line 406
|
TMP2 = MIN( TMP2, -PIVMIN ) |
TMP2 = MIN( TMP2, -PIVMIN ) |
END IF |
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 |
DO 90 J = 2, N |
TMP2 = D( J ) - E2( J-1 ) / TMP2 - TMP1 |
TMP2 = D( J ) - E2( J-1 ) / TMP2 - TMP1 |
IF( TMP2.LE.PIVMIN ) THEN |
IF( TMP2.LE.PIVMIN ) THEN |
Line 487 CIBM PREFER SCALAR
|
Line 471 CIBM PREFER SCALAR
|
100 CONTINUE |
100 CONTINUE |
KL = KLNEW |
KL = KLNEW |
* |
* |
* End of Serial Version of the loop |
|
* |
|
END IF |
END IF |
* |
* |
* Check for convergence |
* Check for convergence |