Function factorial_array 

pub fn factorial_array<D>(
    input: &Array<u32, D>,
    config: &ArrayConfig,
) -> SpecialResult<Array<f64, D>>
where
    D: Dimension,

Enhanced factorial function computation