pub fn compute_cosine_incidence_angle(
    body: &AirlessBody,
    normals: &Matrix3xX<Float>,
    sun_direction: &Vec3
) -> DRVector<Float>