spec_2d=fft.fftshift(tempamp)#1D spectrumleng = len(kstep) spec_1d = npy.zeros(leng) krange = Kh <= kstep[0] spec_1d[0] = spec_2d[krange].sum()forindinrange(1, leng): krange = (kstep[ind-1] < Kh) & (Kh <= kstep[ind]) spec_1d[ind] = spec_2d[krange].sum() ...
:type fourier_E: numpy 2D array, must be of type np.complex128 and must have a shape that is appropriate for a real fourier transform, i.e. (N,N/2 + 1); N should be a power of 2 :param fourier_B: B mode of the shear map in fourier space :type fourier_B: numpy 2D array,...