fiat team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

It seems that the number of components for a discontinuous vector
Lagrange finite element is always set to the dimension of the shape
and not to the specified number of dimensions.

The output of the following code is

    (3,) (5,)

when it should be

    (5,) (5,)

    
from FIAT.DiscontinuousLagrange import DiscontinuousVectorLagrange
from FIAT.Lagrange import VectorLagrange
from FIAT.shapes import *

P0 = DiscontinuousVectorLagrange(TETRAHEDRON, 0, 5)
P1 = VectorLagrange(TETRAHEDRON, 1, 5)

print P0.function_space().tensor_dim()
print P1.function_space().tensor_dim()


/Anders