CHET01(3)
π NAME π
okay bestie, you're reading CHET01(3) from the Library Functions Manual, and she is here to reconstruct a Hermitian indefinite matrix and check the residual π
. The actual source file path is:
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/LIN/chet01.f
β¨ SYNOPSIS β¨
π§ Functions/Subroutines π§
subroutine CHET01 (uplo, n, a, lda, afac, ldafac, ipiv, c, ldc, rwork,
resid)
CHET01
π¦ Function/Subroutine Documentation π¦
subroutine CHET01 (character uplo, integer n, complex, dimension( lda, * )
a, integer lda, complex, dimension( ldafac, * ) afac, integer ldafac,
integer, dimension( * ) ipiv, complex, dimension( ldc, * ) c, integer
ldc, real, dimension( * ) rwork, real resid)
CHET01
π Purpose π
CHET01 reconstructs a Hermitian indefinite matrix A from its block L*D*L' or U*D*U' factorization, then 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. literally iconic math behaviour π
β¨
πΈ Parameters πΈ
UPLOβUPLOisCHARACTER*1. It says which triangle of the Hermitian matrixAis stored:'U'means upper triangular, and'L'means lower triangular. Pick your lane, bestie πNβNisINTEGER. It's the number of rows and columns of the matrixA, andN >= 0. no negative-size nonsense, periodt π«ΆAβAisCOMPLEXarray, dimension(LDA,N). This is the original Hermitian matrixA. main character matrix energy πLDAβLDAisINTEGER. It's the leading dimension of the arrayA, andLDA >= max(1,N). respect the leading dimension or perish πAFACβAFACisCOMPLEXarray, dimension(LDAFAC,N). This is the factored form of the matrixA.AFACcontains the block diagonal matrixDand the multipliers used to obtain the factorLorUfrom the blockL*D*L'orU*D*U'factorization as computed byCHETRF. yes, she contains lore πLDAFACβLDAFACisINTEGER. It's the leading dimension of the arrayAFAC, andLDAFAC >= max(1,N).get that leading dimension in check bestie πIPIVβIPIVisINTEGERarray, dimension(N). These are the pivot indices fromCHETRF. pivot era π₯CβCisCOMPLEXarray, dimension(LDC,N). that's the reconstructed-matrix slot, no extra commentary needed π¦LDCβLDCisINTEGER. It's the leading dimension of the arrayC, andLDC >= max(1,N). another leading dimension, another slay πRWORKβRWORKisREALarray, dimension(N). real workspace, real vibes β¨RESIDβRESIDisREAL. IfUPLO = 'L', it'snorm(L*D*L' - A) / ( N * norm(A) * EPS ). IfUPLO = 'U', it'snorm(U*D*U' - A) / ( N * norm(A) * EPS ). this is the final boss residual check π€
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 124 of file chet01.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
bye bestie π β this page is from LAPACK Version 3.12.0, and the footer path is /home/abui...IN/chet01.f(3). No date is listed in the footer, so consider this your timestamp-free sign-off ππ
π¦ where this page comes from
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms