yade-users team mailing list archive

[Jan Stránský <question406814@xxxxxxxxxxxxxxxxxxxxx>]

Question #406814 on Yade changed:
https://answers.launchpad.net/yade/+question/406814

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hello,
you can check the source code of kineticEnergy() [1] and as
O.energy['kinetic'] [2] is couputed.
The formulas used is the same (after a quick look, you may do more detailed
investigation if you are interested), however the computation itself **IS**
different, so the finite precision of computation and different order of
operations may cause different results [3].
How many particles do you have for this computation?

cheers
Jan

[1] https://github.com/yade/trunk/blob/master/pkg/dem/Shop_01.cpp#L180
[2]
https://github.com/yade/trunk/blob/master/pkg/dem/NewtonIntegrator.cpp#L40
[3] https://yade-dem.org/doc/formulation.html#result-indeterminism

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.