In the case of 1/x[0], I'd use the approach that Martin suggested. It is much more efficient to define functions of coordinates in the UFL form directly. But in general, you can interpolate the function on a space of arbitrary degree, it's up to you. If you want the function to be evaluated directly at the same points where the integration takes place, use QuadratureElements as Garth suggested. Note that these are not supported on facets.
Kristian ---------- Forwarded message ---------- From: Ido Akkerman <[email protected]> Date: 3 June 2013 10:52 Subject: Re: [Bitbucket] Issue #49: Interior/Exterior facet coordinate bug? (fenics-project/dolfin) To: [email protected] ** [image: iakkerman] *Ido Akkerman* commented on issue #49: Interior/Exterior facet coordinate bug?<https://bitbucket.org/fenics-project/dolfin/issue/49/interior-exterior-facet-coordinate-bug> But that will only work because my example function is quadratic, right? What if I have 1/x[0]? Or something more elaborate? View this issue<https://bitbucket.org/fenics-project/dolfin/issue/49/interior-exterior-facet-coordinate-bug>or add a comment by replying to this email. Unsubscribe from issue emails<https://bitbucket.org/fenics-project/dolfin/issue/49/unsubscribe/k.b.oelgaard/3b47a71460570225ef2648df172a47cd1deca23c/>for this repository. [image: Bitbucket] <https://bitbucket.org>
_______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
