SSPT01(3)
π SSPT01(3) β Library Functions Manual β¨
okay bestie, this is SSPT01(3) from the Library Functions Manual, and it's giving "LAPACK test routine with main character energy." No cap, every detail is still here β we're just making it cute. π
πΈ NAME
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/LIN/sspt01.f
That's the file, periodt. It's the SSPT01 routine page, and yes, the path is long enough to need its own emotional support water bottle. π
β¨ SYNOPSIS
π Functions/Subroutines
Functions/Subroutines
subroutine SSPT01 (uplo, n, a, afac, ipiv, c, ldc, rwork, resid)
SSPT01
π§ Function/Subroutine Documentation
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
This is the "reconstruct and check" moment. SSPT01 reconstructs a symmetric indefinite packed matrix A from its block L*D*L' or U*D*U' factorization and computes the residual. The residual formula is exactly:
norm( C - A ) / ( N * norm(A) * EPS ),
where C is the reconstructed matrix and EPS is the machine epsilon. So it's literally comparing the reconstructed glow-up to the original and checking how close it is. π
π Parameters
UPLOβCHARACTER*1. Specifies whether the upper or lower triangular part of the symmetric matrixAis stored:= 'U': upper triangular;= 'L': lower triangular.NβINTEGER. The number of rows and columns of the matrixA.N >= 0.AβREAL array, dimension (N*(N+1)/2). The original symmetric matrixA, stored as a packed triangular matrix.AFACβREAL array, dimension (N*(N+1)/2). The factored form of the matrixA, stored as a packed triangular matrix.AFACcontains the block diagonal matrixDand the multipliers used to obtain the factorLorUfrom the blockL*D*L'orU*D*U'factorization as computed bySSPTRF.IPIVβINTEGER array, dimension (N). The pivot indices fromSSPTRF.CβREAL array, dimension (LDC,N).LDCβINTEGER. The leading dimension of the arrayC.LDC >= max(1,N).RWORKβREAL array, dimension (N).RESIDβREAL. IfUPLO = 'L',norm(L*D*L' - A) / ( N * norm(A) * EPS ); ifUPLO = '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.
π cute sign-off: this whole icon is from LAPACK, Version 3.12.0, and the footer path situation is /home/abui...IN/sspt01.f(3). mwah! π
π¦ where this page comes from
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms