rpl/lapack/lapack/zptrfs.f - diff

Return to zptrfs.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.7, 2010/12/21 13:53:54	version 1.8, 2011/07/22 07:38:19
Line 1	Line 1
SUBROUTINE ZPTRFS( UPLO, N, NRHS, D, E, DF, EF, B, LDB, X, LDX,	SUBROUTINE ZPTRFS( UPLO, N, NRHS, D, E, DF, EF, B, LDB, X, LDX,
$ FERR, BERR, WORK, RWORK, INFO )	$ FERR, BERR, 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 --
*	*
* .. Scalar Arguments ..	* .. Scalar Arguments ..
CHARACTER UPLO	CHARACTER UPLO
Line 327	Line 327
* m(i,j) = abs(A(i,j)), i = j,	* m(i,j) = abs(A(i,j)), i = j,
* m(i,j) = -abs(A(i,j)), i .ne. j,	* m(i,j) = -abs(A(i,j)), i .ne. j,
*	*
* and e = [ 1, 1, ..., 1 ]'. Note M(A) = M(L)DM(L)'.	* and e = [ 1, 1, ..., 1 ]*T. Note M(A) = M(L)DM(L)*H.
*	*
* Solve M(L) * x = e.	* Solve M(L) * x = e.
*	*
Line 336	Line 336
RWORK( I ) = ONE + RWORK( I-1 )*ABS( EF( I-1 ) )	RWORK( I ) = ONE + RWORK( I-1 )*ABS( EF( I-1 ) )
70 CONTINUE	70 CONTINUE
*	*
* Solve D * M(L)' * x = b.	* Solve D * M(L)*H x = b.
*	*
RWORK( N ) = RWORK( N ) / DF( N )	RWORK( N ) = RWORK( N ) / DF( N )
DO 80 I = N - 1, 1, -1	DO 80 I = N - 1, 1, -1