TESTING/LIN/dtrt01.f(3)
π NAME π
Okay bestie, this is the manual page for TESTING/LIN/dtrt01.f(3) from the Library Functions Manual β and the whole page is about TESTING/LIN/dtrt01.f. Lowkey iconic.
β¨ SYNOPSIS β¨
Functions/Subroutines
subroutine dtrt01 (uplo, diag, n, a, lda, ainv, ldainv, rcond, work,
resid)
DTRT01 β that's the subroutine, no cap.
π¦ Function/Subroutine Documentation π¦
Signature
subroutine dtrt01 (character uplo, character diag, integer n, double
precision, dimension( lda, * ) a, integer lda, double precision,
dimension( ldainv, * ) ainv, integer ldainv, double precision rcond,
double precision, dimension( * ) work, double precision resid)
DTRT01
Purpose
DTRT01 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. Periodt.
Parameters
UPLOβCHARACTER*1. Specifies whether the matrix A is upper or lower triangular.= 'U': upper triangular.= 'L': lower triangular.DIAGβCHARACTER*1. Specifies whether or not the matrix A is unit triangular.= 'N': non-unit triangular.= 'U': unit triangular.NβINTEGER. The order of the matrix A.N >= 0.AβDOUBLE PRECISION array, dimension(LDA,N). The triangular matrix A. IfUPLO = 'U', the leadingnbynupper triangular part of the arrayAcontains the upper triangular matrix, and the strictly lower triangular part ofAis not referenced. IfUPLO = 'L', the leadingnbynlower 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 be1.LDAβINTEGER. The leading dimension of the arrayA.LDA >= max(1,N).AINVβDOUBLE PRECISION array, dimension(LDAINV,N). On entry, the (triangular) inverse of the matrix A, in the same storage format as A. On exit, the contents ofAINVare destroyed. π So don't expect your input to survive, bestie.LDAINVβINTEGER. The leading dimension of the arrayAINV.LDAINV >= max(1,N).RCONDβDOUBLE PRECISION. The reciprocal condition number of A, computed as1/(norm(A) * norm(AINV)).WORKβDOUBLE PRECISION array, dimension(N).RESIDβDOUBLE PRECISION.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 122 of file dtrt01.f.
π Author π
This page was generated automatically by Doxygen for LAPACK from the source code.
Cute sign-off: LAPACK β Version 3.12.0 β TESTING/LIN/dtrt01.f(3). That's the footer vibe, bestie; slay and bye! πβ¨
π¦ 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