ffc team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Mon, Apr 12, 2010 at 09:34:38PM +0800, Garth N. Wells wrote:
>
>
> On 12/04/10 21:29, Anders Logg wrote:
> >On Mon, Apr 12, 2010 at 09:21:32PM +0800, Garth N. Wells wrote:
> >>
> >>
> >>On 12/04/10 21:19, Garth N. Wells wrote:
> >>>
> >>>
> >>>On 12/04/10 20:47, Anders Logg wrote:
> >>>>We are doing some work where we need to do run-time quadrature over
> >>>>arbitrary polyhedra.
> >>>
> >>>We (Mehdi and I) do this already (using UFC), so I don't see why a new
> >>>function is required. Can you explain why evaluate_tensor is not enough?
> >>>
> >>
> >>I meant 'tabulate_tensor'.
> >
> >Which function do you call for evaluating the integrand?
> >
>
> We evaluate it inside ufc::tabulate_tensor. We construct our forms
> with an extra argument, say an object "CutCellIntegrator", which can
> provide quadrature schemes which depend on the considered cell.

That would require a special purpose code generator. Having
evaluate_integrand would allow more flexibility for users to implement
their own special quadrature scheme.

I also don't understand how it can work at all since the quadrature
points are not known at compile-time. Or do you a fixed set of reference
polyhedra that you map to?

--
Anders

Attachment: signature.asc
Description: Digital signature