Scalar finiteDifferenceIncrement() const { using math::sqrt; return 2.*sqrt(sqrt(Eigen::NumTraits<Scalar>::epsilon())); }