TESTING/LIN/csyt01_aa.f(3)
πΈ TESTING/LIN/csyt01_aa.f(3) πΈ
okay bestie, welcome to the manual page for TESTING/LIN/csyt01_aa.f β yes, itβs living its best Library Functions Manual life, and we are reading it together. No cap, every detail is staying. π
β¨ NAME β¨
TESTING/LIN/csyt01_aa.f β the file, the icon, the whole manual-page moment. π
π SYNOPSIS π
π« Functions/Subroutines π«
subroutine csyt01_aa (uplo, n, a, lda, afac, ldafac, ipiv, c, ldc,
rwork, resid)
CSYT01 β thatβs the subroutine name, short and iconic. β¨
π¦ Function/Subroutine Documentation π¦
okay so hereβs the full signature, bestie, and it is long but we are not skipping a crumb:
subroutine csyt01_aa (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)
CSYT01 β thatβs the name, periodt. π
π Purpose π
okay so CSYT01 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 and EPS is the machine epsilon. Itβs giving βreconstruct and check,β and honestly? Obsessed. β¨
πΈ Parameters πΈ
UPLOβUPLOisCHARACTER*1. It specifies whether the upper or lower triangular part of the hermitian matrixAis stored:= 'U':Upper triangular= 'L':Lower triangular
NβNisINTEGER. The number of rows and columns of the matrixA.N >= 0.AβAisREALarray, dimension(LDA,N). The original hermitian matrixA.LDAβLDAisINTEGER. The leading dimension of the arrayA.LDA >= max(1,N).AFACβAFACisREALarray, dimension(LDAFAC,N). 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 byCSYTRF.LDAFACβLDAFACisINTEGER. The leading dimension of the arrayAFAC.LDAFAC >= max(1,N).IPIVβIPIVisINTEGERarray, dimension(N). The pivot indices fromCSYTRF.CβCisREALarray, dimension(LDC,N).LDCβLDCisINTEGER. The leading dimension of the arrayC.LDC >= max(1,N).RWORKβRWORKisREALarray, dimension(N).RESIDβRESIDisREAL. IfUPLO = 'L',norm(L*D*L' - A) / ( N * norm(A) * EPS ). IfUPLO = '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 122 of file csyt01_aa.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code. π
π cute sign-off: source is LAPACK, version 3.12.0, manual page TESTING/LIN/csyt01_aa.f(3) β thatβs the footer vibe, bestie! Mwah! π
π¦ 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