pub(super) fn dconv( x: &[f64], d: &[Complex64], fft: &mut Fft3dWorkspace, buf: &mut [Complex64], out: &mut [f64], )
Dconv: out = real(ifftn(D.*fftn(x)))
Dconv
out = real(ifftn(D.*fftn(x)))