Hi,
I’m running the simple_example.py script using the default settings, and everything works as expected.
However, when I change the ngrid value to something different from nparticles (for example, ngrid=256 and nparticles=128), the results become unphysical, there's almost no structure formation, and most particles appear clumped together in a small region.
I would like to use a finer grid size to improve the accuracy, like in COLA documentation suggests a factor 3. Do you know what might be causing this issue when the grid resolution is different than nparticles? Could it be related to the mass assignment scheme or the way forces are computed?
Thanks in advance,
Hi,
I’m running the
simple_example.pyscript using the default settings, and everything works as expected.However, when I change the
ngridvalue to something different fromnparticles(for example,ngrid=256andnparticles=128), the results become unphysical, there's almost no structure formation, and most particles appear clumped together in a small region.I would like to use a finer grid size to improve the accuracy, like in COLA documentation suggests a factor 3. Do you know what might be causing this issue when the grid resolution is different than nparticles? Could it be related to the mass assignment scheme or the way forces are computed?
Thanks in advance,