FFC

Comment 6 for bug 488661

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