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.an...@gmail.com> > The problem seems in quaternions. Will try to fix it. > ______________________________ > > Anton Gladkyy > > > 2010/3/7 Anton Gladky <gladky.an...@gmail.com> > > 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.an...@gmail.com> >> >> Ok, thanks. Will try. >>> ______________________________ >>> >>> Anton Gladkyy >>> >>> >>> 2010/3/2 Václav Šmilauer <eudo...@arcig.cz> >>> >>> >>>>> 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@lists.launchpad.net >>>> Unsubscribe : >>>> https://launchpad.net/~yade-dev<https://launchpad.net/%7Eyade-dev> >>>> More help : https://help.launchpad.net/ListHelp >>>> >>> >>> >> >
_______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : yade-dev@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
