rpl/lapack/lapack/dlatrz.f - diff

Return to dlatrz.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.2, 2010/04/21 13:45:20	version 1.9, 2011/07/22 07:38:08
Line 1	Line 1
SUBROUTINE DLATRZ( M, N, L, A, LDA, TAU, WORK )	SUBROUTINE DLATRZ( M, N, L, A, LDA, TAU, WORK )
*	*
* -- LAPACK routine (version 3.2) --	* -- LAPACK routine (version 3.3.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..--
* November 2006	* -- April 2011 --
*	*
* .. Scalar Arguments ..	* .. Scalar Arguments ..
INTEGER L, LDA, M, N	INTEGER L, LDA, M, N
Line 64	Line 64
*	*
* where	* where
*	*
* T( k ) = I - tauu( k )u( k )', u( k ) = ( 1 ),	* T( k ) = I - tauu( k )u( k )**T, u( k ) = ( 1 ),
* ( 0 )	* ( 0 )
* ( z( k ) )	* ( z( k ) )
*	*
Line 90	Line 90
INTEGER I	INTEGER I
* ..	* ..
* .. External Subroutines ..	* .. External Subroutines ..
EXTERNAL DLARFP, DLARZ	EXTERNAL DLARFG, DLARZ
* ..	* ..
* .. Executable Statements ..	* .. Executable Statements ..
*	*
Line 112	Line 112
* Generate elementary reflector H(i) to annihilate	* Generate elementary reflector H(i) to annihilate
* [ A(i,i) A(i,n-l+1:n) ]	* [ A(i,i) A(i,n-l+1:n) ]
*	*
CALL DLARFP( L+1, A( I, I ), A( I, N-L+1 ), LDA, TAU( I ) )	CALL DLARFG( L+1, A( I, I ), A( I, N-L+1 ), LDA, TAU( I ) )
*	*
* Apply H(i) to A(1:i-1,i:n) from the right	* Apply H(i) to A(1:i-1,i:n) from the right
*	*