version 1.7, 2010/12/21 13:53:40
|
version 1.8, 2011/07/22 07:38:12
|
Line 1
|
Line 1
|
SUBROUTINE DTGEVC( SIDE, HOWMNY, SELECT, N, S, LDS, P, LDP, VL, |
SUBROUTINE DTGEVC( SIDE, HOWMNY, SELECT, N, S, LDS, P, LDP, VL, |
$ LDVL, VR, LDVR, MM, M, WORK, INFO ) |
$ LDVL, VR, LDVR, MM, M, WORK, INFO ) |
* |
* |
* -- LAPACK routine (version 3.2) -- |
* -- LAPACK 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 .. |
CHARACTER HOWMNY, SIDE |
CHARACTER HOWMNY, SIDE |
Line 631
|
Line 631
|
* to underflow. (E.g., less than SMALL.) |
* 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 |
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 |
DO 110 JA = 1, NA |
SUMS( JA, JW ) = ZERO |
SUMS( JA, JW ) = ZERO |
SUMP( JA, JW ) = ZERO |
SUMP( JA, JW ) = ZERO |
Line 675 CIBM PREFER SCALAR
|
Line 647 CIBM PREFER SCALAR
|
110 CONTINUE |
110 CONTINUE |
120 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 |
DO 130 JA = 1, NA |
IF( ILCPLX ) THEN |
IF( ILCPLX ) THEN |
SUM( JA, 1 ) = -ACOEF*SUMS( JA, 1 ) + |
SUM( JA, 1 ) = -ACOEF*SUMS( JA, 1 ) + |