rpl/lapack/lapack/dptrfs.f - diff

Return to dptrfs.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.7, 2010/12/21 13:53:37	version 1.8, 2011/07/22 07:38:10
Line 1	Line 1
SUBROUTINE DPTRFS( N, NRHS, D, E, DF, EF, B, LDB, X, LDX, FERR,	SUBROUTINE DPTRFS( N, NRHS, D, E, DF, EF, B, LDB, X, LDX, FERR,
$ BERR, WORK, INFO )	$ BERR, WORK, 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, LDB, LDX, N, NRHS	INTEGER INFO, LDB, LDX, N, NRHS
Line 263	Line 263
* 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)*T.
*	*
* Solve M(L) * x = e.	* Solve M(L) * x = e.
*	*
Line 272	Line 272
WORK( I ) = ONE + WORK( I-1 )*ABS( EF( I-1 ) )	WORK( I ) = ONE + WORK( I-1 )*ABS( EF( I-1 ) )
60 CONTINUE	60 CONTINUE
*	*
* Solve D * M(L)' * x = b.	* Solve D * M(L)*T x = b.
*	*
WORK( N ) = WORK( N ) / DF( N )	WORK( N ) = WORK( N ) / DF( N )
DO 70 I = N - 1, 1, -1	DO 70 I = N - 1, 1, -1