yade-dev team mailing list archive

[Anton Gladky <gladky.anton@xxxxxxxxx>]

Hi, Vaclav!

For now, I have defined, that nan appears in
Ig2_Sphere_Sphere_Dem3DofGeom::go in Dem3DofGeom_SphereSphere.cpp

176: Vector3r normal=(state2.pos+shift2)-state1.pos;

state1.pos gives nan when particles penetrate. Could you not say, where
should I go next to find the problem? In state.hpp?
I am a little bit lost with that.

Thanks
______________________________

Anton Gladkyy


2010/3/2 Anton Gladky <gladky.anton@xxxxxxxxx>

> Ok, thanks. Will try.
> ______________________________
>
> Anton Gladkyy
>
>
> 2010/3/2 Václav Šmilauer <eudoxos@xxxxxxxx>
>
>
>>> the problem of flying particles is definitely in NewtonIntegrator.
>>> I have parsed the code twice, but everything seems is ok there.
>>> Also I added some verification tests in test2.cpp, but...
>>>
>>> Is there any way to catch the problem, if the crash does not appears?
>>>
>>
>> You have to start with very simple simulation (such as sphere-sphere) and
>> then watch wehyre NaNs start to appear; that is place that will be
>> responsible for disappearing particles; it might not be integrator itself,
>> but if it gets NaN acceleration (for whatever reason), then it will also
>> compute NaN position, since NaNs propagate through computations.
>>
>>
>> Cheers, Vaclav
>>
>>
>> _______________________________________________
>> Mailing list: https://launchpad.net/~yade-dev<https://launchpad.net/%7Eyade-dev>
>> Post to     : yade-dev@xxxxxxxxxxxxxxxxxxx
>> Unsubscribe : https://launchpad.net/~yade-dev<https://launchpad.net/%7Eyade-dev>
>> More help   : https://help.launchpad.net/ListHelp
>>
>
>