rpl/lapack/lapack/dppcon.f - diff

Return to dppcon.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.7, 2010/12/21 13:53:36	version 1.8, 2011/07/22 07:38:10
Line 1	Line 1
SUBROUTINE DPPCON( UPLO, N, AP, ANORM, RCOND, WORK, IWORK, INFO )	SUBROUTINE DPPCON( UPLO, N, AP, ANORM, RCOND, WORK, IWORK, 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 DLACN2 in place of DLACON, 5 Feb 03, SJH.	* Modified to call DLACN2 in place of DLACON, 5 Feb 03, SJH.
*	*
Line 128	Line 128
IF( KASE.NE.0 ) THEN	IF( KASE.NE.0 ) THEN
IF( UPPER ) THEN	IF( UPPER ) THEN
*	*
* Multiply by inv(U').	* Multiply by inv(U**T).
*	*
CALL DLATPS( 'Upper', 'Transpose', 'Non-unit', NORMIN, N,	CALL DLATPS( 'Upper', 'Transpose', 'Non-unit', NORMIN, N,
$ AP, WORK, SCALEL, WORK( 2*N+1 ), INFO )	$ AP, WORK, SCALEL, WORK( 2*N+1 ), INFO )
Line 146	Line 146
$ AP, WORK, SCALEL, WORK( 2*N+1 ), INFO )	$ AP, WORK, SCALEL, WORK( 2*N+1 ), INFO )
NORMIN = 'Y'	NORMIN = 'Y'
*	*
* Multiply by inv(L').	* Multiply by inv(L**T).
*	*
CALL DLATPS( 'Lower', 'Transpose', 'Non-unit', NORMIN, N,	CALL DLATPS( 'Lower', 'Transpose', 'Non-unit', NORMIN, N,
$ AP, WORK, SCALEU, WORK( 2*N+1 ), INFO )	$ AP, WORK, SCALEU, WORK( 2*N+1 ), INFO )