Commit fffb6965 authored by Cristian Lalescu's avatar Cristian Lalescu
Browse files

fix generation of interpolator compute array

parent e8a0dbac
......@@ -51,7 +51,7 @@ rFFTW_interpolator<rnumber, interp_neighbours>::rFFTW_interpolator(
this->compute = new bool[this->descriptor->sizes[0]];
std::fill_n(this->compute, this->descriptor->sizes[0], false);
for (int iz = this->descriptor->starts[0]-interp_neighbours-1;
iz <= this->descriptor->starts[0]+this->descriptor->subsizes[0]+interp_neighbours+1;
iz <= this->descriptor->starts[0]+this->descriptor->subsizes[0]+interp_neighbours;
iz++)
this->compute[((iz + this->descriptor->sizes[0]) % this->descriptor->sizes[0])] = true;
}
......
......@@ -54,8 +54,8 @@ if __name__ == '__main__':
['-n', '32',
'--run',
'--initialize',
'--ncpu', '2',
'--nparticles', '1000',
'--ncpu', '4',
'--nparticles', '10000',
'--niter_todo', '32',
'--precision', 'single',
'--wd', 'data/single'] +
......
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