/// Reduce point modulo to the lattice.
inline gf_mesh<cyclic_lattice>::mesh_point_t gf_mesh<cyclic_lattice>::modulo_reduce(index_t const& r) const {
return mesh_point_t{*this, {_modulo(r[0], 0), _modulo(r[1], 1), _modulo(r[2], 2)}};
}
/// flatten the index
linear_index_t index_to_linear(index_t const& i) const {
return _modulo(i[0], 0) * s2 + _modulo(i[1], 1) * s1 + _modulo(i[2], 2);
}
/// Reduce index modulo to the lattice.
index_t index_modulo(index_t const& r) const {
return index_t{_modulo(r[0], 0), _modulo(r[1], 1), _modulo(r[2], 2)};
}