114 SUBROUTINE cptts2( IUPLO, N, NRHS, D, E, B, LDB )
122 INTEGER iuplo, ldb, n, nrhs
126 COMPLEX b( ldb, * ), e( * )
146 $ CALL
csscal( nrhs, 1. / d( 1 ),
b, ldb )
150 IF( iuplo.EQ.1 )
THEN
162 b( i,
j ) =
b( i,
j ) -
b( i-1,
j )*conjg( e( i-1 ) )
168 b( i,
j ) =
b( i,
j ) / d( i )
170 DO 30 i = n - 1, 1, -1
171 b( i,
j ) =
b( i,
j ) -
b( i+1,
j )*e( i )
183 b( i,
j ) =
b( i,
j ) -
b( i-1,
j )*conjg( e( i-1 ) )
188 b( n,
j ) =
b( n,
j ) / d( n )
189 DO 50 i = n - 1, 1, -1
190 b( i,
j ) =
b( i,
j ) / d( i ) -
b( i+1,
j )*e( i )
206 b( i,
j ) =
b( i,
j ) -
b( i-1,
j )*e( i-1 )
212 b( i,
j ) =
b( i,
j ) / d( i )
214 DO 90 i = n - 1, 1, -1
215 b( i,
j ) =
b( i,
j ) -
b( i+1,
j )*conjg( e( i ) )
227 b( i,
j ) =
b( i,
j ) -
b( i-1,
j )*e( i-1 )
232 b( n,
j ) =
b( n,
j ) / d( n )
233 DO 110 i = n - 1, 1, -1
234 b( i,
j ) =
b( i,
j ) / d( i ) -
235 $
b( i+1,
j )*conjg( e( i ) )
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine cptts2(IUPLO, N, NRHS, D, E, B, LDB)
CPTTS2 solves a tridiagonal system of the form AX=B using the L D LH factorization computed by spttrf...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine csscal(N, SA, CX, INCX)
CSSCAL