Question #697363 on Yade changed: https://answers.launchpad.net/yade/+question/697363
Status: Open => Answered Jan Stránský proposed the following answer: Hello, > 1 + 3 still, please be more specific. What is "the exterior"? Particles forming a convex hull? Particles close enough to the trimming sphere surface? ... ? Is more important center of particles? Or their surface? ... ? The solution is as easy as: ### exterior = [b for b in O.bodies if isExterior(b)] for b in exterior: O.forces.setPermF(b.id,someForce) # or addF ### isExterior is a function depending on your needs. E.g. the most primitive could be: ### limitDistance = 1.23 center = Vector3(1,2,3) def distance(v1,v2): return (v1-v2).norm() def isExterior(b): return distance(b.state.pos,center) < limitDistance ### > 2... efficiency w.r.t. code writing and efficiency w.r.t. performance are often contradicting requirements... Cheers Jan -- You received this question notification because your team yade-users is an answer contact for Yade. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp