pub fn vec3_reflect(v: [f64; 3], n: [f64; 3]) -> [f64; 3]
Reflect vector v about unit normal n: v - 2*(v·n)*n.
v
n
v - 2*(v·n)*n