Function sh_basis 

pub fn sh_basis(dir: Vec3) -> [f32; 9]

Evaluate all 9 SH basis functions at direction dir.

The returned array follows the convention: [Y0_0, Y1_{-1}, Y1_0, Y1_1, Y2_{-2}, Y2_{-1}, Y2_0, Y2_1, Y2_2]