On Tue, Nov 29, 2005 at 11:40:17PM -0600, Anders Logg wrote:
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
_______________________________________________
FIAT-dev mailing list
FIAT-dev@xxxxxxxxxx
http://www.fenics.org/cgi-bin/mailman/listinfo/fiat-dev
Hi!
Garth Wells sent a patch (to this list IIRC) which fixes this
bug. Here's his message:
http://www.fenics.org/pipermail/fiat-dev/2005-August/000060.html
and here's a direct link to the patch:
http://www.fenics.org/pipermail/fiat-dev/attachments/20050817/
ae24c432/P0_patch.obj
I have successfully used this patch for the updated elasticity forms
in DOLFIN, so it seems to work ok.
Please apply this patch ;).
Johan