diff --git a/nifty_gridder.cc b/nifty_gridder.cc index e449368295a321fcb5fac70039d469d6729647c5..e1c8b63b1aff8ed528e575d52360fdfc053be82d 100644 --- a/nifty_gridder.cc +++ b/nifty_gridder.cc @@ -1492,7 +1492,7 @@ double psz=gconf.Pixsize_x(); double dw = psz;// distance between two w planes, FIXME cout << "delta w: " << dw << endl; - double w_eps=1e-7; + double w_eps=1e-7; // FIXME auto w_supp = get_supp(w_eps); auto beta=2.3*w_supp; wmin -= 0.5*w_supp*dw; @@ -1507,15 +1507,16 @@ cout << "nplanes: " << nplanes << endl; 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]; } -cout << "nvis/plane: "; -for (auto nv:nvis_plane) cout << nv << " "; -cout << endl; +//cout << "nvis/plane: "; +//for (auto nv:nvis_plane) cout << nv << " "; +//cout << endl; auto accum_ = makeArray<complex<T>>({nx_dirty, ny_dirty}); auto accum = accum_.mutable_data(); for (ptrdiff_t i=0; i<accum_.size(); ++i) accum[i] = 0.; for (size_t iw=0; iw<nplanes; ++iw) { +//cout << "working on w plane #" << iw << endl; if (nvis_plane[iw]==0) continue; double wcur = wmin+iw*dw; size_t cnt=0; @@ -1540,6 +1541,7 @@ cout << endl; // correct for w gridding { py::gil_scoped_release release; +cout << "applying correction for gridding in w direction" << endl; constexpr double pi = 3.141592653589793238462643383279502884197; double x0 = -0.5*nx_dirty*psx,