← Back to team overview

fiat team mailing list archive

Re: Bug in DiscontinuousVectorLagrange

 

Will patch this afternoon -- the bug came while I was in California for a month this summer...
Rob Kirby

"Mathematical software should be mathematical."



On Nov 30, 2005, at 8:40 AM, Anders Logg wrote:

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

_______________________________________________
FIAT-dev mailing list
FIAT-dev@xxxxxxxxxx
http://www.fenics.org/cgi-bin/mailman/listinfo/fiat-dev




References