Can you build and link against libmesh with METHOD=devel or, preferably, METHOD=dbg? Either of those will turn on asserts in libmesh and hopefully which you will hopefully trip instead of just a default. Right now you're linking against libmesh in opt mode and the segfault may be triggered much later than the actual problem.
> On Sep 25, 2015, at 1:31 PM, Harshad Sahasrabudhe <hsaha...@purdue.edu> wrote: > > Hi, > > I get a segfault in the function FEMap::compute_affine_map() when I try to > fe->reinit() an element. Here is the stack trace: > > #0 0x000000376b632625 in raise () from /lib64/libc.so.6 > #1 0x000000376b633e05 in abort () from /lib64/libc.so.6 > #2 0x000000376b670537 in __libc_message () from /lib64/libc.so.6 > #3 0x000000376b675e66 in malloc_printerr () from /lib64/libc.so.6 > #4 0x000000376b6789ba in _int_free () from /lib64/libc.so.6 > #5 0x00002b3d41a4a7f9 in libMesh::FEMap::compute_affine_map(unsigned int, > std::vector<double, std::allocator<double> > const&, libMesh::Elem const*) > () from > /depot/ncn/apps/conte/PETSc35_Intel15_no-mic_debug/libs/libmesh/libmesh/.libs/libmesh_opt.so.0 > #6 0x00002b3d41a4b35a in libMesh::FEMap::compute_map(unsigned int, > std::vector<double, std::allocator<double> > const&, libMesh::Elem const*) > () from > /depot/ncn/apps/conte/PETSc35_Intel15_no-mic_debug/libs/libmesh/libmesh/.libs/libmesh_opt.so.0 > #7 0x00002b3d4196d1c9 in libMesh::FE<3u, > (libMesh::FEFamily)0>::reinit(libMesh::Elem const*, > std::vector<libMesh::Point, std::allocator<libMesh::Point> > const*, > std::vector<double, std::allocator<double> > const*) () from > /depot/ncn/apps/conte/PETSc35_Intel15_no-mic_debug/libs/libmesh/libmesh/.libs/libmesh_opt.so.0 > > Anyone has any idea where this could be happening in that function? I'm > sorry I can't get the exact line numbers because of the way we compile our > code, which leads to the loss of some debugging information, and I have no > way of getting around it. > > Thanks! > Harshad > ------------------------------------------------------------------------------ > _______________________________________________ > Libmesh-users mailing list > Libmesh-users@lists.sourceforge.net > https://lists.sourceforge.net/lists/listinfo/libmesh-users ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
