fiat team mailing list archive

[Johan Jansson <johanjan@xxxxxxxxxxxxxxxx>]

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