def Mat.gemv_t

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

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

177# y = alpha A^T x + beta y
178def Mat.gemv_t(+rows: Nat, +cols: Nat, alpha: F32, a: Mat<rows, cols>, x: Vec<rows>, beta: F32, y: Vec<cols>) -> IO(Result<&1, &1, U32 & String, Vec<cols>>):
179  MatA{m} = a
180  VecA{xs} = x
181  VecA{ys} = y
182  IO.bind(Result<&1, &1, U32 & String, Array<F32>>, Result<&1, &1, U32 & String, Vec<cols>>, B.Blas.sgemv(1, U32.from_nat(rows), U32.from_nat(cols), alpha, m, xs, beta, ys), r => Vec.wrap(cols, r))