dolfin team mailing list archive

[Kristian Oelgaard <k.b.oelgaard@xxxxxxxxx>]

On Wed, Nov 25, 2009 at 10:56 PM, Anders Logg <logg@xxxxxxxxx> wrote:

> I think I've tracked down the bug but don't know how it should be
> handled. We need input from Kristian on this.
>
> I've added some debug printing in FFC that should make things clear
> when trying to compile the following simple form with FFC:
>
>  P1 = FiniteElement("Lagrange", triangle, 1)
>  Q1 = FiniteElement("Quadrature", triangle, 2)
>
>  v = TestFunction(P1)
>  g = Function(Q1)
>
>  L = v*g*ds
>
> What happens is that FFC creates quadrature elements in compiler.py by
> calling create_element. The constructor of QuadratureElement then
> calls make_quadrature in FIAT just to get the number of quadrature
> points. In this call, it uses the cell, not the facet although it
> should be integrated on the facet.
>
> Then later make_quadrature is also called from quadraturerepresentation.py
> (strange that this is done twhice), but that time with the facet
> argument.
>
> The results is that in the first call to make_quadrature, the number
> of points is 4, but in the second case only 2.
>
> Kristian, could you take a quick look?
>

Yes, I'll take a look.

Kristian


> --
> Anders
>
> -----BEGIN PGP SIGNATURE-----
> Version: GnuPG v1.4.9 (GNU/Linux)
>
> iEYEARECAAYFAksNqAwACgkQTuwUCDsYZdF5gQCeOAqEqfHnGnLDU1VCWpashePN
> DuIAnilry+F1lXpgdKfTG2vQTDWKZDze
> =XmsH
> -----END PGP SIGNATURE-----
>
> _______________________________________________
> Mailing list: https://launchpad.net/~dolfin<https://launchpad.net/%7Edolfin>
> Post to     : dolfin@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~dolfin<https://launchpad.net/%7Edolfin>
> More help   : https://help.launchpad.net/ListHelp
>
>