TESTING/LIN/dpbt01.f(3)
β¨ NAME β¨
okay bestie, you're looking at TESTING/LIN/dpbt01.f(3) from the Library Functions Manual β and it's giving subroutine main character energy π
. This is the page for DPBT01.
π¦ SYNOPSIS π¦
Functions/Subroutines
subroutine dpbt01 (uplo, n, kd, a, lda, afac, ldafac, rwork, resid)
DPBT01
π Function/Subroutine Documentation π
subroutine dpbt01 (character uplo, integer n, integer kd, double precision,
dimension( lda, * ) a, integer lda, double precision, dimension( ldafac,
* ) afac, integer ldafac, double precision, dimension( * ) rwork, double
precision resid)
DPBT01
DPBT01 β yes, that's the name, and it's iconic. β¨
Purpose
Bestie, DPBT01 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. No cap, that's the whole assignment. π
Parameters
UPLOβ isCHARACTER*1. Specifies whether the upper or lower triangular part of the symmetric matrix A is stored:= 'U':Upper triangular= 'L':Lower triangular
Nβ isINTEGER. The number of rows and columns of the matrix A.N >= 0.KDβ isINTEGER. The number of super-diagonals of the matrix A ifUPLO = 'U', or the number of sub-diagonals ifUPLO = 'L'.KD >= 0.Aβ isDOUBLE PRECISIONarray, dimension(LDA,N). The original symmetric band matrix A. IfUPLO = 'U', the upper triangular part of A is stored as a band matrix; ifUPLO = 'L', the lower triangular part of A is stored. The columns of the appropriate triangle are stored in the columns of A and the diagonals of the triangle are stored in the rows of A. SeeDPBTRFfor further details.LDAβ isINTEGER. The leading dimension of the array A.LDA >= max(1,KD+1).AFACβ isDOUBLE PRECISIONarray, dimension(LDAFAC,N). The factored form of the matrix A. AFAC contains the factor L or U from the L*L' or U'*U factorization in band storage format, as computed byDPBTRF.LDAFACβ isINTEGER. The leading dimension of the array AFAC.LDAFAC >= max(1,KD+1).RWORKβ isDOUBLE PRECISIONarray, dimension(N).RESIDβ isDOUBLE PRECISION. 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 dpbt01.f.
π Author π
This page was generated automatically by Doxygen for LAPACK from the source code. Slay, automation. π€β¨
okay that's the whole page, bestie! Sending you off with the source/version tea: LAPACK Version 3.12.0, page TESTING/LIN/dpbt01.f(3) β bye, love you, mwah! πππ¦
π¦ 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