def xx2k(self, xx):
"""
Private: oversampled FFT on CPU
Firstly, zeroing the self.k_Kd array
Second, copy self.x_Nd array to self.k_Kd array by cSelect
Third: inplace FFT
"""
# dd = numpy.size(self.Kd)
output_x = numpy.zeros(self.Kd, dtype=self.dtype,order='C')
# output_x[crop_slice_ind(xx.shape)] = xx
output_x.flat[self.KdCPUorder]=xx.flat[self.NdCPUorder]
k = numpy.fft.fftn(output_x, self.Kd, range(0, self.ndims))
return k
评论列表
文章目录