I will think about it. Not sure when I will come up with a solution. Clearly in the example you sent the bug surface because capillaryForce >> elasticForce, which is also wrong from a physical point of view (I mean, capillary forces will never generate 50% overlap). I hope the problem is similar in your other simulation, else we have two bugs instead of one...
-- You received this bug notification because you are a member of Yade developers, which is the registrant for Yade. https://bugs.launchpad.net/bugs/966186 Title: fusionDetection does not work, when new contacts are created with hertz model Status in Yet Another Dynamic Engine: New Bug description: This bug is related to a question. https://answers.launchpad.net/yade/+question/190529 I attached a script to reproduce the bug, when fusionDetection=True and hertz model is used. The bug occurs, when contacts are created. In the script I created two fixed particles and one above the two. The upper particle falls down on the other two particles. At the beginning of the script you will find these lines: fusion_detection_on = 1 hertz_model_on = 0 All combinations with 0 and 1 are working, except this one: fusion_detection_on = 1 hertz_model_on = 1 The problem is somewhere in line 223 to 301 here: https://github.com/jakob- ifgt/trunk/blob/master/pkg/dem/Law2_ScGeom_CapillaryPhys_Capillarity.cpp Debugging gives: me@debian ~/YADE/my-yade-projects/fusionDetection-bug >/home/me/YADE/YADE-CapTest-Debug/bin/yade-unknown --debug check-fusion-detection-bug.py Welcome to Yade unknown (debug build) TCP python prompt on localhost:9000, auth cookie `ackess' XMLRPC info provider on http://localhost:21000 Running script check-fusion-detection-bug.py press Play button [[ ^L clears screen, ^U kills line. F12 controller, F11 3d view, F10 both, F9 generator, F8 plot. ]] Yade [1]: SIGSEGV/SIGABRT handler called; gdb batch file is `/tmp/yade-pCNX6x/tmp-0' GNU gdb (GDB) 7.0.1-debian Copyright (C) 2009 Free Software Foundation, Inc. License GPLv3+: GNU GPL version 3 or later <http://gnu.org/licenses/gpl.html> This is free software: you are free to change and redistribute it. There is NO WARRANTY, to the extent permitted by law. Type "show copying" and "show warranty" for details. This GDB was configured as "x86_64-linux-gnu". For bug reporting instructions, please see: <http://www.gnu.org/software/gdb/bugs/>. [Thread debugging using libthread_db enabled] [New Thread 0x7eff25eaa700 (LWP 7040)] [New Thread 0x7eff36769700 (LWP 7039)] [New Thread 0x7eff36f6a700 (LWP 7038)] 0x00007eff5a5159d3 in *__GI___poll (fds=<value optimized out>, nfds=<value optimized out>, timeout=47) at ../sysdeps/unix/sysv/linux/poll.c:87 87 ../sysdeps/unix/sysv/linux/poll.c: Datei oder Verzeichnis nicht gefunden. in ../sysdeps/unix/sysv/linux/poll.c Thread ID 0 not known. Thread 4 (Thread 0x7eff36f6a700 (LWP 7038)): #0 0x00007eff5a51a1a3 in select () at ../sysdeps/unix/syscall-template.S:82 #1 0x00000000004f7fd3 in ?? () #2 0x00000000004a7ba5 in PyEval_EvalFrameEx () #3 0x00000000004a95c1 in PyEval_EvalCodeEx () #4 0x00000000004a7752 in PyEval_EvalFrameEx () #5 0x00000000004a95c1 in PyEval_EvalCodeEx () #6 0x0000000000538a10 in ?? () #7 0x000000000041ef47 in PyObject_Call () #8 0x0000000000427c1f in ?? () #9 0x000000000041ef47 in PyObject_Call () #10 0x00000000004a1a53 in PyEval_CallObjectWithKeywords () #11 0x00007eff36f74cab in ?? () from /usr/lib/pymodules/python2.6/sip.so #12 0x00007eff3721b498 in ?? () from /usr/lib/pymodules/python2.6/PyQt4/QtCore.so #13 0x00007eff3724b68a in ?? () from /usr/lib/pymodules/python2.6/PyQt4/QtCore.so #14 0x00007eff565c5e15 in ?? () from /usr/lib/libQtCore.so.4 #15 0x00007eff5b64f8ca in start_thread (arg=<value optimized out>) at pthread_create.c:300 #16 0x00007eff5a52092d in clone () at ../sysdeps/unix/sysv/linux/x86_64/clone.S:112 #17 0x0000000000000000 in ?? () Thread 3 (Thread 0x7eff36769700 (LWP 7039)): #0 0x00007eff5a51a1a3 in select () at ../sysdeps/unix/syscall-template.S:82 #1 0x00000000004f7fd3 in ?? () #2 0x00000000004a7ba5 in PyEval_EvalFrameEx () #3 0x00000000004a95c1 in PyEval_EvalCodeEx () #4 0x00000000004a7752 in PyEval_EvalFrameEx () #5 0x00000000004a95c1 in PyEval_EvalCodeEx () #6 0x0000000000538a10 in ?? () #7 0x000000000041ef47 in PyObject_Call () #8 0x0000000000427c1f in ?? () #9 0x000000000041ef47 in PyObject_Call () #10 0x00000000004a1a53 in PyEval_CallObjectWithKeywords () #11 0x00007eff36f74cab in ?? () from /usr/lib/pymodules/python2.6/sip.so #12 0x00007eff3721b498 in ?? () from /usr/lib/pymodules/python2.6/PyQt4/QtCore.so #13 0x00007eff3724b68a in ?? () from /usr/lib/pymodules/python2.6/PyQt4/QtCore.so #14 0x00007eff565c5e15 in ?? () from /usr/lib/libQtCore.so.4 #15 0x00007eff5b64f8ca in start_thread (arg=<value optimized out>) at pthread_create.c:300 #16 0x00007eff5a52092d in clone () at ../sysdeps/unix/sysv/linux/x86_64/clone.S:112 #17 0x0000000000000000 in ?? () Current language: auto The current source language is "auto; currently asm". Thread 2 (Thread 0x7eff25eaa700 (LWP 7040)): #0 0x00007eff5a4f163d in __libc_waitpid (pid=7041, stat_loc=<value optimized out>, options=0) at ../sysdeps/unix/sysv/linux/waitpid.c:41 #1 0x00007eff5a48f2c9 in do_system (line=<value optimized out>) at ../sysdeps/posix/system.c:149 #2 0x00007eff5a48f600 in __libc_system (line=<value optimized out>) at ../sysdeps/posix/system.c:190 #3 0x00007eff5976488a in crashHandler (sig=11) at core/main/pyboot.cpp:45 #4 <signal handler called> #5 0x00007eff4083d07b in Law2_ScGeom_CapillaryPhys_Capillarity::checkFusion (this=0x7eff3077ac50) at /home/me/YADE/trunk/pkg/dem/Law2_ScGeom_CapillaryPhys_Capillarity.cpp:295 #6 0x00007eff4083c4b7 in Law2_ScGeom_CapillaryPhys_Capillarity::action (this=0x7eff3077ac50) at /home/me/YADE/trunk/pkg/dem/Law2_ScGeom_CapillaryPhys_Capillarity.cpp:178 #7 0x00007eff50f203a5 in Scene::moveToNextTimeStep (this=0xbf3250) at /home/me/YADE/trunk/core/Scene.cpp:97 #8 0x00007eff50f278e5 in SimulationFlow::singleAction (this=0xb70550) at /home/me/YADE/trunk/core/SimulationFlow.cpp:20 #9 0x00007eff50f28e7e in ThreadWorker::callSingleAction (this=0xb70550) at /home/me/YADE/trunk/core/ThreadWorker.cpp:71 #10 0x00007eff50f283e3 in ThreadRunner::call (this=0xbf66a0) at /home/me/YADE/trunk/core/ThreadRunner.cpp:53 #11 0x00007eff50f28125 in ThreadRunner::run (this=0xbf66a0) at /home/me/YADE/trunk/core/ThreadRunner.cpp:27 #12 0x00007eff50ffd23d in boost::_mfi::mf0<void, ThreadRunner>::operator() (this=0x7eff30ea4488, p=0xbf66a0) at /usr/include/boost/bind/mem_fn_template.hpp:49 #13 0x00007eff50fecf18 in boost::_bi::list1<boost::_bi::value<ThreadRunner*> >::operator()<boost::_mfi::mf0<void, ThreadRunner>, boost::_bi::list0> (this=0x7eff30ea4498, f=..., a=...) at /usr/include/boost/bind/bind.hpp:253 #14 0x00007eff50fdbf3b in boost::_bi::bind_t<void, boost::_mfi::mf0<void, ThreadRunner>, boost::_bi::list1<boost::_bi::value<ThreadRunner*> > >::operator() (this=0x7eff30ea4488) at /usr/include/boost/bind/bind_template.hpp:20 #15 0x00007eff50fc8284 in boost::detail::function::void_function_obj_invoker0<boost::_bi::bind_t<void, boost::_mfi::mf0<void, ThreadRunner>, boost::_bi::list1<boost::_bi::value<ThreadRunner*> > >, void>::invoke (function_obj_ptr=...) at /usr/include/boost/function/function_template.hpp:153 #16 0x00007eff510a3ef3 in boost::function0<void>::operator() (this=0x7eff30ea4480) at /usr/include/boost/function/function_template.hpp:1013 #17 0x00007eff5107c9f2 in boost::detail::thread_data<boost::function0<void> >::run (this=0x7eff30ea4350) at /usr/include/boost/thread/detail/thread.hpp:56 #18 0x00007eff58e86200 in thread_proxy () from /usr/lib/libboost_thread.so.1.42.0 #19 0x00007eff5b64f8ca in start_thread (arg=<value optimized out>) at pthread_create.c:300 #20 0x00007eff5a52092d in clone () at ../sysdeps/unix/sysv/linux/x86_64/clone.S:112 #21 0x0000000000000000 in ?? () Thread 1 (Thread 0x7eff5ba5b700 (LWP 7035)): #0 0x00007eff5a5159d3 in *__GI___poll (fds=<value optimized out>, nfds=<value optimized out>, timeout=47) at ../sysdeps/unix/sysv/linux/poll.c:87 #1 0x00007eff51ec02c9 in ?? () from /lib/libglib-2.0.so.0 #2 0x00007eff51ec071c in g_main_context_iteration () from /lib/libglib-2.0.so.0 #3 0x00007eff566de6b3 in QEventDispatcherGlib::processEvents(QFlags<QEventLoop::ProcessEventsFlag>) () from /usr/lib/libQtCore.so.4 #4 0x00007eff56c7519e in ?? () from /usr/lib/libQtGui.so.4 #5 0x00007eff566b39c2 in QEventLoop::processEvents(QFlags<QEventLoop::ProcessEventsFlag>) () from /usr/lib/libQtCore.so.4 #6 0x00007eff566b3d9c in QEventLoop::exec(QFlags<QEventLoop::ProcessEventsFlag>) () from /usr/lib/libQtCore.so.4 #7 0x00007eff566b8a2b in QCoreApplication::exec() () from /usr/lib/libQtCore.so.4 #8 0x00007eff3721dee5 in ?? () from /usr/lib/pymodules/python2.6/PyQt4/QtCore.so #9 0x00007eff3be75991 in ?? () from /usr/lib/python2.6/lib-dynload/readline.so #10 0x0000000000526a52 in PyOS_Readline () #11 0x000000000049da86 in ?? () #12 0x00000000004a7ba5 in PyEval_EvalFrameEx () #13 0x00000000004a95c1 in PyEval_EvalCodeEx () #14 0x00000000004a7752 in PyEval_EvalFrameEx () #15 0x00000000004a95c1 in PyEval_EvalCodeEx () #16 0x00000000004a7752 in PyEval_EvalFrameEx () #17 0x00000000004a95c1 in PyEval_EvalCodeEx () #18 0x00000000004a7752 in PyEval_EvalFrameEx () #19 0x00000000004a95c1 in PyEval_EvalCodeEx () #20 0x0000000000538a10 in ?? () #21 0x000000000041ef47 in PyObject_Call () #22 0x0000000000427c1f in ?? () #23 0x000000000041ef47 in PyObject_Call () #24 0x0000000000426edd in ?? () #25 0x000000000041ef47 in PyObject_Call () #26 0x00000000004a72b8 in PyEval_EvalFrameEx () #27 0x00000000004a95c1 in PyEval_EvalCodeEx () #28 0x00000000004a7752 in PyEval_EvalFrameEx () #29 0x00000000004a95c1 in PyEval_EvalCodeEx () #30 0x00000000004a9692 in PyEval_EvalCode () #31 0x00000000004c98be in PyRun_FileExFlags () #32 0x00000000004c9ad4 in PyRun_SimpleFileExFlags () #33 0x000000000041a6bd in Py_Main () #34 0x00007eff5a46fc8d in __libc_start_main (main=<value optimized out>, argc=<value optimized out>, ubp_av=<value optimized out>, init=<value optimized out>, fini=<value optimized out>, rtld_fini=<value optimized out>, stack_end=0x7fff4de53f28) at libc-start.c:228 #35 0x00000000004198d9 in _start () Current language: auto The current source language is "auto; currently c". segmentation fault To manage notifications about this bug go to: https://bugs.launchpad.net/yade/+bug/966186/+subscriptions _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : yade-dev@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp