def Mat.gemv

bend-blas-lapack@0.0.0.1 · mat.bend:167

import Base
import bend-blas-lapack@0.0.0.1/blas.bend as B

167# y = alpha A x + beta y
168def Mat.gemv(+rows: Nat, +cols: Nat, alpha: F32, a: Mat<rows, cols>, x: Vec<cols>, beta: F32, y: Vec<rows>) -> IO(Result<&1, &1, U32 & String, Vec<rows>>):
169  MatA{m} = a
170  VecA{xs} = x
171  VecA{ys} = y
172  IO.bind(Result<&1, &1, U32 & String, Array<F32>>, Result<&1, &1, U32 & String, Vec<rows>>, B.Blas.sgemv(0, U32.from_nat(rows), U32.from_nat(cols), alpha, m, xs, beta, ys), r => Vec.wrap(rows, r))