Question #689025 on Yade changed: https://answers.launchpad.net/yade/+question/689025
Chareyre proposed the following answer: My dear Sacha, thanks a lot for revealing two issues in a row. Do you know that you can report them on gitlab directly (gitlab issue = bug report)? I hope that one will be as fast as previous one but I have doubts. If you have candidate scripts that could help keeping clumps constantly working they are welcome, we apparently miss a few of them. You can check «checkScripts» in sources for examples. Bruno Le jeu. 27 févr. 2020 15:37, Sacha Duverger < [email protected]> a écrit : > New question #689025 on Yade: > https://answers.launchpad.net/yade/+question/689025 > > Hello, > > While I was performing triaxial tests I noticed that the same script gave > different values for unbalancedForce() depending on which version of yade I > was using: yade 2018.02b or yadedaily 20191105-2658~35899c7~bionic1. It > seems that the problem only occurs with clumps. > > Here is a MWE: > > ################## > > from yade import plot > import numpy as np > > O.bodies.appendClumped([sphere((0,0,0), 1), sphere((1.5,0,0), 1)]) > > wallIds = O.bodies.append(aabbWalls(thickness=1e-10)) > > O.engines=[ForceResetter() > ,InsertionSortCollider([Bo1_Sphere_Aabb(), > Bo1_Box_Aabb()]) > ,InteractionLoop( > [Ig2_Sphere_Sphere_ScGeom(), > Ig2_Box_Sphere_ScGeom6D()], > [Ip2_FrictMat_FrictMat_FrictPhys()], > [Law2_ScGeom_FrictPhys_CundallStrack()]) > ,GlobalStiffnessTimeStepper() > ,NewtonIntegrator(damping=0, gravity=(0,-10,0)) > > ,PyRunner(command='yade.plot.addData(uF=unbalancedForce())', iterPeriod=1) > ] > > O.run(1000, wait=True) > > print(np.mean(plot.data['uF'])) > > ################## > > yade 2018.02b outputs 0.11062286226025234 . > yadedaily 20191105-2658~35899c7~bionic1 outputs 75.5814592275 . > > Which version of Yade should I trust? > > Thanks in advance > > -- > 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 : [email protected] > Unsubscribe : https://launchpad.net/~yade-users > More help : https://help.launchpad.net/ListHelp > > > -- 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 : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
