Function row_yz 

fn row_yz(m: &[[f32; 4]; 4], y: usize, z: usize) -> [f32; 3]

niimath’s per-row partial sums y·m[r][1] + z·m[r][2] + m[r][3], in f32.