Function discrete_gradient 

fn discrete_gradient(
    x: &[f64],
    grid: &Grid,
    axis: usize,
    transpose: bool,
) -> Vec<f64>

Allocating wrapper around discrete_gradient_into.