Developer Reference

  • 2020.2
  • 07/15/2020
  • Public Content
Contents

?sptrs

Solves a system of linear equations with a UDU- or LDL-factored symmetric coefficient matrix using packed storage.

Syntax

call ssptrs
(
uplo
,
n
,
nrhs
,
ap
,
ipiv
,
b
,
ldb
,
info
)
call dsptrs
(
uplo
,
n
,
nrhs
,
ap
,
ipiv
,
b
,
ldb
,
info
)
call csptrs
(
uplo
,
n
,
nrhs
,
ap
,
ipiv
,
b
,
ldb
,
info
)
call zsptrs
(
uplo
,
n
,
nrhs
,
ap
,
ipiv
,
b
,
ldb
,
info
)
call sptrs
(
ap
,
b
,
ipiv
[
,
uplo
]
[
,
info
]
)
Include Files
  • mkl.fi
    ,
    lapack.f90
Description
The routine solves for
X
the system of linear equations
A*X
=
B
with a symmetric matrix
A
, given the Bunch-Kaufman factorization of
A
:
if
uplo
=
'U'
,
A
=
U*D*U
T
if
uplo
=
'L'
,
A
=
L*D*L
T
,
where
U
and
L
are upper and lower packed triangular matrices with unit diagonal and
D
is a symmetric block-diagonal matrix. The system is solved with multiple right-hand sides stored in the columns of the matrix
B
. You must supply the factor
U
(or
L
) and the array
ipiv
returned by the factorization routine
?sptrf
.
Input Parameters
uplo
CHARACTER*1
.
Must be
'U'
or
'L'
.
Indicates how the input matrix
A
has been factored:
If
uplo
=
'U'
, the array
ap
stores the packed factor
U
of the factorization
A
=
U*D*U
T
. If
uplo
=
'L'
, the array
ap
stores the packed factor
L
of the factorization
A
=
L*D*L
T
.
n
INTEGER
.
The order of matrix
A
;
n
0.
nrhs
INTEGER
.
The number of right-hand sides;
nrhs
0.
ipiv
INTEGER
.
Array, size at least
max(1,
n
)
. The
ipiv
array, as returned by
?sptrf
.
ap
REAL
for
ssptrs
DOUBLE PRECISION
for
dsptrs
COMPLEX
for
csptrs
DOUBLE COMPLEX
for
zsptrs
.
The dimension of array
ap
must be at least max(1,
n
(
n
+1)/2). The array
ap
contains the factor
U
or
L
, as specified by
uplo
, in