[vsipl++] [patch] General dispatch, use BLAS for dot- and matrix-matrix products
Mark Mitchell
mark at codesourcery.com
Tue Oct 11 20:16:32 UTC 2005
Jules Bergmann wrote:
> For dot-product, three implementations are provided: a generic
> implementation for all types and blocks, a BLAS implementation for BLAS
> types and direct data access blocks, and a BLAS implementation that
> specializes handling of a conjugated vector (for cvjdot). By using the
> block type to determine if conjugation is necessary, the expressions
> 'cvjdot(x, y)' and 'dot(x, conj(y))' are evaluated identically.
That's really, really cool. The performance numbers look great, too.
Does VSIPL have cvjdot? If not, we can probably show better performance
than any VSIPL implementation on this code.
--
Mark Mitchell
CodeSourcery, LLC
mark at codesourcery.com
(916) 791-8304
More information about the vsipl++
mailing list