On Fri, May 8, 2015 at 2:45 PM, John Peterson <[email protected]> wrote:

>
>
> On Thu, May 7, 2015 at 8:27 AM, John Peterson <[email protected]>
> wrote:
>
>>
>> Can you run this through a debugger and get a stack trace too?  And
>> please send me the mesh directly, the mailing list will eat attachments...
>>
>
> This code hits an assert for me in debug mode during
> EquationSystems::reinit():
>
> Assertion `elem->n_systems() == this->n_systems()' failed.
> elem->n_systems() = 0
> this->n_systems() = 1
>
> * thread #1: tid = 0x3161f7, 0x00000001035aa702
> libmesh_dbg.0.dylib`libMesh::print_trace(out_stream=0x0000000000000000) +
> 18 at print_trace.C:201, queue = 'com.apple.main-thread', stop reason =
> breakpoint 1.1
>   * frame #0: 0x00000001035aa702
> libmesh_dbg.0.dylib`libMesh::print_trace(out_stream=0x0000000000000000) +
> 18 at print_trace.C:201
>     frame #1: 0x00000001035a3830
> libmesh_dbg.0.dylib`libMesh::MacroFunctions::report_error(file=0x0000000104bf0503,
> line=187, date=0x0000000104bbb74d, time=0x0000000104bf0481) + 608 at
> libmesh_common.C:82
>     frame #2: 0x000000010407faa8
> libmesh_dbg.0.dylib`libMesh::EquationSystems::reinit(this=0x00007fff5fbfe3a0)
> + 4168 at equation_systems.C:187
>     frame #3: 0x00000001000023cb vasileios-dbg`main(argc=1,
> argv=0x00007fff5fbfe6f8) + 2123 at vasileios.cc:40
>     frame #4: 0x00007fff908835c9 libdyld.dylib`start + 1
>     frame #5: 0x00007fff908835c9 libdyld.dylib`start + 1
>
> This seems like a legitimate bug... it's like the System doesn't know
> about the new element that you manually added.  I tried
> calling e->set_n_systems(1); but that just led to a different assert
> elsewhere.
>
> We usually only add new elements through MeshRefinement before calling
> reinit(), but in that case they are flagged as JUST_REFINED and have a
> valid parent Elem... I wonder if changing the "level 0" mesh in between
> reinit() calls has just never been tested.
>

This is almost certainly the same bug as
https://github.com/libMesh/libmesh/issues/435, which I will update with
your test code that demonstrates the problem.

-- 
John
------------------------------------------------------------------------------
One dashboard for servers and applications across Physical-Virtual-Cloud 
Widest out-of-the-box monitoring support with 50+ applications
Performance metrics, stats and reports that give you Actionable Insights
Deep dive visibility with transaction tracing using APM Insight.
http://ad.doubleclick.net/ddm/clk/290420510;117567292;y
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to