pub fn mul(a: &[[f32; 4]; 4], b: &[[f32; 4]; 4]) -> [[f32; 4]; 4]
Matrix product a * b (the transform b is applied first, then a).
a * b
b
a