initialized from Vega matrices. The ode rhs
now considers damping only if the coefficients are
positive, non zero. It also accumulates contact forces
if any. The updateAttachedMesh in the VegaVolumetricMesh
class now takes an Eigen vector to do the update.