| Thread Previous • Date Previous • Date Next • Thread Next |
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.
Garth
I suspect/hope this should be easy to add to FFC since the current quadrature code must do something like this, but in addition iterate over the points and compute the weighted sum. This would let a user handle the quadrature manually in cases where that is necessary. Would it be possible to add? -- Anders _______________________________________________ Mailing list: https://launchpad.net/~ffc Post to : ffc@xxxxxxxxxxxxxxxxxxx Unsubscribe : https://launchpad.net/~ffc More help : https://help.launchpad.net/ListHelp
| Thread Previous • Date Previous • Date Next • Thread Next |
