Developer Reference for Intel® oneAPI Math Kernel Library for Fortran

ID 766686
Date 11/07/2023
Public

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

Document Table of Contents

?la_syrcond

Estimates the Skeel condition number for a symmetric indefinite matrix.

Syntax

call sla_syrcond( uplo, n, a, lda, af, ldaf, ipiv, cmode, c, info, work, iwork )

call dla_syrcond( uplo, n, a, lda, af, ldaf, ipiv, cmode, c, info, work, iwork )

Include Files

  • mkl.fi

Description

The function estimates the Skeel condition number of

op(A) * op2(C)

where

the cmode parameter determines op2 as follows:

cmode Value

op2(C)

1

C

0

I

-1

inv(C)

The Skeel condition number

cond(A) = norminf(|inv(A)||A|)

is computed by computing scaling factors R such that

diag(R)*A*op2(C)

is row equilibrated and by computing the standard infinity-norm condition number.

Input Parameters

uplo

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

Specifies the triangle of A to store:

If uplo = 'U', the upper triangle of A is stored,

If uplo = 'L', the lower triangle of A is stored.

n

INTEGER. The number of linear equations, that is, the order of the matrix A; n≥ 0.

a, af, c, work

REAL for sla_syrcond

DOUBLE PRECISION for dla_syrcond

Arrays:

ab (lda,*) contains the n-by-n matrix A.

af (ldaf,*) contains the The block diagonal matrix D and the multipliers used to obtain the factor L or U as computed by ?sytrf.

The second dimension of a and af must be at least max(1, n).

c, DIMENSIONn. The vector C in the formula op(A) * op2(C).

work is a workspace array of DIMENSION (3*n).

lda

INTEGER. The leading dimension of the array ab. lda≥max(1,n).

ldaf

INTEGER. The leading dimension of af. ldaf≥max(1,n).

ipiv

INTEGER.

Array with DIMENSIONn. Details of the interchanges and the block structure of D as determined by ?sytrf.

cmode

INTEGER. Determines op2(C) in the formula op(A) * op2(C) as follows:

If cmode = 1, op2(C) = C.

If cmode = 0, op2(C) = I.

If cmode = -1, op2(C) = inv(C).

iwork

INTEGER. Workspace array with DIMENSIONn.

Output Parameters

info

INTEGER.

If info = 0, the execution is successful.

If i > 0, the i-th parameter is invalid.

See Also