TESTING/LIN/cpot01.f(3)

📦
version 3.12.1 · 3 versions
, in arch rolling

💅 TESTING/LIN/cpot01.f(3) 💅

okay so this page is from the Library Functions Manual, and the star of the show is literally TESTING/LIN/cpot01.f. buckle up bestie, it's giving LAPACK testing routine energy ✨

🎀 NAME 🎀

TESTING/LIN/cpot01.f — no cap, that's the name 💖

✨ SYNOPSIS ✨

🧚 Functions/Subroutines 🧚

subroutine cpot01 (uplo, n, a, lda, afac, ldafac, rwork, resid)
    CPOT01

🦋 Function/Subroutine Documentation 🦋

okay here's the full signature, copied exactly because we do not play with technical details:

subroutine cpot01 (character uplo, integer n, complex, dimension( lda, * )
a, integer lda, complex, dimension( ldafac, * ) afac, integer ldafac,
real, dimension( * ) rwork, real resid)

CPOT01 — that's the routine, periodt 💅

🎯 Purpose 🎯

okay so CPOT01 reconstructs a Hermitian positive definite matrix A from its L*L' or U'*U factorization, and then it computes the residual:

norm( L*L' - A ) / ( N * norm(A) * EPS )
norm( U'*U - A ) / ( N * norm(A) * EPS )

where EPS is the machine epsilon, L' is the conjugate transpose of L, and U' is the conjugate transpose of U. literally that's the whole purpose, and it ate 🔥

🧸 Parameters 🧸

  • UPLO — UPLO is CHARACTER*1. It specifies whether the upper or lower triangular part of the Hermitian matrix A is stored:
    • = 'U': Upper triangular
    • = 'L': Lower triangular
  • 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 Hermitian matrix A.
  • LDA — LDA is INTEGER. The leading dimension of the array A. LDA >= max(1,N)
  • AFAC — AFAC is COMPLEX array, dimension (LDAFAC,N). On entry, the factor L or U from the L * L**H or U**H * U factorization of A. Overwritten with the reconstructed matrix, and then with the difference L * L**H - A (or U**H * U - A).
  • LDAFAC — LDAFAC is INTEGER. The leading dimension of the array AFAC. LDAFAC >= max(1,N).
  • RWORK — RWORK is REAL array, dimension (N)
  • RESID — RESID is REAL. If UPLO = 'L', norm(L * L**H - A) / ( N * norm(A) * EPS ); if UPLO = 'U', norm(U**H * U - A) / ( N * norm(A) * EPS )

🧑‍💻 Author 🧑‍💻

these iconic institutions/companies, fr fr:

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

Definition at line 105 of file cpot01.f. — yup, line 105, bestie 📌

🤖 Author 🤖

Generated automatically by Doxygen for LAPACK from the source code. — the robot did its thing 🫶

💖 okay that's the whole page — source/version: LAPACK, Version 3.12.0; page: TESTING/LIN/cpot01.f(3); and the footer doesn't list a date, so no fake date from me, bestie 😭🫶✨

📦 where this page comes from