dolfin team mailing list archive

[A Navaei <axnavaei@xxxxxxxxxxxxxx>]

Since the backend support in dolfin comes with vectors and matrices,
corresponding to tensors of rank 1 and 2, I was wondering if there is
a standard efficient way of creating higher ranked tensors based on
this.

For instance, if we are interested in creating sparse tensors of rank
4, at least in MTL style, should we create a 2D matrix of 2D matrix
type? That is:

mtl::Matrix< mtl::Matrix< double > >

I don't expect this to be as efficient as a 4D matrix as this cannot
be converted to a 1D array pointer. Am I missing something here?


-Ali