TESTING/LIN/zhet01_aa.f(3)

okay bestie, this page is TESTING/LIN/zhet01_aa.f(3) from the Library Functions Manual 💖 buckle up, it’s giving numerical linear algebra with main character energy.

🌸 NAME 🌸

TESTING/LIN/zhet01_aa.f — that’s the whole name, no cap ✨

✨ SYNOPSIS ✨

🦋 Functions/Subroutines 🦋

subroutine zhet01_aa (uplo, n, a, lda, afac, ldafac, ipiv, c, ldc,
    rwork, resid)
    ZHET01_AA

💅 Function/Subroutine Documentation 💅

🫶 subroutine zhet01_aa

subroutine zhet01_aa (character uplo, integer n, complex*16, dimension( lda,
    * ) a, integer lda, complex*16, dimension( ldafac, * ) afac, integer
    ldafac, integer, dimension( * ) ipiv, complex*16, dimension( ldc, * ) c,
    integer ldc, double precision, dimension( * ) rwork, double precision
    resid)
    ZHET01_AA

ZHET01_AA 👑

🎀 Purpose 🎀

Okay so ZHET01_AA reconstructs a hermitian indefinite matrix A from its block L*D*L' or U*D*U' factorization, and then computes the residual:

norm( C - A ) / ( N * norm(A) * EPS ),

where C is the reconstructed matrix and EPS is the machine epsilon. Literally the whole glow-up: rebuild A, check the residual. 💖

💖 Parameters 💖

• UPLO — CHARACTER*1. Specifies whether the upper or lower triangular part of the hermitian matrix A is stored:
  • = 'U': Upper triangular
  • = 'L': Lower triangular
• N — INTEGER. The number of rows and columns of the matrix A. N >= 0.
• A — COMPLEX*16 array, dimension (LDA,N). The original hermitian matrix A.
• LDA — INTEGER. The leading dimension of the array A. LDA >= max(1,N)
• AFAC — COMPLEX*16 array, dimension (LDAFAC,N). The factored form of the matrix A. AFAC contains the block diagonal matrix D and the multipliers used to obtain the factor L or U from the block L*D*L' or U*D*U' factorization as computed by ZHETRF.
• LDAFAC — INTEGER. The leading dimension of the array AFAC. LDAFAC >= max(1,N).
• IPIV — INTEGER array, dimension (N). The pivot indices from ZHETRF.
• C — COMPLEX*16 array, dimension (LDC,N).
• LDC — INTEGER. The leading dimension of the array C. LDC >= max(1,N).
• RWORK — COMPLEX*16 array, dimension (N).
• RESID — COMPLEX*16.
  • If UPLO = 'L', norm(L*D*L' - A) / ( N * norm(A) * EPS )
  • If UPLO = 'U', norm(U*D*U' - A) / ( N * norm(A) * EPS )

👑 Author 👑

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

Definition at line 122 of file zhet01_aa.f. 📌

---

🌈 Author 🌈

Generated automatically by Doxygen for LAPACK from the source code. 🤖✨

---

okay bye bestie, this little icon is from LAPACK Version 3.12.0, file TESTING/LIN/zhet01_aa.f(3) — slay and stay hermitian 💜🎀