rpl/lapack/lapack/zungbr.f - diff

Return to zungbr.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.7, 2010/12/21 13:53:58	version 1.8, 2011/07/22 07:38:21
Line 202	Line 202
END IF	END IF
ELSE	ELSE
*	*
* Form P', determined by a call to ZGEBRD to reduce a k-by-n	* Form P**H, determined by a call to ZGEBRD to reduce a k-by-n
* matrix	* matrix
*	*
IF( K.LT.N ) THEN	IF( K.LT.N ) THEN
Line 216	Line 216
* If k >= n, assume m = n	* If k >= n, assume m = n
*	*
* Shift the vectors which define the elementary reflectors one	* Shift the vectors which define the elementary reflectors one
* row downward, and set the first row and column of P' to	* row downward, and set the first row and column of P**H to
* those of the unit matrix	* those of the unit matrix
*	*
A( 1, 1 ) = ONE	A( 1, 1 ) = ONE
Line 231	Line 231
60 CONTINUE	60 CONTINUE
IF( N.GT.1 ) THEN	IF( N.GT.1 ) THEN
*	*
* Form P'(2:n,2:n)	* Form P**H(2:n,2:n)
*	*
CALL ZUNGLQ( N-1, N-1, N-1, A( 2, 2 ), LDA, TAU, WORK,	CALL ZUNGLQ( N-1, N-1, N-1, A( 2, 2 ), LDA, TAU, WORK,
$ LWORK, IINFO )	$ LWORK, IINFO )