Commit 8ee92fa7 authored by alexander.dornheim's avatar alexander.dornheim
Browse files

Replaced sqrt by np.sqrt and more

parent dd9d6c2a
......@@ -109,20 +109,20 @@ def Goettingen_data_processor(config):
dqx = 2*np.pi/(Nx*dx)
dqy = 2*np.pi/(Ny*dy)
[Qx,Qy] = np.meshgrid(dqx*(range(Nx)-floor(Nx/2)-1),dqy*(range(Ny)-floor(Ny/2)-1))
[Qx,Qy] = np.meshgrid(dqx*(range(Nx)-np.floor(Nx/2)-1),dqy*(range(Ny)-np.floor(Ny/2)-1))
##################################
# Prepare data from ProxToolbox:
##################################
# Fresnel propagator:
config['use_farfield_formula']=false
config['use_farfield_formula']=False
config['Nx']=Nx
config['Ny']=Ny
config['fresnel_nr'] = 1*2*np.pi*Nx # Fresnel number
config['magn']=1 # magnification factor (should this be M above, or
# is it okay since the demagnified pixel size is used?)
kappa = sqrt(k**2-(Qx**2+Qy**2))
kappa = np.sqrt(k**2-(Qx**2+Qy**2))
config['FT_conv_kernel'] = fftshift(exp(-1j*kappa*z_eff))
config['beam'] = abs(ifft2(fft2(sqrt(WG*7.210116465530644e-04))/config['FT_conv_kernel']))
# the scaling factor of the waveguide array above was
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment