Hi Rodrigo,

Are you using the "f.at(...)" syntax, per https://www.firedrakeproject.org/point-evaluation.html ?

If so, you can pass in an extra argument tolerance, such as tolerance=1e-8

Under the hood, there are some calculations to see if the point is in a particular cell.  If the point requested is on the boundary of the domain, or on the boundary between cells, due to numerical precision issues, it may happen that all nearby cells think the point is not in that cell.  There is a default tolerance (of, I believe, 1e-14) to allow for this, but sometimes this is not enough.

If increasing this to 1e-8 doesn't work, there's probably something else going wrong!

Cheers,

Andrew

On Sat, 13 Jun 2020 at 20:22, Rodrigo Piccinini <rbpiccinini@gmail.com> wrote:

Hello,

I'm getting the PointNotInDomainError when trying to compute a point evaluation in a extruded mesh.

I've tried different coordinates that are inside a cell, but I keep having the same problem.

I'm using a DG scheme of 0 degree, so that values are constant inside a cell.

Is there a workaround to get a cell value?

Currently, I'm getting it from the vtk exported file, but it takes much more time.

Thanks in advance for your help!

--
Rodrigo