yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #22403
Re: [Question #689025]: Difference in unbalancedForce() between yade and yadedaily
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 <
question689025@xxxxxxxxxxxxxxxxxxxxx> 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 : yade-users@xxxxxxxxxxxxxxxxxxx
> 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.