Commit 37304c94 authored by Martin Reinecke's avatar Martin Reinecke
Browse files

fixes

parent 9dbdef7a
...@@ -1462,7 +1462,7 @@ template<typename T> pyarr_c<uint32_t> getIndices(const Baselines<T> &baselines, ...@@ -1462,7 +1462,7 @@ template<typename T> pyarr_c<uint32_t> getIndices(const Baselines<T> &baselines,
return res; return res;
} }
template<typename T> pyarr_c<T> vis2dirty_wstack(const Baselines<T> &baselines, template<typename T> pyarr_c<complex<T>> vis2dirty_wstack(const Baselines<T> &baselines,
const GridderConfig<T> &gconf, const pyarr<uint32_t> &idx_, const GridderConfig<T> &gconf, const pyarr<uint32_t> &idx_,
const pyarr<complex<T>> &vis_) const pyarr<complex<T>> &vis_)
{ {
...@@ -1512,16 +1512,16 @@ cout << "nplanes: " << nplanes << endl; ...@@ -1512,16 +1512,16 @@ cout << "nplanes: " << nplanes << endl;
for (int i=max<int>(0,iplane-w_supp); i<min<int>(nplanes,iplane+w_supp+1); ++i) for (int i=max<int>(0,iplane-w_supp); i<min<int>(nplanes,iplane+w_supp+1); ++i)
if (abs(wval[ipart]-(wmin+i*dw))<dwmax) ++nvis_plane[i]; if (abs(wval[ipart]-(wmin+i*dw))<dwmax) ++nvis_plane[i];
} }
//cout << "nvis/plane: "; cout << "nvis/plane: ";
//for (auto nv:nvis_plane) cout << nv << " "; for (auto nv:nvis_plane) cout << nv << " ";
//cout << endl; cout << endl;
auto accum_ = makeArray<complex<T>>({nx_dirty, ny_dirty}); auto accum_ = makeArray<complex<T>>({nx_dirty, ny_dirty});
auto accum = accum_.mutable_data(); auto accum = accum_.mutable_data();
for (ptrdiff_t i=0; i<accum_.size(); ++i) for (ptrdiff_t i=0; i<accum_.size(); ++i)
accum[i] = 0.; accum[i] = 0.;
for (size_t iw=0; iw<nplanes; ++iw) for (size_t iw=0; iw<nplanes; ++iw)
{ {
//cout << "working on w plane #" << iw << endl; cout << "working on w plane #" << iw << endl;
if (nvis_plane[iw]==0) continue; if (nvis_plane[iw]==0) continue;
double wcur = wmin+iw*dw; double wcur = wmin+iw*dw;
size_t cnt=0; size_t cnt=0;
......
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