Commit 17928e9b authored by Martin Reinecke's avatar Martin Reinecke
Browse files

fixes

parent 37304c94
......@@ -1493,12 +1493,13 @@ double dw;
double x0 = -0.5*nx_dirty*psx,
y0 = -0.5*ny_dirty*psy;
double nmin = sqrt(1.-x0*x0-y0*y0)-1.;
dw = 1./(2*abs(nmin));
dw = 0.25/abs(nmin);
}
cout << "delta w: " << dw << endl;
double w_eps=1e-7; // FIXME
auto w_supp = get_supp(w_eps);
cout << "w_supp: " << w_supp << endl;
auto beta=2.3*w_supp;
wmin -= 0.5*w_supp*dw;
wmax += 0.5*w_supp*dw;
......@@ -1569,11 +1570,13 @@ cout << "applying correction for gridding in w direction" << endl;
{
double fy = y0+j*psy;
fy*=fy;
auto n=sqrt(1.-fx-fy);
auto n=sqrt(1.-fx-fy)-1.;
double fct = 0.;
for (int i=0; i<p; ++i)
fct += wgt[i]*psi[i]*cos(pi*w_supp*n*x[i]);
fct = 1./(w_supp*fct);
if ((i==nx_dirty/2) &&(j==ny_dirty/2))
cout << i << " " << j << " " << fct << endl;
size_t i2 = nx_dirty-i, j2 = ny_dirty-j;
accum[ny_dirty*i+j]*=fct;
if ((i>0)&&(i<i2))
......
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