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