I had also this idea, that Eigen is better. It is here: http://bazaar.launchpad.net/%7Eyade-dev/yade/trunk/annotate/head%3A/pkg/dem/Engine/GlobalEngine/NewtonIntegrator.cpp#L183
Vector3r axis = state->angVel; Real angle = axis.Normalize(); ______________________________ Anton Gladkyy 2010/3/9 Václav Šmilauer <eudo...@arcig.cz> > > 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. > > Definitely Eigen's (since normalized vector should have |v|=1, right?). do > NOT add this to yadeMath, it is correct to fail for such nonsense as > normalizing 0 vector. Please document where it happens, that part of yade > has to be fixed. > > _______________________________________________ > 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
