Second thoughts: Now I understand [1] a bit more but actually it solves a problem created ex-nihilo. There should not be any O.step() between appending bodies and setting unp=un, since obviously there is no reason to time-integrate before the initial conditions are completely defined. In this example time integration is used as a trick to detect interactions if I understand correctly. The trick is not needed, but I understand the idea. If you really want to do that you can just step with O.dt=0. I would not recommend it though, since it is less flexible than defining initial state completely in python script. Bruno
[1] https://answers.launchpad.net/yade/+question/266828 <https://answers.launchpad.net/yade/+question/266828>, #4 especially On 29 July 2015 at 03:51, Bruno Chareyre <bruno.chare...@3sr-grenoble.fr> wrote: > Hi Jérôme, > I don't see what you mean, sorry. > Setting unp=un is explicitely making current state the equilibrium state, > without the need of any additional "trick", else please explain why you > think it is not enough. > Handling modified normals (second specific feature) is the job of Ig2 > functor, not Law functor. > So overall I think the law you mention is really redundant. > > Let me know what you think about this. We should reach consensus about > such things, else users are all confused by different replies. > Thanks > > Bruno > > > On 28 July 2015 at 21:36, Jérôme Duriez < > question269...@answers.launchpad.net> wrote: > >> Question #269724 on Yade changed: >> https://answers.launchpad.net/yade/+question/269724 >> >> Jérôme Duriez posted a new comment: >> Yes Bruno, I know this "trick" too, which may require one or two more >> lines than the one you quote (see >> https://answers.launchpad.net/yade/+question/266828, #4 especially). >> >> Yade suffers probably from having too many laws, but as for >> Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM, I do not think this >> not so interesting discussion may conclude it is not justified, since I >> mentioned earlier a second specific feature about the discontinuity >> contacts..... >> >> -- >> You received this question notification because you are a member of >> yade-users, which is an answer contact for Yade. >> >> _______________________________________________ >> Mailing list: https://launchpad.net/~yade-users >> Post to : yade-us...@lists.launchpad.net >> Unsubscribe : https://launchpad.net/~yade-users >> More help : https://help.launchpad.net/ListHelp >> > > > _______________________________________________ > Mailing list: https://launchpad.net/~yade-dev > Post to : yade-dev@lists.launchpad.net > Unsubscribe : https://launchpad.net/~yade-dev > More help : https://help.launchpad.net/ListHelp > >
_______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : yade-dev@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp