BLAS/SRC/zgbmv.f(3)

✨ zgbmv β€” the band matrix girlie ✨

ok so this is BLAS/SRC/zgbmv.f, bestie πŸ’– she lives in the BLAS folder and she does band matrix-vector math and honestly? she ate.


πŸŽ€ SYNOPSIS πŸŽ€

Functions/Subroutines

subroutine zgbmv (trans, m, n, kl, ku, alpha, a, lda, x, incx, beta, y, incy)

ZGBMV πŸ‘‘


πŸ’œ Function/Subroutine Documentation πŸ’œ

subroutine zgbmv

subroutine zgbmv (character trans, integer m, integer n, integer kl, integer
    ku, complex*16 alpha, complex*16, dimension(lda,*) a, integer lda,
    complex*16, dimension(*) x, integer incx, complex*16 beta, complex*16,
    dimension(*) y, integer incy)

ZGBMV πŸ’…


🌸 Purpose 🌸

ok so ZGBMV performs one of these matrix-vector operations, pick your fighter:

y := alpha*A*x + beta*y,   or   y := alpha*A**T*x + beta*y,   or

y := alpha*A**H*x + beta*y,

where alpha and beta are scalars (just lil numbers, no cap), x and y are vectors, and A is an m by n band matrix with kl sub-diagonals and ku super-diagonals. iconic honestly πŸ¦‹


🌈 Parameters 🌈

  • TRANS β€” TRANS is CHARACTER*1 β€” on entry, TRANS tells you which operation is happening, and it's giving choose-your-own-adventure:
    • TRANS = 'N' or 'n' β†’ y := alpha*A*x + beta*y.
    • TRANS = 'T' or 't' β†’ y := alpha*A**T*x + beta*y.
    • TRANS = 'C' or 'c' β†’ y := alpha*A**H*x + beta*y.
  • M β€” M is INTEGER β€” on entry, M is the number of rows of the matrix A. M must be at least zero. (zero is allowed bestie, we don't judge πŸ’–)
  • N β€” N is INTEGER β€” on entry, N is the number of columns of the matrix A. N must be at least zero.
  • KL β€” KL is INTEGER β€” on entry, KL is the number of sub-diagonals of the matrix A. KL must satisfy 0 .le. KL. (yes the .le. is Fortran-speak for ≀, we love her 🌸)
  • KU β€” KU is INTEGER β€” on entry, KU is the number of super-diagonals of the matrix A. KU must satisfy 0 .le. KU.
  • ALPHA β€” ALPHA is COMPLEX*16 β€” on entry, ALPHA is the scalar alpha. main character scalar energy ✨
  • A β€” A is COMPLEX*16 array, dimension ( LDA, N ) β€” okay this one's a whole situation, pay attention:
    • before entry, the leading ( kl + ku + 1 ) by n part of the array A must contain the matrix of coefficients, supplied column by column, with:
      • the leading diagonal of the matrix in row ( ku + 1 ) of the array,
      • the first super-diagonal starting at position 2 in row ku,
      • the first sub-diagonal starting at position 1 in row ( ku + 2 ),
      • and so on. πŸ”
    • elements in the array A that don't correspond to elements in the band matrix (like the top left ku by ku triangle) are not referenced β€” they're just vibing, ignored, unbothered πŸ’…
    • the following program segment will transfer a band matrix from conventional full matrix storage to band storage:
      DO 20, J = 1, N
         K = KU + 1 - J
         DO 10, I = MAX( 1, J - KU ), MIN( M, J + KL )
            A( K + I, J ) = matrix( I, J )
   10    CONTINUE
   20 CONTINUE
  • LDA β€” LDA is INTEGER β€” on entry, LDA specifies the first dimension of A as declared in the calling (sub) program. LDA must be at least ( kl + ku + 1 ). no shortcuts bestie 🫢
  • X β€” X is COMPLEX*16 array, dimension at least ( 1 + ( n - 1 )*abs( INCX ) ) when TRANS = 'N' or 'n' and at least ( 1 + ( m - 1 )*abs( INCX ) ) otherwise. Before entry, the incremented array X must contain the vector x. πŸ¦‹
  • INCX β€” INCX is INTEGER β€” on entry, INCX specifies the increment for the elements of X. INCX must not be zero. 🚨 (fr fr, do not)
  • BETA β€” BETA is COMPLEX*16 β€” on entry, BETA specifies the scalar beta. When BETA is supplied as zero then Y need not be set on input. little treat for you πŸ“
  • Y β€” Y is COMPLEX*16 array, dimension at least ( 1 + ( m - 1 )*abs( INCY ) ) when TRANS = 'N' or 'n' and at least ( 1 + ( n - 1 )*abs( INCY ) ) otherwise. Before entry, the incremented array Y must contain the vector y. On exit, Y is overwritten by the updated vector y. If either m or n is zero, then Y not referenced and the function performs a quick return. she said "not today" and left 😭
  • INCY β€” INCY is INTEGER β€” on entry, INCY specifies the increment for the elements of Y. INCY must not be zero. πŸ’€

πŸ‘‘ Author πŸ‘‘

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

πŸ”₯ Further Details: πŸ”₯

  • Level 2 Blas routine.
  • The vector and matrix arguments are not referenced when N = 0, or M = 0
  • -- Written on 22-October-1986.
    • Jack Dongarra, Argonne National Lab.
    • Jeremy Du Croz, Nag Central Office.
    • Sven Hammarling, Nag Central Office.
    • Richard Hanson, Sandia National Labs.
Definition at line 188 of file zgbmv.f.

🧚 Author 🧚

Generated automatically by Doxygen for LAPACK from the source code. she's literally generated, she's literally that girl πŸ€–πŸ’–


ok bye bestie, that's the whole page πŸŽ€ this one came from Version 3.12.0, Tue Nov 28 2023 12:08:41, and the file is BLAS/SRC/zgbmv.f(3). slay, go call your band matrix, but like... with the right LDA. mwah πŸ’…βœ¨