Ok, thanks, will try. The previous problem with "disappearing" particles was, that normalization of NULL-vector in WM3 gives (0,0,0), but Eigen gives (nan,nan,nan). So I needed to add such check in yademath. I do not know, what implementation is more correct: WM3 or Eigen. ______________________________
Anton Gladkyy 2010/3/9 Václav Šmilauer <eudo...@arcig.cz> > > Seems, r63 works normally for non-clumped particles. > > > > But some clump members disappear after a first step. Is it somewhere in > > clump.cpp? > There is no mysterious math in clump, there are conjugate quaternions, but > nothing special. Try creating clumps, save it and see if there is difference > in the XML for eigen and non-eigen versions. There is eigenvalue > decomposition used to compute clump inertia matrix and it might be not > correctly implemented in the glue. 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