TESTING/LIN/csyt03.f(3)

✨ NAME ✨

okay so this is the TESTING/LIN/csyt03.f(3) page from the Library Functions Manual, bestie 💖

🌸 SYNOPSIS 🌸

It’s giving “functions and subroutines” energy:

   Functions/Subroutines
       subroutine csyt03 (uplo, n, a, lda, ainv, ldainv, work, ldwork, rwork,
           rcond, resid)
           CSYT03

🎀 Function/Subroutine Documentation 🎀

Here’s the full signature, verbatim bestie:

subroutine csyt03 (character uplo, integer n, complex, dimension( lda, * )
    a, integer lda, complex, dimension( ldainv, * ) ainv, integer ldainv,
    complex, dimension( ldwork, * ) work, integer ldwork, real, dimension( *
    ) rwork, real rcond, real resid)

CSYT03

💅 Purpose 💅

Bestie, CSYT03 computes the residual for a complex symmetric matrix times its inverse, like this:

norm( I - A*AINV ) / ( N * norm(A) * norm(AINV) * EPS )

where EPS is the machine epsilon. slay ✨

🧚 Parameters 🧚

• UPLO — UPLO is CHARACTER*1. It specifies whether the upper or lower triangular part of the complex symmetric matrix A is stored: = 'U': Upper triangular; = 'L': Lower triangular. iykyk 💅
• N — N is INTEGER. The number of rows and columns of the matrix A. N >= 0.
• A — A is COMPLEX array, dimension (LDA,N). The original complex symmetric matrix A.
• 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 inverse of the matrix A, stored as a symmetric matrix in the same format as A. In this version, AINV is expanded into a full matrix and multiplied by A, so the opposing triangle of AINV will be changed; i.e., if the upper triangular part of AINV is stored, the lower triangular part will be used as work space. not me reading that twice 😭
• LDAINV — LDAINV is INTEGER. The leading dimension of the array AINV. LDAINV >= max(1,N).
• WORK — WORK is COMPLEX array, dimension (LDWORK,N).
• LDWORK — LDWORK is INTEGER. The leading dimension of the array WORK. LDWORK >= max(1,N).
• RWORK — RWORK is REAL array, dimension (N).
• RCOND — RCOND is REAL. The reciprocal of the condition number of A, computed as RCOND = 1/ (norm(A) * norm(AINV)).
• RESID — RESID is REAL. norm(I - A*AINV) / ( N * norm(A) * norm(AINV) * EPS ).

👑 Author 👑

• Univ. of Tennessee
• Univ. of California Berkeley
• Univ. of Colorado Denver
• NAG Ltd.

Definition at line 124 of file csyt03.f.

💜 Author 💜

Generated automatically by Doxygen for LAPACK from the source code.

---

okay bestie that’s the whole page — it’s LAPACK Version 3.12.0, footer path TESTING/LIN/csyt03.f(3), and that’s the vibe 💖✨🫶