> > 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. _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
