Alexander Schmolck schrieb:
Fernando Perez Fernando.Perez@colorado.edu writes:
I found something today a bit unpleasant: if you install numeric without any BLAS support, 'matrixmultiply is dot==True', so they are fully interchangeable. However, to my surprise, if you build numeric with the blas optimizations, they are NOT identical.
Oops, my bad (I submitted the patch and while pretty much all the real coding was done by Richard Everson this is my oversight).
No prob. It's been fixed in Numeric 23.5, so no more worries.
Pretty significant difference...
Yup, someone should incorporate optional atlas dot support into numarray if it hasn't happened already (won't be me, IIRC it took some convincing to get this into Numeric and I won't be using numarray for anything real in the near future).
I'll leave that question to the numarray guys, I have no idea where it stands in terms of blas/atlas support. I certainly hope it has it or that this optimization can be brought in, as it makes a huge difference for the large array case.