pub fn multiply(a: [f32; 16], b: [f32; 16]) -> [f32; 16]
Column-major multiplication. Using matrices preserves shear from nonuniform parent scale followed by child rotation instead of approximating it as TRS.