Commit 6bc37688 authored by Martin Reinecke's avatar Martin Reinecke
Browse files

fixes

parent 7575f488
......@@ -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,
......
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