Hi, Folks --
I am getting an assertion failure when calling EquationSystems::init()
which which seems to be related to threading support when I add a
periodic boundary to a simple rectangular mesh. Everything seems to run
fine when I don't add the periodic boundary. I am not using multiple
threads (or I certainly don't think that I am!), so I'm not sure what is
going on here.
Here's the error and stack trace:
8 0x00dbd674 in BoolAcquire (this=0xbfff85cc, b...@0x1b4e7da) at threads.h:64
64 BoolAcquire(bool& b) : _b(b) { libmesh_assert(!_b); _b = true; }
(gdb) where
#0 0x97aeae42 in __kill ()
#1 0x97aeae34 in kill$UNIX2003 ()
#2 0x97b5d23a in raise ()
#3 0x97b69679 in abort ()
#4 0x0334c30f in __gnu_cxx::__verbose_terminate_handler ()
#5 0x03349ec9 in __cxxabiv1::__terminate ()
#6 0x03349f0c in std::terminate ()
#7 0x0334a00b in __cxa_throw ()
#8 0x00dbd674 in BoolAcquire (this=0xbfff85cc, b...@0x1b4e7da) at
threads.h:64
#9 0x01470936 in parallel_reduce<ConstNodeRange, <unnamed>::FindBBox>
(ran...@0xbfff8c40, bo...@0xbfff8bf4) at threads.h:292
#10 0x01466764 in MeshTools::bounding_box (me...@0xbfffd77c) at
src/mesh/mesh_tools.C:290
#11 0x0168afec in PointLocatorTree::init (this=0x35656c0,
build_type=Trees::NODES) at src/utils/point_locator_tree.C:117
#12 0x0168a96b in PointLocatorTree (this=0x35656c0, me...@0xbfffd77c,
master=0x0) at src/utils/point_locator_tree.C:42
#13 0x01686f03 in PointLocatorBase::build (t=MeshEnums::TREE,
me...@0xbfffd77c, master=0x0) at src/utils/point_locator_base.C:64
#14 0x01398e11 in MeshBase::point_locator (this=0xbfffd77c) at
src/mesh/mesh_base.C:293
#15 0x00dbd210 in PeriodicBoundaries::neighbor (this=0x354f764,
boundary_id=3, me...@0xbfffd77c, e=0x3532340, side=3) at
src/base/dof_map_constraints.C:1342
#16 0x00eb2ff3 in FEBase::compute_periodic_constraints
(constrain...@0x354f740, dof_m...@0x354f610, boundari...@0x354f764,
me...@0xbfffd77c, variable_number=0, elem=0x3532340) at
src/fe/fe_base.C:1983
#17 0x00fa3a86 in FEInterface::compute_periodic_constraints
(constrain...@0x354f740, dof_m...@0x354f610, boundari...@0x354f764,
me...@0xbfffd77c, variable_number=0, elem=0x3532340) at
src/fe/fe_interface.C:1995
#18 0x00dab668 in operator() (this=0xbfffc6b4, ran...@0xbfffc600) at
src/base/dof_map_constraints.C:89
#19 0x00dbd28f in parallel_for<ConstElemRange,
<unnamed>::ComputeConstraints> (ran...@0xbfffc600, bo...@0xbfffc6b4) at
threads.h:267
#20 0x00dabdf9 in DofMap::create_dof_constraints (this=0x354f610,
me...@0xbfffd77c) at src/base/dof_map_constraints.C:168
#21 0x01648058 in System::init_data (this=0x354f340) at
src/systems/system.C:205
#22 0x0162dfa0 in ExplicitSystem::init_data (this=0x354f340) at
src/systems/explicit_system.C:72
#23 0x01647fdc in System::init (this=0x354f340) at src/systems/system.C:182
#24 0x0161919c in EquationSystems::init (this=0xbfffd728) at
src/systems/equation_systems.C:116
Thanks for any suggestions,
-- Boyce
------------------------------------------------------------------------------
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
