docs: def get_enclosing_orthorhombic_box(cell):
      """
      Return lower and upper bounds of the orthorhombic box that encloses
      the parallelepiped spanned by the three cell vectors of cell.
      """
      
      # Cell vectors
      cx, cy, cz = cell
      
      # The cell has eight corners, one is at the origin, three at cx, cy, cz
      # and the last ones are...
      c1 = cx+cy
      c2 = cx+cz
      c3 = cy+cz
      c4 = cx+cy+cz
      
      corners = np.array([[0,0,0],cx,cy,cz,c1,c2,c3,c4])
      lower = np.min(corners, axis=0)
      upper = np.max(corners, axis=0)
      
      return lower, upper