CTRT01(3)
Hey bestie π welcome to CTRT01(3), hanging out in the Library Functions Manual. This is a LAPACK testing-lin situation, and the vibe is: check the residual of a triangular matrix times its inverse. No cap.
πΈ NAME πΈ
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/LIN/ctrt01.f
β¨ SYNOPSIS β¨
π Functions/Subroutines π
subroutine CTRT01 (uplo, diag, n, a, lda, ainv, ldainv, rcond, rwork, resid)
CTRT01
π Function/Subroutine Documentation π
π¦ subroutine CTRT01 π¦
subroutine CTRT01 (character uplo, character diag, integer n, complex, dimension( lda, * ) a, integer lda, complex, dimension( ldainv, * ) ainv, integer ldainv, real rcond, real, dimension( * ) rwork, real resid)
CTRT01
β¨ Purpose β¨
Okay so CTRT01 computes the residual for a triangular matrix A times its inverse:
RESID = norm( A*AINV - I ) / ( N * norm(A) * norm(AINV) * EPS ),
where EPS is the machine epsilon. Obsessed with a clean residual, honestly.
π Parameters π
UPLOβUPLOisCHARACTER*1. Specifies whether the matrixAis upper or lower triangular.= 'U': Upper triangular= 'L': Lower triangular
DIAGβDIAGisCHARACTER*1. Specifies whether or not the matrixAis unit triangular.= 'N': Non-unit triangular= 'U': Unit triangular
NβNisINTEGER. The order of the matrixA.N >= 0.AβAisCOMPLEX array, dimension (LDA,N). The triangular matrixA. IfUPLO = 'U', the leading n by n upper triangular part of the arrayAcontains the upper triangular matrix, and the strictly lower triangular part ofAis not referenced. IfUPLO = 'L', the leading n by n lower triangular part of the arrayAcontains the lower triangular matrix, and the strictly upper triangular part ofAis not referenced. IfDIAG = 'U', the diagonal elements ofAare also not referenced and are assumed to be 1.LDAβLDAisINTEGER. The leading dimension of the arrayA.LDA >= max(1,N).AINVβAINVisCOMPLEX array, dimension (LDAINV,N). On entry, the (triangular) inverse of the matrixA, in the same storage format asA. On exit, the contents ofAINVare destroyed.LDAINVβLDAINVisINTEGER. The leading dimension of the arrayAINV.LDAINV >= max(1,N).RCONDβRCONDisREAL. The reciprocal condition number ofA, computed as1/(norm(A) * norm(AINV)).RWORKβRWORKisREAL array, dimension (N).RESIDβRESIDisREAL.norm(A*AINV - I) / ( N * norm(A) * norm(AINV) * EPS )
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 123 of file ctrt01.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
Thatβs the whole page, bestie π
. Footer tea: source/version LAPACK Version 3.12.0 and the path /home/abui...IN/ctrt01.f(3). The original footer doesnβt print a date, so no invented date β periodt π«Άβ¨
π¦ 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