Re: [Yade-dev] velocity==NaN!

2010-08-30 Thread Janek Kozicki
chiara modenese said: (by the date of Thu, 26 Aug 2010 16:27:50 +0100) > On 26 August 2010 13:25, Bruno Chareyre wrote: > > > I bet on contact law returning nan value (division by 0, or something like > > that). > > Then velocity is NaN when this force is integrated. > > > You win the bet,

Re: [Yade-dev] velocity==NaN!

2010-08-27 Thread Václav Šmilauer
> Interesting. Also remember that you can advance line by line and inspect > values with the debugger, with various conditional breakpoints. > I do that with kdevelop but I guess it can as well be done with other > gdb frontends or with gdb alone. Oh yes, of course. This is just much much easier.

Re: [Yade-dev] velocity==NaN!

2010-08-27 Thread Bruno Chareyre
In the last commit, I added sub-stepping functinality where you can advance by a single engine (rather than full iteration at once), that whould help in such cases. (I keep in mind to show force acting on Body in the Inspector) Interesting. Also remember that you can advance line by line an

Re: [Yade-dev] velocity==NaN!

2010-08-26 Thread Václav Šmilauer
> You win the bet, sorry. After a day of debugging, I realized I was > returning a negative tangential stiff, and this was causing the > problem with contact damping. Thxs for suggestions anyway. I see from > the mail archive that Janek had some problems with Nan values > (unrelated to my question

Re: [Yade-dev] velocity==NaN!

2010-08-26 Thread chiara modenese
On 26 August 2010 13:25, Bruno Chareyre wrote: > I bet on contact law returning nan value (division by 0, or something like > that). > Then velocity is NaN when this force is integrated. > You win the bet, sorry. After a day of debugging, I realized I was returning a negative tangential stiff, a

Re: [Yade-dev] velocity==NaN!

2010-08-26 Thread Bruno Chareyre
I bet on contact law returning nan value (division by 0, or something like that). Then velocity is NaN when this force is integrated. Cheers. Bruno On 25/08/10 20:12, chiara modenese wrote: Hi all, what are the reasons for which one should get the following message from a simulation? 9088

Re: [Yade-dev] velocity==NaN!

2010-08-26 Thread Jerome Duriez
Hello Chiara, I have no precise remark but you are indeed not the first one having problems with "Nan" velocities : see http://www.mail-archive.com/yade-dev@lists.launchpad.net/msg04519.html and http://www.mail-archive.com/yade-dev@lists.launchpad.net/msg04766.html (at least). Jerome

Re: [Yade-dev] velocity==NaN!

2010-08-25 Thread Anton Gladky
On 08/25/2010 08:34 PM, chiara modenese wrote: . but since this occurs only when I do a small a change in my law (say in a case I include contact damping in the other I do not) I guess if there could be any other cause related to some other parts of the code. I only guess if any other

Re: [Yade-dev] velocity==NaN!

2010-08-25 Thread chiara modenese
On 25 August 2010 19:12, chiara modenese wrote: > Hi all, > > what are the reasons for which one should get the following message from a > simulation? > > 9088 FATAL yade.ThreadRunner > /home/chia/Documents/yade-r2394-psd2/core/ThreadRunner.cpp:31 run: Exception > occured: > Body #5653 has veloc

[Yade-dev] velocity==NaN!

2010-08-25 Thread chiara modenese
Hi all, what are the reasons for which one should get the following message from a simulation? 9088 FATAL yade.ThreadRunner /home/chia/Documents/yade-r2394-psd2/core/ThreadRunner.cpp:31 run: Exception occured: Body #5653 has velocity==NaN! In my case it happens right after I launch a triaxial t