Function compare_intensity 

pub fn compare_intensity(
    reference: &[f64],
    candidate: &[f64],
    saturation_value: Option<f64>,
) -> Result<IntensityComparisonMetrics>