);     }          return res; }  template < class VoxelType, unsigned int ndim > NumpyAnyArray  pythonGaussianGradientMagnitudeND(NumpyArray<ndim, Multiband<VoxelType> > volume,                                   const ConvolutionOptions<ndim-0> & opt,                                   NumpyArray<ndim