const array<bool, dim> is_periodic{
      periodic || periodic_x,
      periodic || periodic_y,
      periodic || periodic_z,
  };

  const array<array<bool, 3>, 2> is_symmetry{{
      {{
          periodic || periodic_x || reflection_x,
          periodic || periodic_y || reflection_y,
          periodic || periodic_z || reflection_z,
      }},
      {{
          periodic || periodic_x || reflection_upper_x,
          periodic || periodic_y || reflection_upper_y,
          periodic || periodic_z || reflection_upper_z,
      }},
  }};

          vbx[orient(d, f)] == domain[orient(d, f)] && !is_periodic[d];

          vbx[orient(d, f)] == domain[orient(d, f)] && !is_symmetry[f][d];