CSYT01(3)
π CSYT01(3) β like, the whole entire manual page, bestie
Okay so first things first, the title is CSYT01(3) and the manual it lives in is Library Functions Manual β and the source file is this absolute journey of a path, no cap:
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/LIN/csyt01.f
That's it. That's the file. She's long, she's slaying, she's in the LAPACK testing directory. Periodt. π
β¨ SYNOPSIS β¨
Functions/Subroutines β the cast π
subroutine CSYT01 (uplo, n, a, lda, afac, ldafac, ipiv, c, ldc, rwork,
resid)
CSYT01
Literally one subroutine, one name, ten arguments. Main character energy fr fr. π
πΈ Function/Subroutine Documentation πΈ
subroutine CSYT01 π¦
Here's the full signature so you can see exactly what she's working with:
subroutine CSYT01 (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)
Yep, still CSYT01. Iconic.
π― Purpose π―
So what does she do? Bestie, she reconstructs a complex symmetric indefinite matrix A from its block L*D*L' or U*D*U' factorization, and 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 transpose of L, and U' is the transpose of U.
That's the whole vibe check. Rebuild the matrix, measure how far off you are, divide by N * norm(A) * EPS. Math girlies rise up. β¨
π§ Parameters π§ (every single one, don't worry)
UPLOβ isCHARACTER*1. It specifies whether the upper or lower triangular part of the complex symmetric matrixAis stored:= 'U': Upper triangular= 'L': Lower triangular
Nβ isINTEGER. The number of rows and columns of the matrixA. Constraint:N >= 0.Aβ isCOMPLEX array, dimension (LDA,N). The original complex symmetric matrixA. The source material, the blueprint, the original before the glow-up. πLDAβ isINTEGER. The leading dimension of the arrayA. Constraint:LDA >= max(1,N).AFACβ isCOMPLEX array, 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β isINTEGER. The leading dimension of the arrayAFAC. Constraint:LDAFAC >= max(1,N).IPIVβ isINTEGER array, dimension (N). The pivot indices fromCSYTRF. She's not inventing these, she's receiving them. π«ΆCβ isCOMPLEX array, dimension (LDC,N). That's the whole description. She's the reconstructed one, and honestly she speaks for herself.LDCβ isINTEGER. The leading dimension of the arrayC. Constraint:LDC >= max(1,N).RWORKβ isREAL array, dimension (N). That's it. That's the whole entry. Scratch space doing scratch space things. πRESIDβ isREAL. The actual result, and it depends onUPLO:- If
UPLO = 'L', it'snorm(L*D*L' - A) / ( N * norm(A) * EPS ) - If
UPLO = 'U', it'snorm(U*D*U' - A) / ( N * norm(A) * EPS )
- If
Heads up bestie, theRESIDdescriptions in the original are written in terms of the factorization minus A while the Purpose section says the residual isnorm( C - A ) / ( N * norm(A) * EPS )β both are in there, so I'm giving you both, exactly as written. No editorializing, no cap. π
π Author π
The whole squad of institutions behind this:
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Legends, all of them. Collaborating across state lines to bring you linear algebra. We love to see it. π
Definition at line 123 of file csyt01.f. Line one-two-three, look at her go. π₯
π« Author (again, but the meta one) π«
This page was generated automatically by Doxygen for LAPACK from the source code. Not a human, not a group chat, a documentation generator. Respect the machine. β
Alright bestie, that's the page β every flag, every array, every constraint, all accounted for. ππ
Sign-off time: this was LAPACK, Version 3.12.0, dated/home/abui...IN/csyt01.f(3)β okay the date field is genuinely just that truncated path, but she's doing her best. Bye, ily, go callCSYTRFfirst. π§π
π¦ 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