Function gamma_array 

pub fn gamma_array<D>(
    input: &Array<f64, D>,
    config: &ArrayConfig,
) -> SpecialResult<GammaResult<D>>
where
    D: Dimension + Send + Sync + 'static,

Enhanced gamma function computation with backend selection