Generic Decomposition Operations<T>.Hermitian Estimate Condition Method
Estimates the reciprocal of the condition number of a complex Hermitian matrix A using the factorization A = U*D*UH or A = L*D*LH computed by ZHETRF.
Definition
Namespace: Extreme.Mathematics.LinearAlgebra.Implementation
Assembly: Extreme.Numerics.Generic (in Extreme.Numerics.Generic.dll) Version: 8.1.4
C#
Assembly: Extreme.Numerics.Generic (in Extreme.Numerics.Generic.dll) Version: 8.1.4
public override void HermitianEstimateCondition(
MatrixTriangle storedTriangle,
int n,
Array2D<Complex<T>> a,
Array1D<int> ipiv,
T anorm,
out T rcond,
out int info
)
Parameters
- storedTriangle MatrixTriangle
Specifies whether the details of the factorization are stored as an upper or lower triangular matrix. = 'U': Upper triangular, form is A = U*D*UH; = 'L': Lower triangular, form is A = L*D*LH.
- n Int32
The order of the matrix A. N >= 0.
- a Array2D<Complex<T>>
A is complex array, dimension (LDA,N) The block diagonal matrix D and the multipliers used to obtain the factor U or L as computed by ZHETRF.
The leading dimension of the array A. LDA >= max(1,N).
- ipiv Array1D<Int32>
Dimension (N) Details of the interchanges and the block structure of D as determined by ZHETRF.
- anorm T
The 1-norm of the original matrix A.
- rcond T
The reciprocal of the condition number of the matrix A, computed as RCOND = 1/(ANORM * AINVNM), where AINVNM is an estimate of the 1-norm of inv(A) computed in this routine.
- info Int32
info is INTEGER = 0: successful exit < 0: if info = -i, the i-th argument had an illegal value
Remarks
An estimate is obtained for norm(inv(A)), and the reciprocal of the condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))).