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?
updateShearForce doesn't need positions, only velocities, except if
avoidGranularRatcheting=false.
I'll perhaps put a warning telling that you cannot use periodic BCs and
request granular ratcheting at the same time.
Considering that avoidGranularRatcheting=false leads to wrong physics
(energy created in elastic behaviour), it is really not a big deal...
Bruno
_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-dev
More help : https://help.launchpad.net/ListHelp
