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
signature.asc
Description: Digital signature
_______________________________________________ Mailing list: https://launchpad.net/~ffc Post to : ffc@lists.launchpad.net Unsubscribe : https://launchpad.net/~ffc More help : https://help.launchpad.net/ListHelp
