File: src/Interaction/Models/contactForce/linearCF.hpp (present in the v-1.0 tag and current main)
In the constructor that precomputes the damping coefficients, the normal
direction is:
etha_n[i] = -2.0*log(en[i])*sqrt(kn[i])/
sqrt(pow(log(en[i]),static_cast<real>(2.0))+ pow(Pi,static_cast<real>(2.0)));
but the tangential line has sqrt(kt[i]) inside the log():
etha_t[i] = -2.0*log( et[i]*sqrt(kt[i]) )/
sqrt(pow(log(et[i]),static_cast<real>(2.0))+ pow(Pi,static_cast<real>(2.0)));
This looks like a typo: for the restitution-derived linear dashpot the stiffness belongs outside the logarithm, as in the normal direction (the denominator already uses the correct log(et[i])).
File:
src/Interaction/Models/contactForce/linearCF.hpp(present in thev-1.0tag and currentmain)In the constructor that precomputes the damping coefficients, the normal
direction is:
but the tangential line has
sqrt(kt[i])inside thelog():This looks like a typo: for the restitution-derived linear dashpot the stiffness belongs outside the logarithm, as in the normal direction (the denominator already uses the correct
log(et[i])).