Forgot to c/c the mailing list. ---------- Forwarded message --------- De: Rodrigo Piccinini <rbpiccinini@gmail.com> Date: sáb., 13 de jun. de 2020 às 19:40 Subject: Re: [firedrake] Point evaluation / PointNotInDomainError To: Andrew McRae <andrew.mcrae@physics.ox.ac.uk> Hi, Andrew. Thanks for the answer. I've attached a simplified example of what I'm doing. This reproduces the error I still get. [image: image.png] Em sáb., 13 de jun. de 2020 às 17:14, Andrew McRae < andrew.mcrae@physics.ox.ac.uk> escreveu:
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
-- -- Rodrigo -- -- Rodrigo
