yade-users team mailing list archive

[Bokkisa Srinivas Vivek <question680398@xxxxxxxxxxxxxxxxxxxxx>]

New question #680398 on Yade:
https://answers.launchpad.net/yade/+question/680398

Hello community,

I was going through the Source code- Clump.cpp

The inertia tensor at local level for the intersecting clump is calculated from the below equation (line no. 178):

Ig += m*( x.dot(x)*Matrix3r::Identity()-x*x.transpose())+Matrix3r(Vector3r::Constant(dv*pow(dx,2)/6.).asDiagonal())

In the second part- "Matrix3r(Vector3r::Constant(dv*pow(dx,2)/6.).asDiagonal())", 

I believe it should be m*pow(dx,2)/6 but instead it computes dv*pow(dx,2)/6, was it a mistake in code or done purposefully?

Reference:
https://github.com/yade/trunk/blob/master/core/Clump.cpp#L178

Thank you





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