ffc team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Tue, Apr 13, 2010 at 08:27:11PM +0800, 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. For this reason, we would like to generate code
> >for evaluating the integrand of a form at an arbitrary point within
> >the cell.
> >
> >I would propose adding something like this to each of the *_integral
> >classes in the UFC interface:
> >
> >   /// Evaluate integrand at given point
> >   virtual void evaluate_integrand(double* A,
> >                                   const double * const * w,
> >                                   const cell&  c,
> >                                   const double* coordinates) const = 0;
> >
>
> What about
>
>     virtual void tabulate_tensor(double* A,
>                          const double * const * w,
>                          const cell&  c,
>                          uint num_quad_points,
>                          const double * const * coordinates,
>                          const double * weights) const = 0;
>
> instead? This would allow more of the quadrature code optimisations
> to be employed.

I thought about this but wasn't sure it was general enough for our
needs. I'll need to think about it.

--
Anders

Attachment: signature.asc
Description: Digital signature