Investigated compile error with @reister. Changing the type of the variable whose entries are used as indices to Eigen::Vector3i solves the problem. Using Eigen::Matrix<size_t, 6, 1>, which seemed even more suitable, did not work.