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