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