version 1.3, 2010/08/06 15:28:41
|
version 1.7, 2010/12/21 13:48:05
|
Line 2
|
Line 2
|
$ SR, SI, H, LDH, ILOZ, IHIZ, Z, LDZ, V, LDV, U, |
$ SR, SI, H, LDH, ILOZ, IHIZ, Z, LDZ, V, LDV, U, |
$ LDU, NV, WV, LDWV, NH, WH, LDWH ) |
$ LDU, NV, WV, LDWV, NH, WH, LDWH ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.2) -- |
* -- LAPACK auxiliary routine (version 3.3.0) -- |
* -- 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 |
* 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 453
|
Line 453
|
* ==== 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 100 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 100 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 ) - REFSUM*V( 2, M22 ) |
H( J, K+1 ) = H( J, K+1 ) - REFSUM |
100 CONTINUE |
H( J, K+2 ) = H( J, K+2 ) - REFSUM*V( 2, M22 ) |
* |
100 CONTINUE |
IF( ACCUM ) THEN |
* |
KMS = K - INCOL |
IF( ACCUM ) THEN |
DO 110 J = MAX( 1, KTOP-INCOL ), KDU |
KMS = K - INCOL |
REFSUM = V( 1, M22 )*( U( J, KMS+1 )+V( 2, M22 )* |
DO 110 J = MAX( 1, KTOP-INCOL ), KDU |
$ U( J, KMS+2 ) ) |
REFSUM = V( 1, M22 )*( U( J, KMS+1 )+ |
U( J, KMS+1 ) = U( J, KMS+1 ) - REFSUM |
$ V( 2, M22 )*U( J, KMS+2 ) ) |
U( J, KMS+2 ) = U( J, KMS+2 ) - REFSUM*V( 2, M22 ) |
U( J, KMS+1 ) = U( J, KMS+1 ) - REFSUM |
|
U( J, KMS+2 ) = U( J, KMS+2 ) - |
|
$ REFSUM*V( 2, M22 ) |
110 CONTINUE |
110 CONTINUE |
ELSE IF( WANTZ ) THEN |
ELSE IF( WANTZ ) THEN |
DO 120 J = ILOZ, IHIZ |
DO 120 J = ILOZ, IHIZ |
REFSUM = V( 1, M22 )*( Z( J, K+1 )+V( 2, M22 )* |
REFSUM = V( 1, M22 )*( Z( J, K+1 )+V( 2, M22 )* |
$ Z( J, K+2 ) ) |
$ Z( J, K+2 ) ) |
Z( J, K+1 ) = Z( J, K+1 ) - REFSUM |
Z( J, K+1 ) = Z( J, K+1 ) - REFSUM |
Z( J, K+2 ) = Z( J, K+2 ) - REFSUM*V( 2, M22 ) |
Z( J, K+2 ) = Z( J, K+2 ) - REFSUM*V( 2, M22 ) |
120 CONTINUE |
120 CONTINUE |
|
END IF |
END IF |
END IF |
END IF |
END IF |
* |
* |