SRC/slacn2.f(3)
π SRC/slacn2.f π
okay so welcome to the Library Functions Manual page for SRC/slacn2.f(3) β yes bestie, we're deep in LAPACK territory and this subroutine is lowkey iconic. It estimates the 1-norm of a square real matrix using reverse communication for matrix-vector products, and it's giving "you do the matvec, I'll handle the norm estimate" energy. π«Άβ¨
π NAME π
SRC/slacn2.f β that's the file, bestie. It's the home of slacn2, and the whole page is about that subroutine. π
π SYNOPSIS π
π¦ Functions/Subroutines π¦
subroutine slacn2 (n, v, x, isgn, est, kase, isave)
SLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vector products. No cap. π«
π₯ Function/Subroutine Documentation π₯
subroutine slacn2 (integer n, real, dimension( * ) v, real, dimension( * )
x, integer, dimension( * ) isgn, real est, integer kase, integer,
dimension( 3 ) isave)
SLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vector products. That's the whole vibe, and honestly it ate. π
π― Purpose π―
SLACN2 estimates the 1-norm of a square, real matrix A. Reverse communication is used for evaluating matrix-vector products. Periodt. π
πΈ Parameters πΈ
Nβ INTEGER. The order of the matrix.N >= 1. So yes, there's at least one row and one column, bestie. πVβ REAL array, dimension(N). On the final return,V = A*W, whereEST = norm(V)/norm(W)(Wis not returned). It's giving "I saw W but I'm not handing it over." πXβ REAL array, dimension(N). On an intermediate return,Xshould be overwritten by:A * X, ifKASE=1,A**T * X, ifKASE=2,
and SLACN2 must be re-called with all the other parameters unchanged. Do not freestyle with the other params, okay? π
ISGNβ INTEGER array, dimension(N). That's it, that's the whole note. πESTβ REAL. On entry withKASE = 1or2andISAVE(1) = 3,ESTshould be unchanged from the previous call toSLACN2. On exit,ESTis an estimate (a lower bound) fornorm(A). It's a lower bound, so don't @ me. β¨KASEβ INTEGER. On the initial call toSLACN2,KASEshould be0. On an intermediate return,KASEwill be1or2, indicating whetherXshould be overwritten byA * XorA**T * X. On the final return fromSLACN2,KASEwill again be0. The full circle moment. πISAVEβ INTEGER array, dimension(3).ISAVEis used to save variables between calls toSLACN2. It's the little memory keeper. π§
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π₯ Further Details π₯
Originally named SONEST, dated March 16, 1988. So yes, this icon has been serving since the 80s. π
This is a thread safe version of SLACON, which uses the array ISAVE in place of a SAVE statement, as follows:
SLACON SLACN2
JUMP ISAVE(1)
J ISAVE(2)
ITER ISAVE(3)
π Contributors π
Nick Higham, University of Manchester. We stan. π«Ά
π References π
N.J. Higham, "FORTRAN codes for estimating the one-norm of a real or complex matrix, with applications to condition estimation", ACM Trans. Math. Soft., vol. 14, no. 4, pp. 381-396, December 1988.
Definition at line 135 of file slacn2.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code. It's the meta author, bestie. π€π
okay bye bestie β this page is from LAPACK Version 3.12.0, and the file is SRC/slacn2.f(3). 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