yade-dev team mailing list archive

[chiara modenese <c.modenese@xxxxxxxxx>]

  2 - applyForceAtContactPoint cannot be used in periodic + ScGeom, since
> this function needs the wrapped position of each body (in Dem3Dof, wrapped
> position is an
>  attribute, in ScGeom there is no such attribute and we don't want to
> define it as it costs cpu). The function call is replaced by separate
> applyForce/applyTorque using a
>  "branch" computed using radius, normal, and penetration depth.
>
> Hi Bruno,
also the function updateShearForce needs the positions of the particles
(which are accessed throughout the State pointer). Why would be working
there and not in applyForceAtContactPoint?
thanks, Chiara