fn pv_ball(
(nx, ny, nz): (usize, usize, usize),
r: f64,
vs: [f64; 3],
) -> Vec<f64>Expand description
The fast QSM’s smoothing ball: radius r mm on the grid (−n/2 : n/2−1)·vs per axis (half
integers on an odd axis, as the dipole kernel), 1 where the voxel centre is within r − h
(h the half diagonal), and partial-volume weighted by 20³ sub-samples at offsets
(−0.475:0.05:0.475)·vs where it is within r + h; normalised to sum 1. Centred.