yade-dev team mailing list archive

[Václav Šmilauer <eudoxos@xxxxxxxx>]

> > The fix in cohesiveFrictionalContactLaw.cpp:133 is :
> > if (isnan(aa.angle())) aa.angle()=0;
> 
> OMG, guys. we need to either modify eigen, or write a wrapper for
> quaternions. It is not possible to put such stupid lines everywhere
> in the code. Same goes about quaternions losing their norm()==1.0

Go ask upstream eigen. They are very responsive on their forums and
bugtracker. It is NORMAL that quaternions get denormalized, I am
surprised you do not know that??!

v.