TESTING/LIN/zhet01.f(3)
πΈ NAME πΈ
okay bestie, welcome to TESTING/LIN/zhet01.f β that's the page/file we're on, and it's all about the zhet01 subroutine. π«Ά
β¨ SYNOPSIS β¨
Functions/Subroutines
subroutine zhet01 (uplo, n, a, lda, afac, ldafac, ipiv, c, ldc, rwork,
resid)
ZHET01
π Function/Subroutine Documentation π
Here's the full signature, character-for-character, because precision is the vibe:
subroutine zhet01 (character uplo, integer n, complex*16, dimension( lda, *
) a, integer lda, complex*16, dimension( ldafac, * ) afac, integer
ldafac, integer, dimension( * ) ipiv, complex*16, dimension( ldc, * ) c,
integer ldc, double precision, dimension( * ) rwork, double precision
resid)
ZHET01
π Purpose π
ZHET01 reconstructs a Hermitian indefinite matrix A from its block L*D*L' or U*D*U' factorization and computes the residual
norm( C - A ) / ( N * norm(A) * EPS ),
where C is the reconstructed matrix, EPS is the machine epsilon, L' is the conjugate transpose of L, and U' is the conjugate transpose of U. That's the whole glow-up. β¨
π Parameters π
UPLO
UPLO is CHARACTER*1. It specifies whether the upper or lower triangular part of the Hermitian matrix A is stored:
= 'U': Upper triangular= 'L': Lower triangular
N
N is INTEGER. It's the number of rows and columns of the matrix A. N >= 0.
A
A is COMPLEX*16 array, dimension (LDA,N). It's the original Hermitian matrix A.
LDA
LDA is INTEGER. It's the leading dimension of the array A. LDA >= max(1,N).
AFAC
AFAC is COMPLEX*16 array, dimension (LDAFAC,N). It's the factored form of the matrix A. AFAC contains the block diagonal matrix D and the multipliers used to obtain the factor L or U from the block L*D*L' or U*D*U' factorization as computed by ZHETRF.
LDAFAC
LDAFAC is INTEGER. It's the leading dimension of the array AFAC. LDAFAC >= max(1,N).
IPIV
IPIV is INTEGER array, dimension (N). These are the pivot indices from ZHETRF.
C
C is COMPLEX*16 array, dimension (LDC,N).
LDC
LDC is INTEGER. It's the leading dimension of the array C. LDC >= max(1,N).
RWORK
RWORK is DOUBLE PRECISION array, dimension (N).
RESID
RESID is DOUBLE PRECISION. Its value is:
- If
UPLO = 'L':norm(L*D*L' - A) / ( N * norm(A) * EPS ) - If
UPLO = 'U':norm(U*D*U' - A) / ( N * norm(A) * EPS )
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π Definition at line 124 of file zhet01.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
π That's the whole page, bestie. It's LAPACK, Version 3.12.0, and the file is TESTING/LIN/zhet01.f(3) β now go slay those Hermitian matrix checks. π
π¦ where this page comes from
- arch rolling: lapack-doc 3.12.1-2 Β· custom
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms