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

Reply via email to