Re: [Libmesh-users] problem with inverse_map()

[Stogner, R. H <roy...@ic...>]

On Wed, 1 May 2019, Manav Bhatia wrote:

>   I am using h-refinement in my analysis, which uses the mesh function routines to compute the value of the function in the interior of an element.
>
>   All of my elements in the original mesh (before any refinements) are squares (quad4).
>
>   For the most part everything works out fine without any issues.
>   Occasionally, however, I will get an error in the inverse_map()
>   like this. I am particularly perplexed by the hmin() size of
>   10^-11.
>
>   The size of my elements before refinement is hmin() = .015 and I
>   allow a max of 4 refinements in any element. Would there be any
>   reason to expect an hmin of order 10^-11 in this case?

Not even close, but there's definitely *something* going seriously
wrong here.

You have a degenerate element; points 0 and 3 and points 1 and 2
coincide.

You have three nodes with invalid processor ids.

You probably ought to

> Rerun in devel/dbg mode for more details.

(preferably dbg) and see whether it catches any problems earlier.
---
Roy