rpl/lapack/lapack/dorbdb1.f - diff

Return to dorbdb1.f CVS log

Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/dorbdb1.f between versions 1.1 and 1.9

version 1.1, 2014/01/27 09:24:34	version 1.9, 2023/08/07 08:39:01
Line 2	Line 2
*	*
* =========== DOCUMENTATION ===========	* =========== DOCUMENTATION ===========
*	*
* Online html documentation available at	* Online html documentation available at
* http://www.netlib.org/lapack/explore-html/	* http://www.netlib.org/lapack/explore-html/
*	*
*> \htmlonly	*> \htmlonly
*> Download DORBDB1 + dependencies	*> Download DORBDB1 + dependencies
Line 20	Line 20
*	*
* SUBROUTINE DORBDB1( M, P, Q, X11, LDX11, X21, LDX21, THETA, PHI,	* SUBROUTINE DORBDB1( M, P, Q, X11, LDX11, X21, LDX21, THETA, PHI,
* TAUP1, TAUP2, TAUQ1, WORK, LWORK, INFO )	* TAUP1, TAUP2, TAUQ1, WORK, LWORK, INFO )
*	*
* .. Scalar Arguments ..	* .. Scalar Arguments ..
* INTEGER INFO, LWORK, M, P, Q, LDX11, LDX21	* INTEGER INFO, LWORK, M, P, Q, LDX11, LDX21
* ..	* ..
Line 29	Line 29
* DOUBLE PRECISION TAUP1(), TAUP2(), TAUQ1(), WORK(),	* DOUBLE PRECISION TAUP1(), TAUP2(), TAUQ1(), WORK(),
* $ X11(LDX11,), X21(LDX21,)	* $ X11(LDX11,), X21(LDX21,)
* ..	* ..
*	*
*	*
*> \par Purpose:	*> \par Purpose:
*> =============	* =============
*>	*>
*>\verbatim	*>\verbatim
*>	*>
Line 151	Line 151
*> \verbatim	*> \verbatim
*> LWORK is INTEGER	*> LWORK is INTEGER
*> The dimension of the array WORK. LWORK >= M-Q.	*> The dimension of the array WORK. LWORK >= M-Q.
*>	*>
*> If LWORK = -1, then a workspace query is assumed; the routine	*> If LWORK = -1, then a workspace query is assumed; the routine
*> only calculates the optimal size of the WORK array, returns	*> only calculates the optimal size of the WORK array, returns
*> this value as the first entry of the WORK array, and no error	*> this value as the first entry of the WORK array, and no error
Line 169	Line 169
* Authors:	* Authors:
* ========	* ========
*	*
*> \author Univ. of Tennessee	*> \author Univ. of Tennessee
*> \author Univ. of California Berkeley	*> \author Univ. of California Berkeley
*> \author Univ. of Colorado Denver	*> \author Univ. of Colorado Denver
*> \author NAG Ltd.	*> \author NAG Ltd.
*
*> \date July 2012
*	*
*> \ingroup doubleOTHERcomputational	*> \ingroup doubleOTHERcomputational
*	*
Line 203	Line 201
SUBROUTINE DORBDB1( M, P, Q, X11, LDX11, X21, LDX21, THETA, PHI,	SUBROUTINE DORBDB1( 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 --
* -- 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
*	*
* .. Scalar Arguments ..	* .. Scalar Arguments ..
INTEGER INFO, LWORK, M, P, Q, LDX11, LDX21	INTEGER INFO, LWORK, M, P, Q, LDX11, LDX21
Line 304	Line 301
$ X11(I+1,I+1), LDX11, WORK(ILARF) )	$ X11(I+1,I+1), LDX11, WORK(ILARF) )
CALL DLARF( 'R', M-P-I, Q-I, X21(I,I+1), LDX21, TAUQ1(I),	CALL DLARF( '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) )
C = SQRT( DNRM2( P-I, X11(I+1,I+1), 1, X11(I+1,I+1),	C = SQRT( DNRM2( P-I, X11(I+1,I+1), 1 )**2
$ 1 )**2 + DNRM2( M-P-I, X21(I+1,I+1), 1, X21(I+1,I+1),	$ + DNRM2( M-P-I, X21(I+1,I+1), 1 )**2 )
$ 1 )**2 )
PHI(I) = ATAN2( S, C )	PHI(I) = ATAN2( S, C )
CALL DORBDB5( P-I, M-P-I, Q-I-1, X11(I+1,I+1), 1,	CALL DORBDB5( 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,

CVSweb interface <joel.bertrand@systella.fr>

Removed from v.1.1
changed lines
	Added in v.1.9