Commit 11e0ccd5 authored by Philipp Arras's avatar Philipp Arras
Browse files

Increase accuracy of w term

parent ee287715
......@@ -304,10 +304,13 @@ def test_against_wdft(nxdirty, nydirty, epsilon, nchan, nrow):
x *= conf.Pixsize_x()
y *= conf.Pixsize_y()
res1 = np.zeros_like(res0)
n = np.cos(np.sqrt(x**2 + y**2))
eps = np.sqrt(x**2+y**2)
s = np.sin(eps)
nm1 = -s*s/(1+np.cos(eps))
n = nm1+1
for ii in idx:
phase = x*uvw[ii, 0] + y*uvw[ii, 1] + uvw[ii, 2]*(n-1)
print(np.max(x*uvw[ii, 0]), np.max(y*uvw[ii, 1]),np.max(uvw[ii, 2]*(n-1)) )
phase = x*uvw[ii, 0] + y*uvw[ii, 1] + uvw[ii, 2]*nm1
res1 += (vis[ii]*np.exp(2j*np.pi*phase)).real
res1 /= n
......
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