zgbmv zgemv zgerc zgeru zhbmv zhemv zher zher2 zhpmv zhpr zhpr2 zsbmv zspmv zspr zspr2 zsymv zsyr zsyr2 ztbmv ztbsv ztpmv ztpsv ztrmv ztrsv

zgbmv() void zgbmv ( char  trans, int  m, int  n, int  kl, int  ku, doublecomplex  alpha, int  ldab, doublecomplex  ab[], doublecomplex  x[], int  incx, doublecomplex  beta, doublecomplex  y[], int  incy  ) y <- αAx + βy, y <- αATx + βy or y <- αAHx + βy (complex band matrices) (BLAS 2) Purpose This routine performs one of the matrix-vector operations y <- alpha*A*x + beta*y, y <- alpha*A^T*x + beta*y or y <- alpha*A^H*x + beta*y where alpha and beta are scalars, x and y are vectors and A is an m x n band matrix, with kl sub-diagonals and ku super-diagonals. Parameters [in] trans Specifies the operation to be performed as follows: = 'N': y <- alpha*A*x + beta*y. = 'T': y <- alpha*A^T*x + beta*y. = 'C': y <- alpha*A^H*x + beta*y. [in] m Number of rows of the matrix A. (m >= 0) (If m = 0, returns without computation) [in] n Number of columns of the matrix A. (n >= 0) (If n = 0, returns without computation) [in] kl Number of sub-diagonals of the matrix A. (kl >= 0) [in] ku Number of super-diagonals of the matrix A. (ku >= 0) [in] alpha Scalar alpha. [in] ldab Leading dimension of the two dimensional array ab[][]. (ldab >= kl + ku + 1) [in] ab[][] Array ab[lab][ldab] (lab >= n) m x n band matrix A in kl+ku+1 x n band matrix form. [in] x[] Array x[lx] (lx >= 1 + (n - 1)*abs(incx) if trans = 'N', lx >= 1 + (m - 1)*abs(incx) otherwise) Input vector x. [in] incx Storage spacing between elements of x. (incx != 0) [in] beta Scalar beta. When beta is supplied as zero then y[] need not be set on input. [in,out] y[] Array y[ly] (ly >= 1 + (m - 1)*abs(incy) if trans = 'N', ly >= 1 + (n - 1)*abs(incy) otherwise) [in] Input vector y. [out] Output vector. (= alpha*A*x + beta*y, alpha*A^T*x + beta*y or alpha*A^H*x + beta*y) [in] incy Storage spacing between elements of y. (incy != 0) Reference BLAS