Commit e578b2f5 authored by Martin Reinecke's avatar Martin Reinecke

even more nitpicky fix

parent 9bd38fd8
......@@ -397,7 +397,7 @@ class ES_Kernel
ES_Kernel(size_t supp_, size_t nthreads)
: ES_Kernel(supp_, 2., nthreads){}
double operator()(double v) const { return (v*v>=1.) ? 0 : exp(beta*(sqrt(1.-v*v)-1.)); }
double operator()(double v) const { return (v*v>1.) ? 0 : exp(beta*(sqrt(1.-v*v)-1.)); }
/* Compute correction factors for the ES gridding kernel
This implementation follows eqs. (3.8) to (3.10) of Barnett et al. 2018 */
double corfac(double v) const
......@@ -965,7 +965,7 @@ template<typename T, typename T2=complex<T>> class Helper
for (size_t i=nexp; i<nvecs; ++i)
kernel[i]=0;
for (size_t i=0; i<nvecs; ++i)
kernel[i] = (kernel[i]*kernel[i]>=1) ? T(0) : exp(beta*(sqrt(T(1)-kernel[i]*kernel[i])-T(1)));
kernel[i] = (kernel[i]*kernel[i]>1) ? T(0) : exp(beta*(sqrt(T(1)-kernel[i]*kernel[i])-T(1)));
if ((iu0<bu0) || (iv0<bv0) || (iu0+supp>bu0+su) || (iv0+supp>bv0+sv))
{
if (grid_w) { dump(); fill(wbuf.begin(), wbuf.end(), T(0)); }
......
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