rpl/lapack/lapack/zgeqr2.f - diff

Return to zgeqr2.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.8, 2010/12/21 13:53:43	version 1.9, 2011/07/22 07:38:14
Line 1	Line 1
SUBROUTINE ZGEQR2( M, N, A, LDA, TAU, WORK, INFO )	SUBROUTINE ZGEQR2( M, N, A, LDA, TAU, WORK, INFO )
*	*
* -- LAPACK routine (version 3.2.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..--
* June 2010	* -- April 2011 --
*	*
* .. Scalar Arguments ..	* .. Scalar Arguments ..
INTEGER INFO, LDA, M, N	INTEGER INFO, LDA, M, N
Line 57	Line 57
*	*
* Each H(i) has the form	* Each H(i) has the form
*	*
* H(i) = I - tau * v * v'	* H(i) = I - tau * v * v**H
*	*
* where tau is a complex scalar, and v is a complex vector with	* where tau is a complex scalar, and v is a complex vector with
* v(1:i-1) = 0 and v(i) = 1; v(i+1:m) is stored on exit in A(i+1:m,i),	* v(1:i-1) = 0 and v(i) = 1; v(i+1:m) is stored on exit in A(i+1:m,i),
Line 106	Line 106
$ TAU( I ) )	$ TAU( I ) )
IF( I.LT.N ) THEN	IF( I.LT.N ) THEN
*	*
* Apply H(i)' to A(i:m,i+1:n) from the left	* Apply H(i)**H to A(i:m,i+1:n) from the left
*	*
ALPHA = A( I, I )	ALPHA = A( I, I )
A( I, I ) = ONE	A( I, I ) = ONE