fn catmull(t: f32) -> [f32; 4]
Catmull-Rom weights for taps at offsets -1, 0, 1, 2 for fractional position t.
t