/// Rotate an r-matrix about the z-axis.
pub fn rz(psi: f64, r: &mut [[f64; 3]; 3]) {
let s = psi.sin();
let c = psi.cos();
let a00 = c * r[0][0] + s * r[1][0];
let a01 = c * r[0][1] + s * r[1][1];
let a02 = c * r[0][2] + s * r[1][2];
let a10 = -s * r[0][0] + c * r[1][0];
let a11 = -s * r[0][1] + c * r[1][1];
let a12 = -s * r[0][2] + c * r[1][2];
r[0][0] = a00;
r[0][1] = a01;
r[0][2] = a02;
r[1][0] = a10;
r[1][1] = a11;
r[1][2] = a12;
}