Hi, I am getting started in libmesh - I've been using it for a month now. I used libmesh introductory example 3 (libmesh_ex3) as a reference for my problem. The idea is to replace the mesh by one of my one and start from there.
I noticed that libmesh_ex3 uses TRI6 or QUAD elements. In my application I was using TRI3 and it did not work. I realized that the example does not work with TRI3 either. Is there a reason for that? Should I use TRI6 instead? Any advice? The stack trace: Assertion `oldvar == static_cast<Told>(static_cast<Tnew>(oldvar))' failed. oldvar = 4294967295 static_cast<Told>(static_cast<Tnew>(oldvar)) = 255 Stack frames: 14 0: libMesh::print_trace(std::ostream&) 1: libMesh::MacroFunctions::report_error(char const*, int, char const*, char const*) 2: obj/bin/libmesh_ex3() [0x4f36b1] 3: libMesh::Elem::set_p_level(unsigned int) 4: libMesh::DofMap::reinit(libMesh::MeshBase&) 5: libMesh::DofMap::distribute_dofs(libMesh::MeshBase&) 6: libMesh::System::init_data() 7: libMesh::ImplicitSystem::init_data() 8: libMesh::LinearImplicitSystem::init_data() 9: libMesh::System::init() 10: libMesh::EquationSystems::init() 11: obj/bin/libmesh_ex3() [0x50a5e5] 12: __libc_start_main 13: obj/bin/libmesh_ex3() [0x41f539] [0] /usr/local/include/libmesh/libmesh_common.h, line 554, compiled Aug 7 2017 at 23:25:17 application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0 [unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=1 : Thanks Renato ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
