I tried another script without facets (using aabbWalls) and without changing particle radii, it crashes after 5 or 6 steps (when neverErase=True):
5468 FATAL yade.ThreadRunner /home/me/YADE/yade/core/ThreadRunner.cpp:31 run: Exception occured: Body #3084 has velocity==NaN! I think the problem is the changes, that Bruno implemented with neverErase. All contacts remain active with this option. But we need a rupture distance, at which liquid bridges break and the contact gets inactive. Something like if (uN<0) && (uN>rupture_distance){ if (neverErase) { phys->kn = 0.0; phys->ks = 0.0; #...maybe more content here... phys->shearForce = phys->normalForce = Vector3r::Zero(); } else {scene->interactions->requestErase(id1,id2); return;} } if (uN<rupture_distance) {scene->interactions->requestErase(id1,id2); return;} -- You received this bug notification because you are a member of Yade developers, which is the registrant for Yade. https://bugs.launchpad.net/bugs/897237 Title: bug in neverErase?! Status in Yet Another Dynamic Engine: New Bug description: If I set neverErase=True Yade crashes with following output: O.run(3000,True) RuntimeError: std::exception and also: Traceback (most recent call last): File "/home/me/YADE/YADE2970/lib/yade-bzr2970/py/yade/qt/__init__.py", line 192, in refreshEvent 15:01:13 self.refreshValues() 15:01:30 File "/home/me/YADE/YADE2970/lib/yade-bzr2970/py/yade/qt/__init__.py", line 265, in refreshValues 15:01:34 s=int(t); ms=int(t*1000)%1000; us=int(t*1000000)%1000; ns=int(t*1000000000)%1000 ValueError: cannot convert float NaN to integer You can reproduce the bug with the attached script (similar to the one i submitted with version 2960). First try to set Law2_ScGeom_MindlinPhys_Mindlin(neverErase=False), this will work. Then set neverErase=True and yade crashes. To manage notifications about this bug go to: https://bugs.launchpad.net/yade/+bug/897237/+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