Developer Reference

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

?sytrs

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

Syntax

call ssytrs
(
uplo
,
n
,
nrhs
,
a
,
lda
,
ipiv
,
b
,
ldb
,
info
)
call dsytrs
(
uplo
,
n
,
nrhs
,
a
,
lda
,
ipiv
,
b
,
ldb
,
info
)
call csytrs
(
uplo
,
n
,
nrhs
,
a
,
lda
,
ipiv
,
b
,
ldb
,
info
)
call zsytrs
(
uplo
,
n
,
nrhs
,
a
,
lda
,
ipiv
,
b
,
ldb
,
info
)
call sytrs
(
a
,
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 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 to this routine the factor
U
(or
L
) and the array
ipiv
returned by the factorization routine
?sytrf
.
Input Parameters
uplo
CHARACTER*1
.
Must be
'U'
or
'L'
.
Indicates how the input matrix
A
has been factored:
If
uplo
=
'U'
, the array
a
stores the upper triangular factor
U
of the factorization
A
=
U*D*U
T
.
If
uplo
=
'L'
, the array
a
stores the lower triangular 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
?sytrf
.
a
,
b
REAL
for
ssytrs
DOUBLE PRECISION
for
dsytrs
COMPLEX
for
csytrs
DOUBLE COMPLEX
for