SRC/sgetc2.f(3)

🌸 NAME 🌸

SRC/sgetc2.f β€” okay bestie, this is the page from the Library Functions Manual, and it's all about the sgetc2 subroutine. Iconic Unix/LAPACK behavior ahead. πŸ’…


✨ SYNOPSIS ✨

🧚 Functions/Subroutines 🧚

  • subroutine sgetc2 (n, a, lda, ipiv, jpiv, info) β€” SGETC2 computes the LU factorization with complete pivoting of the general n-by-n matrix. Slay.

πŸ’– Function/Subroutine Documentation πŸ’–

πŸ“ž subroutine sgetc2 (integer n, real, dimension( lda, * ) a, integer lda, integer, dimension( * ) ipiv, integer, dimension( * ) jpiv, integer info)

SGETC2 computes the LU factorization with complete pivoting of the general n-by-n matrix. Main character energy fr fr.

πŸ“ Purpose πŸ“

So the purpose block is basically this: SGETC2 computes an LU factorization with complete pivoting of the n-by-n matrix A. The factorization has the form A = P * L * U * Q, where P and Q are permutation matrices, L is lower triangular with unit diagonal elements, and U is upper triangular. Math girlies rise.

This is the Level 2 BLAS algorithm. Noted, bestie. 🫢

🌈 Parameters 🌈

  • N β€” N is INTEGER. It's the order of the matrix A. Constraint: N >= 0. No negative N, periodt.
  • A β€” A is REAL array, dimension (LDA, N). On entry, the n-by-n matrix A to be factored. On exit, the factors L and U from the factorization A = P*L*U*Q; the unit diagonal elements of L are not stored. If U(k, k) appears to be less than SMIN, U(k, k) is given the value of SMIN, i.e., giving a nonsingular perturbed system. Petition to store the unit diagonals? Denied, apparently.
  • LDA β€” LDA is INTEGER. It's the leading dimension of the array A. Constraint: LDA >= max(1,N). Respect the leading dimension or suffer.
  • IPIV β€” IPIV is INTEGER array, dimension(N). These are the pivot indices; for 1 <= i <= N, row i of the matrix has been interchanged with row IPIV(i). Row swaps, we love to see it.
  • JPIV β€” JPIV is INTEGER array, dimension(N). These are the pivot indices; for 1 <= j <= N, column j of the matrix has been interchanged with column JPIV(j). Column swaps too? Complete pivoting ate.
  • INFO β€” INFO is INTEGER.
    • = 0: successful exit. We love a clean exit.
    • > 0: if INFO = k, U(k, k) is likely to produce overflow if we try to solve for x in Ax = b. So U is perturbed to avoid the overflow. Overflow warning, bestie β€” do not ignore.

πŸ‘‘ Author πŸ‘‘

  • Univ. of Tennessee
  • Univ. of California Berkeley
  • Univ. of Colorado Denver
  • NAG Ltd.

🫢 Contributors: 🫢

Bo Kagstrom and Peter Poromaa, Department of Computing Science, Umea University, S-901 87 Umea, Sweden. Shoutout to the contributors, truly.

Definition at line 110 of file sgetc2.f.


πŸ‘‘ Author πŸ‘‘

Generated automatically by Doxygen for LAPACK from the source code. Automation ate, no cap.


okay bye bestie, this iconic little page is from LAPACK, Version 3.12.0, and its footer tag is SRC/sgetc2.f(3). πŸ’…βœ¨πŸŽ€

πŸ“¦ where this page comes from