dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Thu, Nov 26, 2009 at 03:43:51PM +0100, Kristian Oelgaard wrote:
>
>
> 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

Great!

--
Anders

Attachment: signature.asc
Description: Digital signature