From: Ian McCulloch (ianmcc_at_[hidden])
Date: 2007-07-02 18:59:52
Jeremy Bruestle wrote:
> I have developed a C++ tensor library using template expressions which
> allows Einstein summation convention to be used to describe
> mathematical operations on tensors in C++. This allows a great
> simplification in syntax, while retaining execution speeds identical
> to hand optimized C code in most cases. An example of the syntax is:
I might be interested in this, but the application I'm thinking of would
require primarily that it is *fast*, especially for large data sets where
attention to the cache heirachy is important. ie., it would need to use
BLAS as a back-end, while (implicitly) reordering indices in memory as
appropriate (perhaps the FLAME library would be useful here?
http://www.cs.utexas.edu/users/flame/). This might be a bit beyond what is
realistically possible for a generic tensor library, but I thought I'd ask.