yade-dev team mailing list archive

[Anton Gladky <gladky.anton@xxxxxxxxx>]

The problem seems in quaternions. Will try to fix it.
______________________________

Anton Gladkyy


2010/3/7 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
>>>
>>
>>
>