← Back to team overview

ffc team mailing list archive

Re: 4th order tensors

 

Jake Ostien wrote:
Hi,

I'd like to construct a 4th order Identity tensor. I am having some trouble trying to figure out if/how FFC handles tensors of higher order than 2. If anyone has successfully accomplished this, and would like to provide some guidance, it would be appreciated. Currently I am trying to hard-code in 81 terms and I am getting lost in the indices.

As far as I know, numpy doesn't support higher order identity tensors. Even if it did, I would be careful using it since useful higher-order tensors typically possess many symmetries. If you don't take this into account, you might have to wait a long time for FFC to compute anything.

I know that supporting tensor-valued functions is on the TODO list.

Garth


A
similar question would be how numpy handles arrays with more than two indices?

Thanks
Jake
_______________________________________________
FFC-dev mailing list
FFC-dev@xxxxxxxxxx
http://www.fenics.org/mailman/listinfo/ffc-dev





Follow ups

References