Эффективное суммирование функции по поводу целочисленных смены решетки в собственном собственном (C ++ 23) ⇐ C++

[Anonymous]

У меня есть функция p Возвращение real_type и взяв собственное :: vector . Математически мне нужно вычислить сумму p (x + k) для каждого D-мерного целочисленного вектора k , компоненты которых находятся внутри [-K_MAX, k_max] . (На практике k_max очень мал, как 3 или 4). Мне нужно вычислить эту сумму миллион раз для различного x (но с тем же k_max ). Как наиболее показан способ достичь этого?template
auto sum_p(P p, Eigen::Vector const& x, std::size_t const k_max = 10)
{
// sum p(x + k) over all $k\in\mathbb Z^D$ with $|k_i|\le k_max$
auto acc = p(x);
for (std::size_t k = 1; k

Подробнее здесь: https://stackoverflow.com/questions/796 ... -eigen-c23