Jeremy Bruestle wrote:
Greetings,
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. Regards, Ian McCulloch