← Back to team overview

fiat team mailing list archive

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