Dear libMesh users.
I was trying to run this very simple test case where a 1D mesh is loaded
initially, and after inserting a new element in the FE-grid (and probably
some new nodes) the EquationSystems reinit() routine seems to crash. I
copy/paste below a test code to demonstrate the case:
//================================================
int main (int argc, char* argv[]) {
LibMeshInit init(argc, argv);
Mesh msh(init.comm(), 1);
EquationSystems es(msh);
es.add_system<System>("a_system");
es.get_system<System>("a_system").add_variable("a_variable");
GmshIO(msh).read("a_1D_mesh.msh");
msh.prepare_for_use();
msh.print_info();
msh.boundary_info->print_info();
es.init();
es.print_info();
{
Elem* e = Elem::build(EDGE2).release();
e->set_id(msh.max_elem_id());
e->set_node(0) = msh.node_ptr(2);
e->set_node(1) = msh.node_ptr(8);
msh.add_elem(e);
msh.prepare_for_use();
msh.print_info();
msh.boundary_info->print_info();
es.reinit();
es.print_info();
}
return 0;
}
//================================================
the error message I get from valgrind is:
==8907== Invalid read of size 4
==8907== at 0x5AA886: libMesh::DofMap::reinit(libMesh::MeshBase&) (in
/home/vvavoura/Work/A_TEST/a_test.dbg.exe)
==8907== by 0x5B09F0:
libMesh::DofMap::distribute_dofs(libMesh::MeshBase&) (in
/home/vvavoura/Work/A_TEST/a_test.dbg.exe)
==8907== by 0x4FEEAC: libMesh::EquationSystems::reinit() (in
/home/vvavoura/Work/A_TEST/a_test.dbg.exe)
==8907== by 0x41143E: main (in /home/vvavoura/Work/A_TEST/a_test.dbg.exe)
==8907== Address 0x4 is not stack'd, malloc'd or (recently) free'd
Cannot figure out what I did wrong :( Also I get the same error if I just
add a new node in the mesh along with a new 1D element.
Many thanks in advance.
cheers,
Vas
------------------------------------------------------------------------------
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
