Function cool 

pub fn cool(temperature: f32, cooling_rate: f32) -> f32

Geometric cooling schedule T_{k+1} = T_k · cooling_rate (cooling_rate ∈ (0,1)).