FFC

Comment 6 for bug 488661

Revision history for this message
Kristian B. Ølgaard (k.b.oelgaard) wrote : Re: [Bug 488661] Re: Compilation fails with quadrature element on boundary

OK, notice that the fix is not straightforward and we need to agree on
certain things to move on.
But since the bug has been inactive for 2-3 years I guess no one needs the fix.

On 3 January 2013 14:21, Marie Rognes <email address hidden> wrote:
> Since this has not been fixed in 2-3 years, I'm reducing the importance
> of this one.
>
> ** Changed in: ffc
> Milestone: 0.7.1 => None
>
> ** Changed in: ffc
> Importance: High => Low
>
> --
> You received this bug notification because you are a member of FFC Core
> Team, which is subscribed to FFC.
> https://bugs.launchpad.net/bugs/488661
>
> Title:
> Compilation fails with quadrature element on boundary
>
> Status in FEniCS Form Compiler:
> Confirmed
>
> Bug description:
> FFC fails to compile the following form:
>
> P1 = FiniteElement("Lagrange", triangle, 1)
> Q1 = FiniteElement("Quadrature", triangle, 2)
>
> v = TestFunction(P1)
> g = Function(Q1)
>
> L = v*g*ds
>
> The problem is that the constructor of QuadratureElement computes the number of quadrature points to be 4,
> but later quadraturerepresentation.py computes new points for the facet integral and then there are only 2
> points. As a result, the following test fails:
>
> if not len(points) == self.__num_quad_points:
> error("Points must be equal to coordinates of quadrature points")
>
> To manage notifications about this bug go to:
> https://bugs.launchpad.net/ffc/+bug/488661/+subscriptions