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