def unit_cell_to_lattice_array(cryst):
a_, b_, c_, alpha, beta, gamma = cryst<EOL>r_alpha = np.deg2rad(alpha)<EOL>r_beta = np.deg2rad(beta)<EOL>r_gamma = np.deg2rad(gamma)<EOL>volume = a_ * b_ * c_ * (<EOL><NUM_LIT:1> - np.cos(r_alpha)**<NUM_LIT:2> - np.cos(r_beta)**<NUM_LIT:2> - np.cos(r_gamma)**<NUM_LIT:2> + <NUM_LIT:2> *<EOL>np.cos(r_alpha) * np.cos(r_beta) * np.cos(r_gamma))**<NUM_LIT:0.5><EOL>a_x = a_<EOL>a_y = b_ * np.cos(r_gamma)<EOL>a_z = c_ * np.cos(r_beta)<EOL>b_x = <NUM_LIT:0><EOL>b_y = b_ * np.sin(r_gamma)<EOL>b_z = c_ * (<EOL>np.cos(r_alpha) - np.cos(r_beta) * np.cos(r_gamma)) / np.sin(r_gamma)<EOL>c_x = <NUM_LIT:0><EOL>c_y = <NUM_LIT:0><EOL>c_z = volume / (a_ * b_ * np.sin(r_gamma))<EOL>lattice_array = np.array(<EOL>[[a_x, a_y, a_z], [b_x, b_y, b_z], [c_x, c_y, c_z]])<EOL>return lattice_array<EOL>
Return parallelpiped unit cell lattice matrix.