pub fn bilinear(tl: f64, tr: f64, bl: f64, br: f64, tx: f64, ty: f64) -> f64
Bilinear interpolation on a unit square. tl=top-left, tr=top-right, bl=bottom-left, br=bottom-right. tx, ty in [0,1].