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 _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
