Developer Reference

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

?ppcon

Estimates the reciprocal of the condition number of a packed symmetric (Hermitian) positive-definite matrix.

Syntax

call sppcon
(
uplo
,
n
,
ap
,
anorm
,
rcond
,
work
,
iwork
,
info
)
call dppcon
(
uplo
,
n
,
ap
,
anorm
,
rcond
,
work
,
iwork
,
info
)
call cppcon
(
uplo
,
n
,
ap
,
anorm
,
rcond
,
work
,
rwork
,
info
)
call zppcon
(
uplo
,
n
,
ap
,
anorm
,
rcond
,
work
,
rwork
,
info
)
call ppcon
(
ap
,
anorm
,
rcond
[
,
uplo
]
[
,
info
]
)
Include Files
  • mkl.fi
    ,
    lapack.f90
Description
The routine estimates the reciprocal of the condition number of a packed symmetric (Hermitian) positive-definite matrix
A
:
κ
1
(
A
) = ||
A
||
1
||
A
-1
||
1
(since
A
is symmetric or Hermitian,
κ
(
A
) =
κ
1
(
A
)
).
An estimate is obtained for
||
A
-1
||
, and the reciprocal of the condition number is computed as
rcond
= 1 / (||
A
|| ||
A
-1
||)
.
Before calling this routine:
  • compute
    anorm
    (either
    ||
    A
    ||
    1
    = max
    j
    Σ
    i
    |
    a
    i
    j
    |
    or
    ||
    A
    ||
    = max
    i
    Σ
    j
    |
    a
    i
    j
    |)
  • call
    ?pptrf
    to compute the Cholesky factorization of
    A
    .
Input Parameters
n
INTEGER
.
The order of the matrix
A
;
n
0.
ap
,
work
REAL
for
sppcon
DOUBLE PRECISION
for
dppcon
COMPLEX
for
cppcon
DOUBLE COMPLEX
for
zppcon
.
Arrays:
ap
(*)
,
work
(*)
.
The array
ap
contains the packed factored matrix
A
, as returned by
?pptrf
. The dimension of
ap
must be at least max(1,
n
(
n
+1)/2).
The array
work
is a workspace for the routine. The dimension of
work
must be at least
max(1, 3*
n
)
for real flavors and
max(1, 2*
n
)
for complex flavors.
anorm
REAL
for single precision flavors
DOUBLE PRECISION
for double precision flavors.
The norm of the