SRC/cla_lin_berr.f(3)
π NAME π
SRC/cla_lin_berr.f β bestie, this is the Fortran file that computes a component-wise relative backward error, and honestly? It's kind of the moment. β¨
π¦ SYNOPSIS π¦
Functions/Subroutines
subroutine cla_lin_berr (n, nz, nrhs, res, ayb, berr)
CLA_LIN_BERR computes a component-wise relative backward error. That's the whole vibe. π
πΈ Function/Subroutine Documentation πΈ
subroutine cla_lin_berr (integer n, integer nz, integer nrhs, complex,
dimension( n, nrhs ) res, real, dimension( n, nrhs ) ayb, real,
dimension( nrhs ) berr)
CLA_LIN_BERR computes a component-wise relative backward error. Say it with me now β component-wise relative backward error. Iconic. π
β¨ Purpose β¨
Okay so here's the actual formula, and yes you're gonna want to screenshot this:
!>
!> CLA_LIN_BERR computes componentwise relative backward error from
!> the formula
!> max(i) ( abs(R(i)) / ( abs(op(A_s))*abs(Y) + abs(B_s) )(i) )
!> where abs(Z) is the componentwise absolute value of the matrix
!> or vector Z.
!>
So abs(Z)? That's just the component-wise absolute value of the matrix or vector Z. Easy. No cap. π§
π«Ά Parameters π«Ά
NβN is INTEGER. The number of linear equations, i.e., the order of the matrixA. And it's gotta beN >= 0, periodt.NZβNZ is INTEGER. Lowkey the safety net of this whole operation: we add(NZ+1)*SLAMCH( 'Safe minimum' )toR(i)in the numerator to guard against spuriously zero residuals. Default value isN. πNRHSβNRHS is INTEGER. The number of right hand sides, i.e., the number of columns of the matricesAYB,RES, andBERR. Constraint:NRHS >= 0.RESβRES is COMPLEX array, dimension (N,NRHS). The residual matrix, i.e., the matrixRin the relative backward error formula above. Main character matrix, honestly. πAYBβAYB is REAL array, dimension (N, NRHS). The denominator in the relative backward error formula above, i.e., the matrixabs(op(A_s))*abs(Y) + abs(B_s). The matricesA,Y, andBare from iterative refinement (seecla_gerfsx_extended.f). Iykyk. π₯BERRβBERR is REAL array, dimension (NRHS). The componentwise relative backward error from the formula above. The final boss output. πβ¨
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line100of filecla_lin_berr.f.
π Author (again, but Doxygen this time) π
Generated automatically by Doxygen for LAPACK from the source code. Not me, not you β Doxygen did that. π
Okay bestie, that's the whole page, every last flag and formula intact, no details dropped, we ate and left zero crumbs. π
Stay safe out there, guard those spuriously zero residuals, and remember: only trust a backward error that's component-wise. π
β¨ LAPACK, Version 3.12.0 Β· SRC/cla_lin_berr.f(3) Β· mwah, love u, bye β¨ ππ¦π
π¦ 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