Function fgrad 

pub(crate) fn fgrad(
    x: &[f64],
    dims: (usize, usize, usize),
    vs: [f64; 3],
    g: &mut [Vec<f64>; 3],
)

fgrad.m: forward differences, zero at the last voxel along each axis