pub fn compute_cosine_phase_angle(
    body: &AirlessBody,
    camera_direction: &Vec3,
    sun_direction: &Vec3
) -> DRVector<Float>