search

scijs / ndarray-gemm

General matrix multiplication for ndarrays
MIT License
18 stars 2 forks source link

Out-of-place matrix-matrix multiplication only? #2

Open fasiha opened 8 years ago

fasiha commented 8 years ago

Am I right in concluding, from some testing, that this works only for out-of-place operations, i.e., gemm(res, a, b) where a≠res≠b?

When I try something like gemm(c, c, a), the result is full of zeros.

If this is the case, it might be worth mentioning on the README. Are there any workarounds to enabling in-place multiplication, or would that ruin performance?

fasiha commented 8 years ago

Might be worth noting a simple implementation of dot over at scijs/ndarray#26 that creates intermediate ndarrays for each product it has to compute.