2008/3/14, Anders Logg <logg@xxxxxxxxx>:
On Thu, Mar 13, 2008 at 05:46:03PM -0400, Marie Rognes wrote:
'> > Just a comment: I typically want to write something like
>
> VectorElement("Nedelec", ...)
>
> in order to have a vector of Nedelec elements...!
So, does that represent a tensor field?
I've implemented value_rank()
in the VectorElement and TensorElement that way, so the value rank of
a VectorElement is the value rank of the subelement (Nedelec) plus
one. Thus each row of the tensor is represented by a Nedelec element.
If you want a set of vector elements instead:
fe = FiniteElement("Nedelec",...)
me = MixedElement(fe, fe, fe)
# or
me = fe+fe+fe
should also work. The value rank of a MixedElement constructed one of
these ways is undefined, since Functions should be defined over its
subelements anyway.