virtual double bondFunction(const Vector3d &vec1, const Vector3d &vec2) const =0
virtual double operator()(const Vector3d &vec) const =0
MethodUser(double averagingDomainSize)
virtual ~MethodUser()=default
double averagingDomainSize
Eigen::Matrix< double, 1, DIM, Eigen::RowMajor > Vector3d