pub fn phase_hg(cos_theta: Float, g: Float) -> Float