rpl/lapack/lapack/dgtts2.f - diff

Return to dgtts2.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.7, 2010/12/21 13:53:27	version 1.8, 2011/07/22 07:38:05
Line 17	Line 17
* =======	* =======
*	*
* DGTTS2 solves one of the systems of equations	* DGTTS2 solves one of the systems of equations
* AX = B or A'X = B,	* AX = B or ATX = B,
* with a tridiagonal matrix A using the LU factorization computed	* with a tridiagonal matrix A using the LU factorization computed
* by DGTTRF.	* by DGTTRF.
*	*
Line 27	Line 27
* ITRANS (input) INTEGER	* ITRANS (input) INTEGER
* Specifies the form of the system of equations.	* Specifies the form of the system of equations.
* = 0: A * X = B (No transpose)	* = 0: A * X = B (No transpose)
* = 1: A'* X = B (Transpose)	* = 1: A*T X = B (Transpose)
* = 2: A'* X = B (Conjugate transpose = Transpose)	* = 2: A*T X = B (Conjugate transpose = Transpose)
*	*
* N (input) INTEGER	* N (input) INTEGER
* The order of the matrix A.	* The order of the matrix A.
Line 138	Line 138
END IF	END IF
ELSE	ELSE
*	*
* Solve A' * X = B.	* Solve A*T X = B.
*	*
IF( NRHS.LE.1 ) THEN	IF( NRHS.LE.1 ) THEN
*	*
* Solve U'*x = b.	* Solve U*Tx = b.
*	*
J = 1	J = 1
70 CONTINUE	70 CONTINUE
Line 154	Line 154
$ B( I-2, J ) ) / D( I )	$ B( I-2, J ) ) / D( I )
80 CONTINUE	80 CONTINUE
*	*
* Solve L'*x = b.	* Solve L*Tx = b.
*	*
DO 90 I = N - 1, 1, -1	DO 90 I = N - 1, 1, -1
IP = IPIV( I )	IP = IPIV( I )
Line 170	Line 170
ELSE	ELSE
DO 120 J = 1, NRHS	DO 120 J = 1, NRHS
*	*
* Solve U'*x = b.	* Solve U*Tx = b.
*	*
B( 1, J ) = B( 1, J ) / D( 1 )	B( 1, J ) = B( 1, J ) / D( 1 )
IF( N.GT.1 )	IF( N.GT.1 )