version 1.3, 2010/08/06 15:28:58
|
version 1.7, 2010/12/21 13:48:06
|
Line 2
|
Line 2
|
$ H, LDH, ILOZ, IHIZ, Z, LDZ, V, LDV, U, LDU, NV, |
$ H, LDH, ILOZ, IHIZ, Z, LDZ, V, LDV, U, LDU, NV, |
$ WV, LDWV, NH, WH, LDWH ) |
$ WV, LDWV, NH, WH, LDWH ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.2) -- |
* -- LAPACK auxiliary routine (version 3.3.0) -- |
* Univ. of Tennessee, Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd.. |
* Univ. of Tennessee, Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd.. |
* November 2006 |
* November 2010 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
INTEGER IHIZ, ILOZ, KACC22, KBOT, KTOP, LDH, LDU, LDV, |
INTEGER IHIZ, ILOZ, KACC22, KBOT, KTOP, LDH, LDU, LDV, |
Line 441
|
Line 441
|
* ==== Special case: 2-by-2 reflection (if needed) ==== |
* ==== Special case: 2-by-2 reflection (if needed) ==== |
* |
* |
K = KRCOL + 3*( M22-1 ) |
K = KRCOL + 3*( M22-1 ) |
IF( BMP22 .AND. ( V( 1, M22 ).NE.ZERO ) ) THEN |
IF( BMP22 ) THEN |
DO 90 J = JTOP, MIN( KBOT, K+3 ) |
IF ( V( 1, M22 ).NE.ZERO ) THEN |
REFSUM = V( 1, M22 )*( H( J, K+1 )+V( 2, M22 )* |
DO 90 J = JTOP, MIN( KBOT, K+3 ) |
$ H( J, K+2 ) ) |
REFSUM = V( 1, M22 )*( H( J, K+1 )+V( 2, M22 )* |
H( J, K+1 ) = H( J, K+1 ) - REFSUM |
$ H( J, K+2 ) ) |
H( J, K+2 ) = H( J, K+2 ) - |
H( J, K+1 ) = H( J, K+1 ) - REFSUM |
$ REFSUM*DCONJG( V( 2, M22 ) ) |
H( J, K+2 ) = H( J, K+2 ) - |
90 CONTINUE |
|
* |
|
IF( ACCUM ) THEN |
|
KMS = K - INCOL |
|
DO 100 J = MAX( 1, KTOP-INCOL ), KDU |
|
REFSUM = V( 1, M22 )*( U( J, KMS+1 )+V( 2, M22 )* |
|
$ U( J, KMS+2 ) ) |
|
U( J, KMS+1 ) = U( J, KMS+1 ) - REFSUM |
|
U( J, KMS+2 ) = U( J, KMS+2 ) - |
|
$ REFSUM*DCONJG( V( 2, M22 ) ) |
|
100 CONTINUE |
|
ELSE IF( WANTZ ) THEN |
|
DO 110 J = ILOZ, IHIZ |
|
REFSUM = V( 1, M22 )*( Z( J, K+1 )+V( 2, M22 )* |
|
$ Z( J, K+2 ) ) |
|
Z( J, K+1 ) = Z( J, K+1 ) - REFSUM |
|
Z( J, K+2 ) = Z( J, K+2 ) - |
|
$ REFSUM*DCONJG( V( 2, M22 ) ) |
$ REFSUM*DCONJG( V( 2, M22 ) ) |
110 CONTINUE |
90 CONTINUE |
|
* |
|
IF( ACCUM ) THEN |
|
KMS = K - INCOL |
|
DO 100 J = MAX( 1, KTOP-INCOL ), KDU |
|
REFSUM = V( 1, M22 )*( U( J, KMS+1 )+ |
|
$ V( 2, M22 )*U( J, KMS+2 ) ) |
|
U( J, KMS+1 ) = U( J, KMS+1 ) - REFSUM |
|
U( J, KMS+2 ) = U( J, KMS+2 ) - |
|
$ REFSUM*DCONJG( V( 2, M22 ) ) |
|
100 CONTINUE |
|
ELSE IF( WANTZ ) THEN |
|
DO 110 J = ILOZ, IHIZ |
|
REFSUM = V( 1, M22 )*( Z( J, K+1 )+V( 2, M22 )* |
|
$ Z( J, K+2 ) ) |
|
Z( J, K+1 ) = Z( J, K+1 ) - REFSUM |
|
Z( J, K+2 ) = Z( J, K+2 ) - |
|
$ REFSUM*DCONJG( V( 2, M22 ) ) |
|
110 CONTINUE |
|
END IF |
END IF |
END IF |
END IF |
END IF |
* |
* |