У меня есть следующий скрипт [code]#include #include
int main() { using namespace Eigen; Tensor xi1st(3, 3, 3, 3); for (int i = 0; i < 3; ++i) { for (int j = 0; j < 3; ++j) { for (int k = 0; k < 3; ++k) { for (int l = 0; l < 3; ++l) { xi1st(i, j, k, l) = std::complex(i + j + k + l, i - j + k - l); } } } } int i = 1, j = 1, k = 1; auto xi1st_slice = xi1st.slice(Eigen::array({0, i, j, k}), Eigen::array({3, 1, 1, 1})); auto xi1st_conj = xi1st_slice.conjugate(); std::complex xi1st_sum = xi1st_conj.sum(); std::cout