Function psum 

fn psum<T>(on: Option<&Support>, n: usize, f: impl Fn(usize) -> T + Sync) -> T
where
    T: Copy + Send + Sync + Add<Output = T> + Default,

det_sum of f over a vector of n, or over its voxels in on (the rest zero)