On Mon, Aug 17, 2009 at 5:24 PM, Rutger ter Borg <rutger@terborg.net> wrote:

Jesse Perla wrote:

.. and it would be nice if operator* could be made a call to e.g. gemv (I
still want to check out the performance of CUDA).

My time for being able to work on the bindings have been a bit low the past
couple of months, I hope to be able to pick up where we left soon.

Very good idea. I find that people frequently want to have 2 rules for *, etc. depending on the matrix size. Perhaps this could be a compile time option in ublas? You set a minimum size threshold (#define UBLAS_LAPACK_THRESHOLD 200), and then a preprocessor instructions could say which implementation to use for the large matrices? Those implementations could be based on the bindings... Alternatively, this could be added into a policy inside of the matrix class itself, which would be more elegant, but tougher to write and could cause all sorts of regressions.

-Jesse