rpl/lapack/lapack/zunbdb1.f - diff

Return to zunbdb1.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 203	Line 203
SUBROUTINE ZUNBDB1( M, P, Q, X11, LDX11, X21, LDX21, THETA, PHI,	SUBROUTINE ZUNBDB1( 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 308	Line 308
CALL ZLARF( 'R', M-P-I, Q-I, X21(I,I+1), LDX21, TAUQ1(I),	CALL ZLARF( 'R', M-P-I, Q-I, X21(I,I+1), LDX21, TAUQ1(I),
$ X21(I+1,I+1), LDX21, WORK(ILARF) )	$ X21(I+1,I+1), LDX21, WORK(ILARF) )
CALL ZLACGV( Q-I, X21(I,I+1), LDX21 )	CALL ZLACGV( Q-I, X21(I,I+1), LDX21 )
C = SQRT( DZNRM2( P-I, X11(I+1,I+1), 1, X11(I+1,I+1),	C = SQRT( DZNRM2( P-I, X11(I+1,I+1), 1 )**2
$ 1 )**2 + DZNRM2( M-P-I, X21(I+1,I+1), 1, X21(I+1,I+1),	$ + DZNRM2( M-P-I, X21(I+1,I+1), 1 )**2 )
$ 1 )**2 )
PHI(I) = ATAN2( S, C )	PHI(I) = ATAN2( S, C )
CALL ZUNBDB5( P-I, M-P-I, Q-I-1, X11(I+1,I+1), 1,	CALL ZUNBDB5( P-I, M-P-I, Q-I-1, X11(I+1,I+1), 1,
$ X21(I+1,I+1), 1, X11(I+1,I+2), LDX11,	$ X21(I+1,I+1), 1, X11(I+1,I+2), LDX11,