rpl/lapack/lapack/dorgbr.f - diff

Return to dorgbr.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.7, 2010/12/21 13:53:34	version 1.8, 2011/07/22 07:38:08
Line 201	Line 201
END IF	END IF
ELSE	ELSE
*	*
* Form P', determined by a call to DGEBRD to reduce a k-by-n	* Form P**T, determined by a call to DGEBRD to reduce a k-by-n
* matrix	* matrix
*	*
IF( K.LT.N ) THEN	IF( K.LT.N ) THEN
Line 215	Line 215
* 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**T to
* those of the unit matrix	* those of the unit matrix
*	*
A( 1, 1 ) = ONE	A( 1, 1 ) = ONE
Line 230	Line 230
60 CONTINUE	60 CONTINUE
IF( N.GT.1 ) THEN	IF( N.GT.1 ) THEN
*	*
* Form P'(2:n,2:n)	* Form P**T(2:n,2:n)
*	*
CALL DORGLQ( N-1, N-1, N-1, A( 2, 2 ), LDA, TAU, WORK,	CALL DORGLQ( N-1, N-1, N-1, A( 2, 2 ), LDA, TAU, WORK,
$ LWORK, IINFO )	$ LWORK, IINFO )