Yeah, I think that's where most of the performance gains came from. I also added a re-write rule for unsafeGet dot (since it doesn't matter if the arguments are conjugated or not if the vectors are real) that shaved off about a tenth of a second.

Patrick

On Jul 24, 2008, at 4:26 PM, Don Stewart wrote:

patperry:
Last month Anatoly Yakovenko published some disturbing numbers about
the Haskell BLAS bindings I wrote being significantly slower than
using plain C.  I wanted to let everyone know that I've closed the
performance gap, and now for doing ten million dot products, the
overhead for using Haskell instead of C is about 0.6 seconds on my
machine, regardless of the size of the vectors. The next version will
incorporate the changes.  If you can't wait for a formal release, the
darcs repository is at http://www-stat.stanford.edu/~patperry/code/blas/

Anyone interested in more details can check out my blog:
http://quantile95.com/2008/07/24/addressing-haskell-blas-performance-issues/

Thanks everyone for the input on this (especially Anatoly).  If any
else finds any performance discrepancies, please let me know and I
will do whatever I can to fix them.


Great work, Patrick!

So if I read correctly, the main change was to flatten the
representation (and thus in loops the vector's structure will be
unpacked and kept in registers, which isn't possible for sum types).

-- Don

_______________________________________________
Haskell-Cafe mailing list
Haskell-Cafe@haskell.org
http://www.haskell.org/mailman/listinfo/haskell-cafe

Reply via email to