Commit c70bd81d authored by Martin Reinecke's avatar Martin Reinecke
Browse files

fixes

parent f009cc0c
......@@ -873,6 +873,7 @@ template<typename T> void apply_holo
for (size_t ipart=0; ipart<nvis; ++ipart)
{
UVW<T> coord = baselines.effectiveCoord(idx(ipart));
coord.FixW();
hlp.prep(coord.u, coord.v);
complex<T> r = 0;
const auto * ptr = hlp.p0r;
......@@ -945,6 +946,7 @@ template<typename T> void get_correlations
for (size_t ipart=0; ipart<nvis; ++ipart)
{
UVW<T> coord = baselines.effectiveCoord(idx(ipart));
coord.FixW();
hlp.prep(coord.u, coord.v);
auto * wptr = hlp.p0w + u0*jump;
auto f0 = emb*emb;
......
......@@ -270,8 +270,8 @@ def test_wstack_against_wdft(nxdirty, nydirty, nchan, nrow, fov):
ms = np.random.rand(nrow, nchan)-0.5 + 1j*(np.random.rand(nrow, nchan)-0.5)
vis = bl.ms2vis(ms, idx)
uvw = bl.effectiveuvw(idx)
res0 = np.zeros((nxdirty, nydirty), dtype=np.complex128)
mi, ma = np.min(uvw[:, 2]), np.max(uvw[:, 2])
res0 = np.zeros((nxdirty, nydirty), dtype=np.float64)
mi, ma = np.min(np.abs(uvw[:, 2])), np.max(np.abs(uvw[:, 2]))
nplanes = 10000
ws = mi + np.arange(nplanes)*(ma-mi)/(nplanes-1)
for ii in range(len(ws)-1):
......
Markdown is supported
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