lacn2(3)
π
bestie's guide to lacn2(3) π
Okay so this is lacn2(3) from the Library Functions Manual, and ngl it's a whole family situation β four subroutines, all doing the same iconic thing in different precisions. Let's get into it π«Ά
β¨ NAME β¨
lacn2 β it's the lacn2: a 1-norm estimate, e.g., || A^{-1} ||_1 in gecon. Periodt. π
πΈ SYNOPSIS πΈ
π Functions
subroutine clacn2 (n, v, x, est, kase, isave)β CLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vector products.subroutine dlacn2 (n, v, x, isgn, est, kase, isave)β DLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vector products.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.subroutine zlacn2 (n, v, x, est, kase, isave)β ZLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vector products.
π¦ Detailed Description π¦
Not me checking this section three times β bestie there is literally nothing here. It's a heading with no body. The actual tea is all in Function Documentation below, so keep scrolling π
π₯ Function Documentation π₯
π clacn2
subroutine clacn2 (integer n, complex, dimension( * ) v, complex, dimension( * ) x, real est, integer kase, integer, dimension( 3 ) isave)
CLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vector products.
β Purpose
CLACN2estimates the 1-norm of a square, complex matrixA. Reverse communication is used for evaluating matrix-vector products.
π Parameters
NβNisINTEGER. The order of the matrix.N >= 1. No zero-sized matrices bestie, the math said so.VβVisCOMPLEXarray, dimension (N). On the final return,V = A*W, whereEST = norm(V)/norm(W)(Wis not returned). Yeah she's not giving youWback, don't ask πXβXisCOMPLEXarray, dimension (N). On an intermediate return,Xshould be overwritten byA * XifKASE=1, orA**H * XifKASE=2, whereA**His the conjugate transpose ofA, andCLACN2must be re-called with all the other parameters unchanged. All of them. Don't touch. πESTβESTisREAL. On entry withKASE = 1or2andISAVE(1) = 3,ESTshould be unchanged from the previous call toCLACN2. On exit,ESTis an estimate (a lower bound) fornorm(A).KASEβKASEisINTEGER. On the initial call toCLACN2,KASEshould be0. On an intermediate return,KASEwill be1or2, indicating whetherXshould be overwritten byA * XorA**H * X. On the final return fromCLACN2,KASEwill again be0. It's giving β¨full circleβ¨ISAVEβISAVEisINTEGERarray, dimension (3).ISAVEis used to save variables between calls toSLACN2.
Yeah bestie, that last line saysSLACN2on theCLACN2page π that's just how the doc is written, I'm reporting it exactly as it lives.
π Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
π Further Details:
- Originally named
CONEST, dated March 16, 1988. - Last modified: April, 1999
- This is a thread safe version of
CLACON, which uses the arrayISAVEin place of aSAVEstatement, as follows:
CLACON CLACN2
JUMP ISAVE(1)
J ISAVE(2)
ITER ISAVE(3)
π«Ά Contributors:
Nick Higham, University of Manchester
π 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 132 of file clacn2.f.
π§ dlacn2
subroutine dlacn2 (integer n, double precision, dimension( * ) v, double precision, dimension( * ) x, integer, dimension( * ) isgn, double precision est, integer kase, integer, dimension( 3 ) isave)
DLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vector products.
β Purpose
DLACN2estimates the 1-norm of a square, real matrixA. Reverse communication is used for evaluating matrix-vector products.
π Parameters
NβNisINTEGER. The order of the matrix.N >= 1.VβVisDOUBLE PRECISIONarray, dimension (N). On the final return,V = A*W, whereEST = norm(V)/norm(W)(Wis not returned).XβXisDOUBLE PRECISIONarray, dimension (N). On an intermediate return,Xshould be overwritten byA * XifKASE=1, orA**T * XifKASE=2, andDLACN2must be re-called with all the other parameters unchanged.ISGNβISGNisINTEGERarray, dimension (N). That's it. That's the whole description. Iconic minimalism πESTβESTisDOUBLE PRECISION. On entry withKASE = 1or2andISAVE(1) = 3,ESTshould be unchanged from the previous call toDLACN2. On exit,ESTis an estimate (a lower bound) fornorm(A).KASEβKASEisINTEGER. On the initial call toDLACN2,KASEshould be0. On an intermediate return,KASEwill be1or2, indicating whetherXshould be overwritten byA * XorA**T * X. On the final return fromDLACN2,KASEwill again be0.ISAVEβISAVEisINTEGERarray, dimension (3).ISAVEis used to save variables between calls toDLACN2.
π Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
π Further Details:
- Originally named
SONEST, dated March 16, 1988. - This is a thread safe version of
DLACON, which uses the arrayISAVEin place of aSAVEstatement, as follows:
DLACON DLACN2
JUMP ISAVE(1)
J ISAVE(2)
ITER ISAVE(3)
π«Ά Contributors:
Nick Higham, University of Manchester
π 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 dlacn2.f.
π slacn2
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.
β Purpose
SLACN2estimates the 1-norm of a square, real matrixA. Reverse communication is used for evaluating matrix-vector products.
π Parameters
NβNisINTEGER. The order of the matrix.N >= 1.VβVisREALarray, dimension (N). On the final return,V = A*W, whereEST = norm(V)/norm(W)(Wis not returned).XβXisREALarray, dimension (N). On an intermediate return,Xshould be overwritten byA * XifKASE=1, orA**T * XifKASE=2, andSLACN2must be re-called with all the other parameters unchanged.ISGNβISGNisINTEGERarray, dimension (N). Again, that's the whole thing. She said "no notes" πESTβESTisREAL. On entry withKASE = 1or2andISAVE(1) = 3,ESTshould be unchanged from the previous call toSLACN2. On exit,ESTis an estimate (a lower bound) fornorm(A).KASEβKASEisINTEGER. 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.ISAVEβISAVEisINTEGERarray, dimension (3).ISAVEis used to save variables between calls toSLACN2.
π Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
π Further Details:
- Originally named
SONEST, dated March 16, 1988. - This is a thread safe version of
SLACON, which uses the arrayISAVEin place of aSAVEstatement, as follows:
SLACON SLACN2
JUMP ISAVE(1)
J ISAVE(2)
ITER ISAVE(3)
π«Ά Contributors:
Nick Higham, University of Manchester
π 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.
π¦ zlacn2
subroutine zlacn2 (integer n, complex*16, dimension( * ) v, complex*16, dimension( * ) x, double precision est, integer kase, integer, dimension( 3 ) isave)
ZLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vector products.
β Purpose
ZLACN2estimates the 1-norm of a square, complex matrixA. Reverse communication is used for evaluating matrix-vector products.
π Parameters
NβNisINTEGER. The order of the matrix.N >= 1.VβVisCOMPLEX*16array, dimension (N). On the final return,V = A*W, whereEST = norm(V)/norm(W)(Wis not returned).XβXisCOMPLEX*16array, dimension (N). On an intermediate return,Xshould be overwritten byA * XifKASE=1, orA**H * XifKASE=2, whereA**His the conjugate transpose ofA, andZLACN2must be re-called with all the other parameters unchanged.ESTβESTisDOUBLE PRECISION. On entry withKASE = 1or2andISAVE(1) = 3,ESTshould be unchanged from the previous call toZLACN2. On exit,ESTis an estimate (a lower bound) fornorm(A).KASEβ π okay bestie this is exactly where part 1 cuts off βKASE, plus the rest ofzlacn2's details, roll right into part 2. Don't move, we're mid-sentence on purpose πKASEβ okay soKASEisINTEGERπ . On the initial call toZLACN2,KASEshould be0. On an intermediate return,KASEwill be1or2, indicating whetherXshould be overwritten byA * XorA**H * X. On the final return fromZLACN2,KASEwill again be0. π«ΆISAVEβISAVEisINTEGER array, dimension (3)β¨.ISAVEis used to save variables between calls toZLACN2. Lowkey the memory keeper of the whole situation. π
πΈ Author πΈ
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π§ Further Details:
Originally named CONEST, dated March 16, 1988. β¨
Last modified: April, 1999 π°οΈ
This is a thread safe version of ZLACON, which uses the array ISAVE in place of a SAVE statement, as follows:
ZLACON ZLACN2
JUMP ISAVE(1)
J ISAVE(2)
ITER ISAVE(3)
π Contributors:
- Nick Higham, University of Manchester
π 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 132 of file zlacn2.f. π«Ά
π Author π
Generated automatically by Doxygen for LAPACK from the source code. π€π
Cute sign-off: this is LAPACK Version 3.12.0, page lacn2(3) β and the footer didnβt give us a date, so thatβs the whole story, bestie. π
β¨
π¦ 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