TESTING/LIN/strt01.f(3)
π TESTING/LIN/strt01.f π
Manual name tea: TESTING/LIN/strt01.f(3) from the Library Functions Manual. This page is literally about TESTING/LIN/strt01.f, no cap. π
πΈ NAME πΈ
TESTING/LIN/strt01.f β bestie, thatβs the name. This manual page is TESTING/LIN/strt01.f(3), and itβs here to explain the strt01 subroutine in all its triangular-inverse-residual glory. β¨
β¨ SYNOPSIS β¨
π«Ά Functions/Subroutines π«Ά
subroutine strt01 (uplo, diag, n, a, lda, ainv, ldainv, rcond, work, resid)βSTRT01
π Function/Subroutine Documentation π
subroutine strt01 (character uplo, character diag, integer n, real, dimension( lda, * ) a, integer lda, real, dimension( ldainv, * ) ainv, integer ldainv, real rcond, real, dimension( * ) work, real resid)
STRT01
π Purpose π
Okay so STRT01 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. Slay. π«
π§ Parameters π§
UPLOβ isCHARACTER*1. Specifies whether the matrix A is upper or lower triangular.= 'U': Upper triangular.= 'L': Lower triangular.DIAGβ isCHARACTER*1. Specifies whether or not the matrix A is unit triangular.= 'N': Non-unit triangular.= 'U': Unit triangular.Nβ isINTEGER. The order of the matrix A.N >= 0.Aβ isREAL array, dimension (LDA,N). The triangular matrix A. IfUPLO = 'U', the leading n by n upper triangular part of the array A contains the upper triangular matrix, and the strictly lower triangular part of A is not referenced. IfUPLO = 'L', the leading n by n lower triangular part of the array A contains the lower triangular matrix, and the strictly upper triangular part of A is not referenced. IfDIAG = 'U', the diagonal elements of A are also not referenced and are assumed to be1.LDAβ isINTEGER. The leading dimension of the array A.LDA >= max(1,N).AINVβ isREAL array, dimension (LDAINV,N). On entry, the (triangular) inverse of the matrix A, in the same storage format as A. On exit, the contents of AINV are destroyed.LDAINVβ isINTEGER. The leading dimension of the array AINV.LDAINV >= max(1,N).RCONDβ isREAL. The reciprocal condition number of A, computed as1/(norm(A) * norm(AINV)).WORKβ isREAL array, dimension (N).RESIDβ isREAL. Its value isnorm(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 122 of file strt01.f.
π Author π
This page was Generated automatically by Doxygen for LAPACK from the source code. β¨
Okay bestie, final footer tea: source/version is LAPACK, Version 3.12.0, and the page id is TESTING/LIN/strt01.f(3). No date is shown on this footer, just so you know. πππ¦
π¦ 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