Improving dot product for standard multidimensional D arrays

[jmh530]

On Monday, 2 March 2020 at 13:35:15 UTC, p.shkadzko wrote:
> [snip]

Thanks. I don't have time right now to review this thoroughly. My 
recollection is that the dot product of two matrices is actually 
matrix multiplication, correct? It generally makes sense to defer 
to other people's implementation of this. I recommend trying 
lubeck's version against numpy. It uses a blas/lapack 
implementation. mir-glas, I believe, also has a version.

Also, I'm not sure if the fastmath attribute would do anything 
here, but something worth looking into.