I found a way to get it in a stable condition with calm() function. After that no calm function is active and the model behaves like this, when running it three times (100k steps):
me@debian ~/YADE/my-yade-projects/09-converted-PFC-model >yadetest -j2 test.py <skip> actual porosity: 45.478219 percent, time step is 6.952821e-08 ball with id 52 has maximum velocity of 8.949591e-04 unbalanced force: 0.105333 actual porosity: 45.478214 percent, time step is 7.086720e-08 ball with id 569 has maximum velocity of 2.067789e-04 unbalanced force: 0.007260 exit() Do you really want to exit ([y]/n)? me@debian ~/YADE/my-yade-projects/09-converted-PFC-model >yadetest -j2 test.py <skip> actual porosity: 45.478219 percent, time step is 6.946380e-08 ball with id 52 has maximum velocity of 8.448021e-04 unbalanced force: 0.104514 actual porosity: 45.478214 percent, time step is 7.089416e-08 ball with id 569 has maximum velocity of 1.872520e-04 unbalanced force: 0.005917 exit() Do you really want to exit ([y]/n)? me@debian ~/YADE/my-yade-projects/09-converted-PFC-model >yadetest -j2 test.py <skip> actual porosity: 45.478220 percent, time step is 6.947337e-08 ball with id 52 has maximum velocity of 8.661632e-04 unbalanced force: 0.106689 actual porosity: 45.478214 percent, time step is 7.084662e-08 ball with id 569 has maximum velocity of 2.275644e-04 unbalanced force: 0.007067 Is this acceptable rate of indeterminism? I see there are unexpected high differences. It is always stable this time. -- You received this bug notification because you are a member of Yade developers, which is the registrant for Yade. https://bugs.launchpad.net/bugs/987164 Title: bug with multithreading Status in Yet Another Dynamic Engine: Invalid Bug description: This bug is related to a question: https://answers.launchpad.net/yade/+question/191904 Using multi-threading option in yade leads to instabilities, that are caused by indeterminism. To reproduce the bug, download attached file, extract and run 1-generation.py. You will get a save-file called save1-test.yade. With this file run 2-relaxation.py and wait until its ready (will take a few minutes, be patient). Then a new save file called save2-test.yade will be created. With this save file you can run 3-settling.py to reproduce the bug, see below: 1.) run on 1 core: me@debian ~/YADE/my-yade-projects/multi-threading-bug >yade 3-settling.py Welcome to Yade unknown TCP python prompt on localhost:9001, auth cookie `sauces' XMLRPC info provider on http://localhost:21001 Running script 3-settling.py actual porosity: 46.090125 percent ball with id 2 has maximum velocity of 0.002235 actual porosity: 41.579271 percent ball with id 7 has maximum velocity of 0.354973 actual porosity: 42.846368 percent ball with id 3 has maximum velocity of 0.180586 Yade: normal exit. Beendet 2.) run on more cores: me@debian ~/YADE/my-yade-projects/multi-threading-bug >yade -j4 3-settling.py Welcome to Yade unknown TCP python prompt on localhost:9001, auth cookie `ekcdsa' XMLRPC info provider on http://localhost:21001 Running script 3-settling.py actual porosity: 46.090125 percent ball with id 2 has maximum velocity of 0.002235 actual porosity: 43.217713 percent ball with id 0 has maximum velocity of 0.103208 actual porosity: 43.106072 percent ball with id 20 has maximum velocity of 0.164024 Yade: normal exit. Beendet 2.) run on more cores again, until instability occurs: me@debian ~/YADE/my-yade-projects/multi-threading-bug >yade -j4 3-settling.py Welcome to Yade unknown TCP python prompt on localhost:9000, auth cookie `ssacue' XMLRPC info provider on http://localhost:21000 Running script 3-settling.py actual porosity: 46.090125 percent ball with id 2 has maximum velocity of 0.002235 actual porosity: 42.783835 percent ball with id 5 has maximum velocity of 0.170237 actual porosity: 57.019197 percent !!!!!!!!!!!!!!!BUG LEADS TO INSTABILITY OF THE MODEL HERE !!!!!!!!!!! ball with id 10 has maximum velocity of 3.735641 pure virtual method called terminate called without an active exception Abgebrochen To manage notifications about this bug go to: https://bugs.launchpad.net/yade/+bug/987164/+subscriptions _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
