rpl/lapack/lapack/dorbdb3.f - diff

Return to dorbdb3.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.2, 2014/01/27 09:28:24	version 1.3, 2016/08/27 15:27:10
Line 201	Line 201
SUBROUTINE DORBDB3( M, P, Q, X11, LDX11, X21, LDX21, THETA, PHI,	SUBROUTINE DORBDB3( 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 292	Line 292
$ X11(I,I), LDX11, WORK(ILARF) )	$ X11(I,I), LDX11, WORK(ILARF) )
CALL DLARF( 'R', M-P-I, Q-I+1, X21(I,I), LDX21, TAUQ1(I),	CALL DLARF( 'R', M-P-I, Q-I+1, X21(I,I), LDX21, TAUQ1(I),
$ X21(I+1,I), LDX21, WORK(ILARF) )	$ X21(I+1,I), LDX21, WORK(ILARF) )
C = SQRT( DNRM2( P-I+1, X11(I,I), 1, X11(I,I),	C = SQRT( DNRM2( P-I+1, X11(I,I), 1 )**2
$ 1 )2 + DNRM2( M-P-I, X21(I+1,I), 1, X21(I+1,I), 1 )2 )	$ + DNRM2( M-P-I, X21(I+1,I), 1 )**2 )
THETA(I) = ATAN2( S, C )	THETA(I) = ATAN2( S, C )
*	*
CALL DORBDB5( P-I+1, M-P-I, Q-I, X11(I,I), 1, X21(I+1,I), 1,	CALL DORBDB5( P-I+1, M-P-I, Q-I, X11(I,I), 1, X21(I+1,I), 1,