Dear all,
We are having trouble with the FE_DGQ<3> element (no problem occurs in
2d). Take as an example step-12 run in 3d, then the program gets stuck
in the constructor. Here is the output from gdb:
#0 0x00002aaab6ae64c4 in __lll_lock_wait () from /lib64/libpthread.so.0
#1 0x00002aaab6ae1e1a in _L_lock_1034 () from /lib64/libpthread.so.0
#2 0x00002aaab6ae1cdc in pthread_mutex_lock ()
from /lib64/libpthread.so.0
#3 0x00000000004a6120 in __gthread_mutex_lock (__mutex=0x729da0)
at
/bubo/sw/comp/gcc/4.5.0/lib/gcc/x86_64-unknown-linux-gnu/4.5.0/../../../../include/c++/4.5.0/x86_64-unknown-linux-gnu/bits/gthr-default.h:758
#4 0x00000000004ac570 in std::mutex::lock (this=0x729da0)
at
/bubo/sw/comp/gcc/4.5.0/lib/gcc/x86_64-unknown-linux-gnu/4.5.0/../../../../include/c++/4.5.0/mutex:88
#5 0x00000000004ac5d2 in dealii::Threads::Mutex::acquire
(this=0x729da0)
at /bubo/home/h8/martinkr/deal.II/include/deal.II/base/thread_management.h:390
#6 0x00000000004d136d in
dealii::Threads::internal::TaskDescriptor<void>::join (this=0x729d70)
at /bubo/home/h8/martinkr/deal.II/include/deal.II/base/thread_management.h:3981
#7 0x00000000004cd1f0 in dealii::Threads::Task<void>::join
(this=0x729000)
at /bubo/home/h8/martinkr/deal.II/include/deal.II/base/thread_management.h:4089
#8 0x00000000004c8502 in dealii::Threads::TaskGroup<void>::join_all
(this=0x7fffffffb510)
at /bubo/home/h8/martinkr/deal.II/include/deal.II/base/thread_management.h:5116
#9 0x00002aaaae391d6e in dealii::FETools::compute_embedding_matrices<3,
double, 3> (fe=...,
matrices=std::vector of length 7, capacity 7 = {...},
isotropic_only=false)
at /bubo/home/h8/martinkr/deal.II/source/fe/fe_tools.cc:810
#10 0x00002aaaae21f698 in dealii::FE_DGQ<3, 3>::FE_DGQ
(this=0x7fffffffda00, degree=1)
at /bubo/home/h8/martinkr/deal.II/source/fe/fe_dgq.cc:147
#11 0x00000000004adce0 in Step12<3>::Step12 (this=0x7fffffffb850) at
step-12.cc:247
#12 0x00000000004a6274 in main () at step-12.cc:810
deal.II is compiled with threads enabled and the compiler uses the
std::thread (since the flag -std=c++0x is set). If threads are disabled,
the problem does not occur.
Does anyone who recently edited the relevant sections in FETools
(compute_embedding_matrices) have an idea why this could happen? Any
guess?
Thanks!
Best regards,
Katharina and Martin
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
