rpl/lapack/lapack/zunbdb2.f - diff

Return to zunbdb2.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.2, 2014/01/27 09:28:44	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,