dolfin team mailing list archive

["Garth N. Wells" <g.n.wells@xxxxxxxxxx>]

I'd like to resurrect some functionality for dealing with small dense
matrices (addition, multiplication, eigenvalues, inverse, etc.).
Versions of DOLFIN prior to the transition to PETSc had some of this.

I'm open to suggestions as to the best data structure for this. The code
in /src/kerenel/common/dolfin/Tensor.h could provide the underlying data
structure, although perhaps this code needs to be readdressed. I 'm
thinking we could

1. Keep Tensor.h as it is.

2. Derive Tensor.h from src/kernel/common/dolfin/Array.h, which in turn
is based on std::vector.

3. Make Tensor.h (and perhaps Array.h also) a wrapper for another
library that handles multi-dimensional arrays (for this Blitz++ looks
nice, http://www.oonumerics.org/blitz/, and there is also
http://www.boost.org/libs/multi_array/doc/index.html). Any suggestions?

Garth