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

SourceForge Headquarters 225 Broadway Suite 1600 San Diego, CA 92101 +1 (858) 422-6466

I ran the debug version and this is where it crashed. 

[2] src/fe/fe_map.C, line 1866, compiled Apr 29 2019 at 13:59:12
WARNING: Newton scheme has not converged in 11 iterations:
   physical_point=(x,y,z)=(    0.27,   0.2775,        0)   physical_guess=(x,y,z)=(    0.27,   0.2775,        0)   dp=(x,y,z)=(-1.21944e-13, -1.12163e-06,        0)   p=(x,y,z)=(      -1, -0.57735,        0)   error=1.12163e-06   in element 16514
  Elem Information
   id()=16514, processor_id()=2
   type()=QUAD4
   dim()=2
   n_nodes()=4
    0  Node id()=13234, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
    1  Node id()=604, processor_id()=2, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=(0/0/13470) (0/1/13471) (0/2/13472) (0/3/13473) (0/4/13474) (0/5/13475) (2/0/2245) (3/0/2245) 
    2  Node id()=13235, processor_id()=4294967295, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=
    3  Node id()=13233, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
   n_sides()=4
    neighbor(0)=nullptr
    neighbor(1)=nullptr
    neighbor(2)=nullptr
    neighbor(3)=nullptr
   hmin()=9.89829e-11, hmax()=0.001875
   volume()=1.85593e-13
   active()=1, ancestor()=0, subactive()=0, has_children()=0
   parent()=672
   level()=4, p_level()=0
   refinement_flag()=DO_NOTHING
   p_refinement_flag()=DO_NOTHING
   DoFs=
[2] src/fe/fe_map.C, line 1866, compiled Apr 29 2019 at 13:59:12
WARNING: Newton scheme has not converged in 12 iterations:
   physical_point=(x,y,z)=(    0.27,   0.2775,        0)   physical_guess=(x,y,z)=(    0.27,   0.2775,        0)   dp=(x,y,z)=(-1.21944e-13, -1.12163e-06,        0)   p=(x,y,z)=(      -1, -0.577351,        0)   error=1.12163e-06   in element 16514
  Elem Information
   id()=16514, processor_id()=2
   type()=QUAD4
   dim()=2
   n_nodes()=4
    0  Node id()=13234, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
    1  Node id()=604, processor_id()=2, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=(0/0/13470) (0/1/13471) (0/2/13472) (0/3/13473) (0/4/13474) (0/5/13475) (2/0/2245) (3/0/2245) 
    2  Node id()=13235, processor_id()=4294967295, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=
    3  Node id()=13233, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
   n_sides()=4
    neighbor(0)=nullptr
    neighbor(1)=nullptr
    neighbor(2)=nullptr
    neighbor(3)=nullptr
   hmin()=9.89829e-11, hmax()=0.001875
   volume()=1.85593e-13
   active()=1, ancestor()=0, subactive()=0, has_children()=0
   parent()=672
   level()=4, p_level()=0
   refinement_flag()=DO_NOTHING
   p_refinement_flag()=DO_NOTHING
   DoFs=
[2] src/fe/fe_map.C, line 1866, compiled Apr 29 2019 at 13:59:12
WARNING: Newton scheme has not converged in 13 iterations:
   physical_point=(x,y,z)=(    0.27,   0.2775,        0)   physical_guess=(x,y,z)=(    0.27,   0.2775,        0)   dp=(x,y,z)=(2.43888e-13, 2.24326e-06,        0)   p=(x,y,z)=(      -1, -0.577349,        0)   error=2.24326e-06   in element 16514
  Elem Information
   id()=16514, processor_id()=2
   type()=QUAD4
   dim()=2
   n_nodes()=4
    0  Node id()=13234, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
    1  Node id()=604, processor_id()=2, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=(0/0/13470) (0/1/13471) (0/2/13472) (0/3/13473) (0/4/13474) (0/5/13475) (2/0/2245) (3/0/2245) 
    2  Node id()=13235, processor_id()=4294967295, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=
    3  Node id()=13233, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
   n_sides()=4
    neighbor(0)=nullptr
    neighbor(1)=nullptr
    neighbor(2)=nullptr
    neighbor(3)=nullptr
   hmin()=9.89829e-11, hmax()=0.001875
   volume()=1.85593e-13
   active()=1, ancestor()=0, subactive()=0, has_children()=0
   parent()=672
   level()=4, p_level()=0
   refinement_flag()=DO_NOTHING
   p_refinement_flag()=DO_NOTHING
   DoFs=
[2] src/fe/fe_map.C, line 1866, compiled Apr 29 2019 at 13:59:12
WARNING: Newton scheme has not converged in 14 iterations:
   physical_point=(x,y,z)=(    0.27,   0.2775,        0)   physical_guess=(x,y,z)=(    0.27,   0.2775,        0)   dp=(x,y,z)=(-1.21944e-13, -1.12163e-06,        0)   p=(x,y,z)=(      -1, -0.57735,        0)   error=1.12163e-06   in element 16514
  Elem Information
   id()=16514, processor_id()=2
   type()=QUAD4
   dim()=2
   n_nodes()=4
    0  Node id()=13234, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
    1  Node id()=604, processor_id()=2, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=(0/0/13470) (0/1/13471) (0/2/13472) (0/3/13473) (0/4/13474) (0/5/13475) (2/0/2245) (3/0/2245) 
    2  Node id()=13235, processor_id()=4294967295, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=
    3  Node id()=13233, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
   n_sides()=4
    neighbor(0)=nullptr
    neighbor(1)=nullptr
    neighbor(2)=nullptr
    neighbor(3)=nullptr
   hmin()=9.89829e-11, hmax()=0.001875
   volume()=1.85593e-13
   active()=1, ancestor()=0, subactive()=0, has_children()=0
   parent()=672
   level()=4, p_level()=0
   refinement_flag()=DO_NOTHING
   p_refinement_flag()=DO_NOTHING
   DoFs=
[2] src/fe/fe_map.C, line 1866, compiled Apr 29 2019 at 13:59:12
WARNING: Newton scheme has not converged in 15 iterations:
   physical_point=(x,y,z)=(    0.27,   0.2775,        0)   physical_guess=(x,y,z)=(    0.27,   0.2775,        0)   dp=(x,y,z)=(-1.21944e-13, -1.12163e-06,        0)   p=(x,y,z)=(      -1, -0.577351,        0)   error=1.12163e-06   in element 16514
  Elem Information
   id()=16514, processor_id()=2
   type()=QUAD4
   dim()=2
   n_nodes()=4
    0  Node id()=13234, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
    1  Node id()=604, processor_id()=2, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=(0/0/13470) (0/1/13471) (0/2/13472) (0/3/13473) (0/4/13474) (0/5/13475) (2/0/2245) (3/0/2245) 
    2  Node id()=13235, processor_id()=4294967295, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=
    3  Node id()=13233, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
   n_sides()=4
    neighbor(0)=nullptr
    neighbor(1)=nullptr
    neighbor(2)=nullptr
    neighbor(3)=nullptr
   hmin()=9.89829e-11, hmax()=0.001875
   volume()=1.85593e-13
   active()=1, ancestor()=0, subactive()=0, has_children()=0
   parent()=672
   level()=4, p_level()=0
   refinement_flag()=DO_NOTHING
   p_refinement_flag()=DO_NOTHING
   DoFs=
[2] src/fe/fe_map.C, line 1866, compiled Apr 29 2019 at 13:59:12
WARNING: Newton scheme has not converged in 16 iterations:
   physical_point=(x,y,z)=(    0.27,   0.2775,        0)   physical_guess=(x,y,z)=(    0.27,   0.2775,        0)   dp=(x,y,z)=(2.43888e-13, 2.24326e-06,        0)   p=(x,y,z)=(      -1, -0.577349,        0)   error=2.24326e-06   in element 16514
  Elem Information
   id()=16514, processor_id()=2
   type()=QUAD4
   dim()=2
   n_nodes()=4
    0  Node id()=13234, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
    1  Node id()=604, processor_id()=2, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=(0/0/13470) (0/1/13471) (0/2/13472) (0/3/13473) (0/4/13474) (0/5/13475) (2/0/2245) (3/0/2245) 
    2  Node id()=13235, processor_id()=4294967295, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=
    3  Node id()=13233, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
   n_sides()=4
    neighbor(0)=nullptr
    neighbor(1)=nullptr
    neighbor(2)=nullptr
    neighbor(3)=nullptr
   hmin()=9.89829e-11, hmax()=0.001875
   volume()=1.85593e-13
   active()=1, ancestor()=0, subactive()=0, has_children()=0
   parent()=672
   level()=4, p_level()=0
   refinement_flag()=DO_NOTHING
   p_refinement_flag()=DO_NOTHING
   DoFs=
[2] src/fe/fe_map.C, line 1866, compiled Apr 29 2019 at 13:59:12
WARNING: Newton scheme has not converged in 17 iterations:
   physical_point=(x,y,z)=(    0.27,   0.2775,        0)   physical_guess=(x,y,z)=(    0.27,   0.2775,        0)   dp=(x,y,z)=(-1.21944e-13, -1.12163e-06,        0)   p=(x,y,z)=(      -1, -0.57735,        0)   error=1.12163e-06   in element 16514
  Elem Information
   id()=16514, processor_id()=2
   type()=QUAD4
   dim()=2
   n_nodes()=4
    0  Node id()=13234, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
    1  Node id()=604, processor_id()=2, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=(0/0/13470) (0/1/13471) (0/2/13472) (0/3/13473) (0/4/13474) (0/5/13475) (2/0/2245) (3/0/2245) 
    2  Node id()=13235, processor_id()=4294967295, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=
    3  Node id()=13233, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
   n_sides()=4
    neighbor(0)=nullptr
    neighbor(1)=nullptr
    neighbor(2)=nullptr
    neighbor(3)=nullptr
   hmin()=9.89829e-11, hmax()=0.001875
   volume()=1.85593e-13
   active()=1, ancestor()=0, subactive()=0, has_children()=0
   parent()=672
   level()=4, p_level()=0
   refinement_flag()=DO_NOTHING
   p_refinement_flag()=DO_NOTHING
   DoFs=
WARNING: Newton scheme has not converged in 21 iterations:
   physical_point=(x,y,z)=(    0.27,   0.2775,        0)   physical_guess=(x,y,z)=(    0.27,   0.2775,        0)   dp=(x,y,z)=(-1.21944e-13, -1.12163e-06,        0)   p=(x,y,z)=(      -1, -0.577351,        0)   error=1.12163e-06   in element 16514
  Elem Information
   id()=16514, processor_id()=2
   type()=QUAD4
   dim()=2
   n_nodes()=4
    0  Node id()=13234, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
    1  Node id()=604, processor_id()=2, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=(0/0/13470) (0/1/13471) (0/2/13472) (0/3/13473) (0/4/13474) (0/5/13475) (2/0/2245) (3/0/2245) 
    2  Node id()=13235, processor_id()=4294967295, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=
    3  Node id()=13233, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
   n_sides()=4
    neighbor(0)=nullptr
    neighbor(1)=nullptr
    neighbor(2)=nullptr
    neighbor(3)=nullptr
   hmin()=9.89829e-11, hmax()=0.001875
   volume()=1.85593e-13
   active()=1, ancestor()=0, subactive()=0, has_children()=0
   parent()=672
   level()=4, p_level()=0
   refinement_flag()=DO_NOTHING
   p_refinement_flag()=DO_NOTHING
   DoFs=
ERROR: Newton scheme FAILED to converge in 21 iterations in element 16514 for physical point = (x,y,z)=(    0.27,   0.2775,        0)
  Elem Information
   id()=16514, processor_id()=2
   type()=QUAD4
   dim()=2
   n_nodes()=4
    0  Node id()=13234, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
    1  Node id()=604, processor_id()=2, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=(0/0/13470) (0/1/13471) (0/2/13472) (0/3/13473) (0/4/13474) (0/5/13475) (2/0/2245) (3/0/2245) 
    2  Node id()=13235, processor_id()=4294967295, Point=(x,y,z)=(0.271875,   0.2775,        0)
    DoFs=
    3  Node id()=13233, processor_id()=4294967295, Point=(x,y,z)=(    0.27,   0.2775,        0)
    DoFs=
   n_sides()=4
    neighbor(0)=nullptr
    neighbor(1)=nullptr
    neighbor(2)=nullptr
    neighbor(3)=nullptr
   hmin()=9.89829e-11, hmax()=0.001875
   volume()=1.85593e-13
   active()=1, ancestor()=0, subactive()=0, has_children()=0
   parent()=672
   level()=4, p_level()=0
   refinement_flag()=DO_NOTHING
   p_refinement_flag()=DO_NOTHING
   DoFs=
Exiting...
[2] src/fe/fe_map.C, line 1905, compiled Apr 29 2019 at 13:59:12

> On May 1, 2019, at 9:29 PM, Stogner, Roy H <roy...@ic...> wrote:
> 
> 
> 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