rpl/lapack/lapack/zgbcon.f - diff

Return to zgbcon.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.7, 2010/12/21 13:53:42	version 1.8, 2011/07/22 07:38:13
Line 1	Line 1
SUBROUTINE ZGBCON( NORM, N, KL, KU, AB, LDAB, IPIV, ANORM, RCOND,	SUBROUTINE ZGBCON( NORM, N, KL, KU, AB, LDAB, IPIV, ANORM, RCOND,
$ WORK, RWORK, INFO )	$ WORK, RWORK, 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 --
*	*
* Modified to call ZLACN2 in place of ZLACON, 10 Feb 03, SJH.	* Modified to call ZLACN2 in place of ZLACON, 10 Feb 03, SJH.
*	*
Line 187	Line 187
$ KL+KU, AB, LDAB, WORK, SCALE, RWORK, INFO )	$ KL+KU, AB, LDAB, WORK, SCALE, RWORK, INFO )
ELSE	ELSE
*	*
* Multiply by inv(U').	* Multiply by inv(U**H).
*	*
CALL ZLATBS( 'Upper', 'Conjugate transpose', 'Non-unit',	CALL ZLATBS( 'Upper', 'Conjugate transpose', 'Non-unit',
$ NORMIN, N, KL+KU, AB, LDAB, WORK, SCALE, RWORK,	$ NORMIN, N, KL+KU, AB, LDAB, WORK, SCALE, RWORK,
$ INFO )	$ INFO )
*	*
* Multiply by inv(L').	* Multiply by inv(L**H).
*	*
IF( LNOTI ) THEN	IF( LNOTI ) THEN
DO 30 J = N - 1, 1, -1	DO 30 J = N - 1, 1, -1