rpl/lapack/lapack/dorglq.f - diff

Return to dorglq.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 1	Line 1
SUBROUTINE DORGLQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )	SUBROUTINE DORGLQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
*	*
* -- 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 INFO, K, LDA, LWORK, M, N	INTEGER INFO, K, LDA, LWORK, M, N
Line 185	Line 185
CALL DLARFT( 'Forward', 'Rowwise', N-I+1, IB, A( I, I ),	CALL DLARFT( 'Forward', 'Rowwise', N-I+1, IB, A( I, I ),
$ LDA, TAU( I ), WORK, LDWORK )	$ LDA, TAU( I ), WORK, LDWORK )
*	*
* Apply H' to A(i+ib:m,i:n) from the right	* Apply H**T to A(i+ib:m,i:n) from the right
*	*
CALL DLARFB( 'Right', 'Transpose', 'Forward', 'Rowwise',	CALL DLARFB( 'Right', 'Transpose', 'Forward', 'Rowwise',
$ M-I-IB+1, N-I+1, IB, A( I, I ), LDA, WORK,	$ M-I-IB+1, N-I+1, IB, A( I, I ), LDA, WORK,
Line 193	Line 193
$ LDWORK )	$ LDWORK )
END IF	END IF
*	*
* Apply H' to columns i:n of current block	* Apply H**T to columns i:n of current block
*	*
CALL DORGL2( IB, N-I+1, IB, A( I, I ), LDA, TAU( I ), WORK,	CALL DORGL2( IB, N-I+1, IB, A( I, I ), LDA, TAU( I ), WORK,
$ IINFO )	$ IINFO )