/// Rotate an r-matrix about the x-axis.
pub fn rx(phi: f64, r: &mut [[f64; 3]; 3]) {
let s = phi.sin();
let c = phi.cos();
let a10 = c * r[1][0] + s * r[2][0];
let a11 = c * r[1][1] + s * r[2][1];
let a12 = c * r[1][2] + s * r[2][2];
let a20 = -s * r[1][0] + c * r[2][0];
let a21 = -s * r[1][1] + c * r[2][1];
let a22 = -s * r[1][2] + c * r[2][2];
r[1][0] = a10;
r[1][1] = a11;
r[1][2] = a12;
r[2][0] = a20;
r[2][1] = a21;
r[2][2] = a22;
}