Not handling matrix dimensions correctly in gemm calls.
raineszm opened this issue · 0 comments
raineszm commented
When calling gemm, the arguments M, N, and K must describe the dimensions of the transposed matrices (as determined by the transa and transb) instead of the original matrices (http://www.netlib.org/lapack/explore-html/d7/d2b/dgemm_8f.html).
As far as I can tell, gemm is the only routine where this matters.