Performs one of the matrix-vector operations
y := alpha*A*x + beta*y, or y := alpha*AT*x + beta*y, or
y := alpha*AH*x + beta*y,
where alpha and beta are scalars, x and y are vectors and A is an
m by n band matrix, with kl sub-diagonals and ku super-diagonals.