void write_output (const VectorType& data, const vector<vector<int> >& mask_indices, ImageType& image) { for (size_t i = 0; i < mask_indices.size(); i++) { for (size_t dim = 0; dim < image.ndim(); dim++) image.index(dim) = mask_indices[i][dim]; image.value() = data[i]; } }