ElementInertia.h

#pragma once
 
#include <Eigen/Eigen>
 
namespace RbCable
{
inline Eigen::Matrix3d TranslationalInverseMass(double mass, double addedMass, const Eigen::Vector3d& axis)
{
   const Eigen::Matrix3d along = axis * axis.transpose();
   return along / mass + (Eigen::Matrix3d::Identity() - along) / (mass + addedMass);
}
} // namespace RbCable