yade-dev team mailing list archive

[Chareyre <897237@xxxxxxxxxxxxxxxxxx>]

Christian,
Yes. GSTS is fine after with kn/ks, but there are also junk values comming from the contact law itself apparently (after #7), maybe because of the other stiffness values.

Chiara,
Stable = !unstable ;)
If you don't see an explosion and NaN values in a few timesteps, then it is stable. There is a very narrow range of dt value where it is not exploding while not completely stable, something like dt +/- 10%.
So, if you get NaN velocities for dt>=T, set dt=0.8*T and you are on the safe side.
If you have, let's say, dt=0.01*T, it is stable too but you are wasting cpu 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/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