TESTING/LIN/spbt01.f(3)
π©° TESTING/LIN/spbt01.f
okay bestie, welcome to the Library Functions Manual page for TESTING/LIN/spbt01.f β itβs giving super-specific LAPACK testing routine energy and Iβm lowkey obsessed. Letβs unpack it π
πΈ NAME
TESTING/LIN/spbt01.f
β¨ SYNOPSIS
π₯ Functions/Subroutines
subroutine spbt01 (uplo, n, kd, a, lda, afac, ldafac, rwork, resid)
SPBT01
π Function/Subroutine Documentation
subroutine spbt01 (character uplo, integer n, integer kd, real, dimension(
lda, * ) a, integer lda, real, dimension( ldafac, * ) afac, integer
ldafac, real, dimension( * ) rwork, real resid)
SPBT01
π Purpose
SPBT01 reconstructs a symmetric positive definite band matrix A from its L*L' or U'*U factorization and computes the residual:
norm( L*L' - A ) / ( N * norm(A) * EPS ) or
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. β¨
π Parameters
UPLOβUPLO is CHARACTER*1. It tells you whether the upper or lower triangular part of the symmetric matrixAis stored:= 'U': Upper triangular= 'L': Lower triangular
NβN is INTEGER. The number of rows and columns of the matrixA.N >= 0.KDβKD is INTEGER. The number of super-diagonals of the matrixAifUPLO = 'U', or the number of sub-diagonals ifUPLO = 'L'.KD >= 0.AβA is REAL array, dimension (LDA,N). The original symmetric band matrixA. IfUPLO = 'U', the upper triangular part ofAis stored as a band matrix; ifUPLO = 'L', the lower triangular part ofAis stored. The columns of the appropriate triangle are stored in the columns ofAand the diagonals of the triangle are stored in the rows ofA. SeeSPBTRFfor further details.LDAβLDA is INTEGER.The leading dimension of the arrayA.LDA >= max(1,KD+1).AFACβAFAC is REAL array, dimension (LDAFAC,N). The factored form of the matrixA.AFACcontains the factorLorUfrom theL*L'orU'*Ufactorization in band storage format, as computed bySPBTRF.LDAFACβLDAFAC is INTEGER. The leading dimension of the arrayAFAC.LDAFAC >= max(1,KD+1).RWORKβRWORK is REAL array, dimension (N).RESIDβRESID is REAL. IfUPLO = 'L',norm(L*L' - A) / ( N * norm(A) * EPS ). IfUPLO = 'U',norm(U'*U - A) / ( N * norm(A) * EPS ).
π Author
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 117 of file spbt01.f.
π§ Author
Generated automatically by Doxygen for LAPACK from the source code.
mwah π catch you in the source, bestie β this footer is serving LAPACK Version 3.12.0 and the path TESTING/LIN/spbt01.f(3) π«Άβ¨
π¦ where this page comes from
- arch rolling: lapack-doc 3.12.1-2 Β· custom
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms