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.
--
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
