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
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
