PNFFT: init f and grad_f with zeros for ik-differentiation