ffc team mailing list archive

[Martin Sandve Alnæs <martinal@xxxxxxxxx>]

On Mon, Jun 22, 2009 at 9:10 PM, Kristian
Oelgaard<k.b.oelgaard@xxxxxxxxxx> wrote:
> Quoting Shawn Walker <walker@xxxxxxxxxxxxxxx>:
>
>> Actually, I am still a little confused.  Whenever you compute an integral
>> by quadrature, you need to evaluate the integrand at the quadrature
>> points.  Those are the points I need.
>
> If you have a quadrature element in your form
>
> element = FiniteElement("Quadrature", triangle, 2)
>
> you can call the
> tabulate_coordinates() function of the dof_map associated with the quadrature
> element. If you don't have a quadrature element there is no direct way of
> getting the quadrature points.
>
>> I just need the quad rule that is used in a particular case.  I see the
>> weights are defined in the tabulate_tensor routine.  Are the points listed
>> in a comment somewhere?
>
> No, currently not. But it might be a natural thing to have. I'll look into it.

It could be an idea to make it available through UFC,
but then again each ufc::*_integral may use different
quadrature rules for different terms...
However, an exception can be thrown if there is no unique rule in use.

Martin