dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Fri, Mar 06, 2009 at 03:40:46PM +0000, A Navaei wrote:
> 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?

There is no standard way to create higher ranked tensors, but the
GenericTensor interface (and the assembler) supports it so if you
should happen to have a data structure for rank 10 tensors, you can
assemble forms of arity 10 into it.

-- 
Anders

Attachment: signature.asc
Description: Digital signature