?pttrsv

Developer Reference for Intel® oneAPI Math Kernel Library for Fortran

Download PDF

ID 766686

Date 3/31/2023

Version

Public

A newer version of this document is available. Customers should click here to go to the newest version.

Visible to Intel only — GUID: GUID-ED5CFCFB-DD13-4AD5-95DA-8AADBAD6F708

View Details

?pttrsv

Solves a symmetric (Hermitian) positive-definite tridiagonal system of linear equations, using the L*D*L^H factorization computed by ?pttrf.

Syntax

call spttrsv(trans, n, nrhs, d, e, b, ldb, info)

call dpttrsv(trans, n, nrhs, d, e, b, ldb, info)

call cpttrsv(uplo, trans, n, nrhs, d, e, b, ldb, info)

call zpttrsv(uplo, trans, n, nrhs, d, e, b, ldb, info)

Description

The ?pttrsvroutine solves one of the triangular systems:

L^T*X = B, or L*X = B for real flavors,

L*X = B, or L^H*X = B,

U*X = B, or U^H*X = B for complex flavors,

where L (or U for complex flavors) is the Cholesky factor of a Hermitian positive-definite tridiagonal matrix A such that

A = L*D*L^H (computed by spttrf/dpttrf)

A = U^H*D*U or A = L*D*L^H (computed by cpttrf/zpttrf).

Input Parameters

uplo

CHARACTER*1. Must be 'U' or 'L'.

Specifies whether the superdiagonal or the subdiagonal of the tridiagonal matrix A is stored and the form of the factorization:

If uplo = 'U', e is the superdiagonal of U, and A = U^H*D*U or A = L*D*L^H;

if uplo = 'L', e is the subdiagonal of L, and A = L*D*L^H.

The two forms are equivalent, if A is real.

trans

CHARACTER.

Specifies the form of the system of equations:

for real flavors:

if trans = 'N': L*X = B (no transpose)

if trans = 'T': L^T*X = B (transpose)

for complex flavors:

if trans = 'N': U*X = B or L*X = B (no transpose)

if trans = 'C': U^H*X = B or L^H*X = B (conjugate transpose).

n

INTEGER. The order of the tridiagonal matrix A. n ≥ 0.

nrhs

INTEGER. The number of right hand sides, that is, the number of columns of the matrix B. nrhs ≥ 0.

d

REAL array of size n. The n diagonal elements of the diagonal matrix D from the factorization computed by ?pttrf.

e

COMPLEX array of size (n-1). The (n-1) off-diagonal elements of the unit bidiagonal factor U or L from the factorization computed by ?pttrf. See uplo.

b

COMPLEX array of size ldbby nrhs.

On entry, the right hand side matrix B.

ldb

INTEGER.

The leading dimension of the array b.

ldb ≥ max(1, n).

Output Parameters

b

On exit, the solution matrix X.

info

INTEGER.

= 0: successful exit

< 0: if info = -i, the i-th argument had an illegal value.

Parent topic: ScaLAPACK Auxiliary Routines

Select Your Language

Using Intel.com Search

Quick Links

Recent Searches

Advanced Search

Only search in

Developer Reference for Intel® oneAPI Math Kernel Library for Fortran

?pttrsv

See Also