yade-dev team mailing list archive

[Anton Gladky <gladky.anton@xxxxxxxxx>]

Hi, Vaclav,

Seems, r63 works normally for non-clumped particles.

But some clump members disappear after a first step. Is it somewhere in
clump.cpp?

Thanks
______________________________

Anton Gladkyy


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