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