Hi all,
I'm using libMesh as part of a tool that solves some E+M problems. One
task we have is to find the current distribution through different
conductors. Overall, this task is pretty similar to example 14 from the
website.
Many of the conductors we're studying are small -- on the order of a 10s
of microns wide. We had been performing our calculations in units of
microns, but recently changed to meters and started to see crashes in
equation_systems.reinit() when running in parallel. These crashes happen
in both our code and in example 14, though in slightly different places
(tracebacks follow). The attached input file will cause example 14 to
crash.
Any thoughts on the best approach to take here?
Thanks in advance,
..david
Traceback 1: ex14
Program received signal SIGSEGV, Segmentation fault.
...
#0 0x00007ffff746105f in DofObject::n_comp (this=0x917f20, s=0, var=0)
at libmesh/libmesh-0.6.4/libmesh/include/base/dof_object.h:558
#1 0x00007ffff7456df7 in DofMap::reinit (this=0x8c7c70, mesh=...)
at src/base/dof_map.C:489
#2 0x00007ffff74574b1 in DofMap::distribute_dofs (this=0x8c7c70, mesh=...)
at src/base/dof_map.C:669
#3 0x00007ffff78b1482 in EquationSystems::reinit (this=0x7fffffffdfa0)
at src/solvers/equation_systems.C:166
#4 0x00000000004308f7 in main (argc=2, argv=0x7fffffffe5a8)
at /home/vanwie/libmesh_test/ex14/ex14.cpp:404
Traceback 2: our code
Program received signal SIGSEGV, Segmentation fault.
...
#0 0x00007ffff0be3bc0 in VecScatterCreate_PtoP ()
from /usr/lib/libpetscvec.so.3.0.0
#1 0x00007ffff0bf00cf in VecScatterCreate ()
from /usr/lib/libpetscvec.so.3.0.0
#2 0x00007ffff5ee5274 in PetscVector<double>::localize (this=0x95ad60,
v_local_in=..., send_list=...) at src/numerics/petsc_vector.C:768
#3 0x00007ffff5f66768 in System::restrict_vectors (this=0x94dc20)
at src/solvers/system.C:246
#4 0x00007ffff5f66799 in System::prolong_vectors (this=0x94dc20)
at src/solvers/system.C:257
#5 0x00007ffff5f3b4fc in EquationSystems::reinit (this=0x94de20)
at src/solvers/equation_systems.C:180
...
DEAL 003:003
3 # Num. Elements
8 # Num. Nodes
12 # Sum of Element Weights (= 3*4)
0 # Num. Boundary Conds.
65536 # String Size (ignore)
1 # Num. Element Blocks.
5 # Element types in each block.
3 # Num. of elements in each block.
Id String
Title String
3 4 1 0
4 5 2 1
6 7 4 3
-0.0001 0.0001 0.
0. 0.0001 0.
0.0001 0.0001 0.
-0.0001 0. 0.
0. 0. 0.
0.0001 0. 0.
-0.0001 -0.0001 0.
0. -0.0001 0.
------------------------------------------------------------------------------
The Palm PDK Hot Apps Program offers developers who use the
Plug-In Development Kit to bring their C/C++ apps to Palm for a share
of $1 Million in cash or HP Products. Visit us here for more details:
http://ad.doubleclick.net/clk;226879339;13503038;l?
http://clk.atdmt.com/CRS/go/247765532/direct/01/
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
