Shameless plug: if you need matrix multiplication, you may want to have a look
at my `linalg <https://github.com/unicredit/linear-algebra>`_ library, that
delegates to (various implementations of) BLAS or cuBLAS and provides a Nim DSL
on top of it, which includes dimension checking