search

JuliaLinearAlgebra / Octavian.jl

Multi-threaded BLAS-like library that provides pure Julia matrix multiplication
https://julialinearalgebra.github.io/Octavian.jl/stable/
Other
230 stars 18 forks source link

matmul! does not size-check #109

Open jaakkor2 opened 3 years ago

jaakkor2 commented 3 years ago 
using Octavian
n,m=100,10
y, A, x = rand(n), rand(n,m), rand(m+1)
matmul!(y,A,x)

whereas

using LinearAlgebra
n,m=100,10
y, A, x = rand(n), rand(n,m), rand(m+1)
mul!(y,A,x)

errors ERROR: DimensionMismatch("second dimension of A, 10, does not match length of x, 11")

Octavian v0.3.1