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
signature.asc
Description: Digital signature
_______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
