Jens
The "DISTRIB=..." options only determine on which process the particles are held. This is a purely computational property and is independent of any physics. You cannot change this without updating the time integration algorithm.
The particles are set up via some kind of initial condition, and there might be a run-time parameter that you can set to modify this.
-erik
On Wed, May 10, 2017 at 8:42 AM, Jens Mahlmann jensmahl@alumni.uv.es wrote:
Dear all,
I am currently working on the implementation of tracer particles into an evolution code. As for this, my orientation is the particle implementation of GRHydro. The number of particles is, in this case, is distributed equally across the processors by the DISTRIB=DEFAULT option. Is there any experience out there with distributing arrays according to a transfer function (e.g., some sort of density), such that higher density regions have more traced particles? In the documentation I only found the DISTRIB=CONSTANT option - which would not serve the purpose.
Thank you very much and my best!
Jens
Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users