Commit 4893a233 authored by Martin Reinecke's avatar Martin Reinecke
Browse files

cleanup

parent 004be60c
......@@ -204,10 +204,8 @@ class Helper
iu0 = int(u-w*0.5 + 1 + nu) - nu;
if (iu0+w>nu+nsafe) iu0 = nu+nsafe-w;
auto v = fmodulo(v_in, 1.)*nv;
// cout << "uv0 " << u << " " << v << endl;
iv0 = int(v-w*0.5 + 1 + nv) - nv;
if (iv0+w>nv+nsafe) iv0 = nv+nsafe-w;
// cout << "iuv0 " << iu0 << " " << iv0 << endl;
double xw=2./w;
for (int i=0; i<w; ++i)
{
......@@ -215,10 +213,7 @@ class Helper
kernel[i+w] = xw*(iv0+i-v);
}
for (int i=0; i<2*w; ++i)
{
kernel[i] = exp(beta*(sqrt(1-kernel[i]*kernel[i])-1.));
// cout << kernel[i] << endl;
}
}
bool need_to_move() const
......@@ -228,7 +223,6 @@ class Helper
{
bu0=max(-nsafe, min(nu+nsafe-su, iu0+nsafe-su/2));
bv0=max(-nsafe, min(nv+nsafe-sv, iv0+nsafe-sv/2));
//cout << "buv0 " << bu0 << " " << bv0 << endl;
}
protected:
......@@ -240,9 +234,6 @@ class Helper
bu0(-1000000), bv0(-1000000)
{
if (min(nu,nv)<2*nsafe) throw runtime_error("grid dimensions too small");
// cout << "beta=" << beta << endl;
// cout << "w=" << w << endl;
// cout << "suv=" << su << " " << sv << endl;
}
};
......
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