CTRT01(3)

πŸ“¦ version 3.12.0, in opensuse-leap 16.0

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 β€” UPLO is CHARACTER*1. Specifies whether the matrix A is upper or lower triangular.
    • = 'U': Upper triangular
    • = 'L': Lower triangular
  • DIAG β€” DIAG is CHARACTER*1. Specifies whether or not the matrix A is unit triangular.
    • = 'N': Non-unit triangular
    • = 'U': Unit triangular
  • N β€” N is INTEGER. The order of the matrix A. N >= 0.
  • A β€” A is COMPLEX array, dimension (LDA,N). The triangular matrix A. If UPLO = '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. If UPLO = '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. If DIAG = 'U', the diagonal elements of A are also not referenced and are assumed to be 1.
  • LDA β€” LDA is INTEGER. The leading dimension of the array A. LDA >= max(1,N).
  • AINV β€” AINV is COMPLEX 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 β€” LDAINV is INTEGER. The leading dimension of the array AINV. LDAINV >= max(1,N).
  • RCOND β€” RCOND is REAL. The reciprocal condition number of A, computed as 1/(norm(A) * norm(AINV)).
  • RWORK β€” RWORK is REAL array, dimension (N).
  • RESID β€” RESID is REAL. 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