rpl/lapack/lapack/zgglse.f - diff

Return to zgglse.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.17, 2018/05/29 07:18:17	version 1.18, 2023/08/07 08:39:21
Line 172	Line 172
*> \author Univ. of Colorado Denver	*> \author Univ. of Colorado Denver
*> \author NAG Ltd.	*> \author NAG Ltd.
*	*
*> \date December 2016
*
*> \ingroup complex16OTHERsolve	*> \ingroup complex16OTHERsolve
*	*
* =====================================================================	* =====================================================================
SUBROUTINE ZGGLSE( M, N, P, A, LDA, B, LDB, C, D, X, WORK, LWORK,	SUBROUTINE ZGGLSE( M, N, P, A, LDA, B, LDB, C, D, X, WORK, LWORK,
$ INFO )	$ INFO )
*	*
* -- LAPACK driver routine (version 3.7.0) --	* -- LAPACK driver routine --
* -- 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..--
* December 2016
*	*
* .. Scalar Arguments ..	* .. Scalar Arguments ..
INTEGER INFO, LDA, LDB, LWORK, M, N, P	INTEGER INFO, LDA, LDB, LWORK, M, N, P
Line 279	Line 276
*	*
CALL ZGGRQF( P, M, N, B, LDB, WORK, A, LDA, WORK( P+1 ),	CALL ZGGRQF( P, M, N, B, LDB, WORK, A, LDA, WORK( P+1 ),
$ WORK( P+MN+1 ), LWORK-P-MN, INFO )	$ WORK( P+MN+1 ), LWORK-P-MN, INFO )
LOPT = WORK( P+MN+1 )	LOPT = INT( WORK( P+MN+1 ) )
*	*
* Update c = Z*H c = ( c1 ) N-P	* Update c = Z*H c = ( c1 ) N-P
* ( c2 ) M+P-N	* ( c2 ) M+P-N