TESTING/LIN/sspt01.f(3)

πŸ“¦
version 3.12.1 Β· 3 versions
, in arch rolling

πŸŽ€ NAME πŸŽ€

okay bestie, this is TESTING/LIN/sspt01.f β€” straight from the Library Functions Manual, page TESTING/LIN/sspt01.f(3). It’s giving β€œLAPACK testing routine,” and the main character is SSPT01 in all its packed-matrix-reconstruction glory. πŸ’…

✨ SYNOPSIS ✨

πŸ¦‹ Functions/Subroutines πŸ¦‹

subroutine sspt01 (uplo, n, a, afac, ipiv, c, ldc, rwork, resid)

SSPT01

🫢 Function/Subroutine Documentation 🫢

Here’s the full signature, no cap:

subroutine sspt01 (character uplo, integer n, real, dimension( * ) a, real,
     dimension( * ) afac, integer, dimension( * ) ipiv, real, dimension( ldc,
     * ) c, integer ldc, real, dimension( * ) rwork, real resid)

SSPT01

✨ Purpose ✨

SSPT01 reconstructs a symmetric indefinite packed matrix A from its block L*D*L' or U*D*U' factorization and computes the residual, bestie:

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

where C is the reconstructed matrix and EPS is the machine epsilon. 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 symmetric 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 REAL array, dimension (N*(N+1)/2). The original symmetric matrix A, stored as a packed triangular matrix.

🌸 AFAC

AFAC is REAL array, dimension (N*(N+1)/2). The factored form of the matrix A, stored as a packed triangular matrix. 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 SSPTRF.

🌸 IPIV

IPIV is INTEGER array, dimension (N). The pivot indices from SSPTRF.

🌸 C

C is REAL array, dimension (LDC,N).

🌸 LDC

LDC is INTEGER. The leading dimension of the array C. LDC >= max(1,N).

🌸 RWORK

RWORK is REAL array, dimension (N).

🌸 RESID

RESID is REAL.

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 109 of file sspt01.f.

πŸ‘‘ Author πŸ‘‘

Generated automatically by Doxygen for LAPACK from the source code.


okay that’s the whole page, bestie πŸ’– β€” footer vibes are source/version <purple>LAPACK</purple>, <purple>Version 3.12.0</purple>, page <purple>TESTING/LIN/sspt01.f(3)</purple> β€” and the footer didn’t hand me a separate date, so no invented date era for us πŸ™ˆ. Mwah! 🫢✨

πŸ“¦ where this page comes from