fn square_to_quad(q: &[[f32; 2]; 4]) -> [[f32; 3]; 3]
Homography mapping the unit square (0,0)(1,0)(1,1)(0,1) onto q (Heckbert).
q