Function poisson_deviance 

pub fn poisson_deviance<T>(
    reference: ArrayView2<'_, T>,
    estimate: ArrayView2<'_, T>,
    valid_mask: Option<ArrayView2<'_, bool>>,
    epsilon: f64,
) -> Result<f64, IntensityMetricError>
where
    T: ToPrimitive,

Return summed Poisson deviance, flooring estimate intensities at epsilon.