fiat team mailing list archive
-
fiat team
-
Mailing list archive
-
Message #00077
Re: Bug in DiscontinuousVectorLagrange
On Wed, Nov 30, 2005 at 02:55:28PM +0100, Johan Jansson wrote:
> 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
Aha. It's been so long I forgot about it.
Strange that this hasn't been in the news: "Critical security bug in
FIAT remains unpatched for months"... :-)
/Anders
Follow ups
References