def get_enclosing_orthorhombic_box(cell):
<EOL>cx, cy, cz = cell<EOL>c1 = cx+cy<EOL>c2 = cx+cz<EOL>c3 = cy+cz<EOL>c4 = cx+cy+cz<EOL>corners = np.array([[<NUM_LIT:0>,<NUM_LIT:0>,<NUM_LIT:0>],cx,cy,cz,c1,c2,c3,c4])<EOL>lower = np.min(corners, axis=<NUM_LIT:0>)<EOL>upper = np.max(corners, axis=<NUM_LIT:0>)<EOL>return lower, upper<EOL>
Return lower and upper bounds of the orthorhombic box that encloses
the parallelepiped spanned by the three cell vectors of cell.