yade-users team mailing list archive

["Sergei D." <sj2001@xxxxxxxxx>]

Hi, all.
Is anybody uses clumps for now?
Seems that clumps incorrectly implemented now.
For example, in NewtonDampedLaw.cpp line 28:
00028         const Vector3r& m=ncb->bex.getTorque(memberId); const
Vector3r& f=ncb->bex.getForce(memberId);

so m and f in global coord. But line 31:

00031         Vector3r
diffClumpAngularAccel=diagDiv(m,clumpRBP->inertia)+diagDiv((rb->se3.position-clumpRBP->se3.position).Cross(f),clumpRBP->inertia);

clumpRBP in local clumps coord.