 <p>I just wanted to put that up on SO, but this pull request should be helpful and remove the need for a separate function for numpy.dot <a href="https://github.com/numpy/numpy/pull/2730" rel="nofollow">https://github.com/numpy/numpy/pull/2730</a> This should be available in numpy 1.7</p> <p>In the meantime, I used the example above to write a function that can replace numpy dot, whatever the order of your arrays are, and make the right call to fblas.dgemm. <a href="http://pastebin.com/M8TfbURi" rel="nofollow">http://pastebin.com/M8TfbURi</a></p> <p>Hope this helps,</p>