Force computation works correctly, it natoms must still be pow of 2
and must not exceed 1024.
for pbc use for_each instead of loop