version 1.1, 2014/01/27 09:24:37
|
version 1.3, 2016/08/27 15:27:15
|
Line 201
|
Line 201
|
SUBROUTINE ZUNBDB2( M, P, Q, X11, LDX11, X21, LDX21, THETA, PHI, |
SUBROUTINE ZUNBDB2( M, P, Q, X11, LDX11, X21, LDX21, THETA, PHI, |
$ TAUP1, TAUP2, TAUQ1, WORK, LWORK, INFO ) |
$ TAUP1, TAUP2, TAUQ1, WORK, LWORK, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.5.0) -- |
* -- LAPACK computational routine (version 3.6.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..-- |
* July 2012 |
* July 2012 |
Line 295
|
Line 295
|
CALL ZLARF( 'R', M-P-I+1, Q-I+1, X11(I,I), LDX11, TAUQ1(I), |
CALL ZLARF( 'R', M-P-I+1, Q-I+1, X11(I,I), LDX11, TAUQ1(I), |
$ X21(I,I), LDX21, WORK(ILARF) ) |
$ X21(I,I), LDX21, WORK(ILARF) ) |
CALL ZLACGV( Q-I+1, X11(I,I), LDX11 ) |
CALL ZLACGV( Q-I+1, X11(I,I), LDX11 ) |
S = SQRT( DZNRM2( P-I, X11(I+1,I), 1, X11(I+1,I), |
S = SQRT( DZNRM2( P-I, X11(I+1,I), 1 )**2 |
$ 1 )**2 + DZNRM2( M-P-I+1, X21(I,I), 1, X21(I,I), 1 )**2 ) |
$ + DZNRM2( M-P-I+1, X21(I,I), 1 )**2 ) |
THETA(I) = ATAN2( S, C ) |
THETA(I) = ATAN2( S, C ) |
* |
* |
CALL ZUNBDB5( P-I, M-P-I+1, Q-I, X11(I+1,I), 1, X21(I,I), 1, |
CALL ZUNBDB5( P-I, M-P-I+1, Q-I, X11(I+1,I), 1, X21(I,I), 1, |