[petsc-users] Load imbalance in DMSwarm & DMDA simulations

[Miguel Molinos]

Dear all,

I am currently running a PIC-like implementation based on DMSwarm and DMDA. In my setup, a particle (DMSwarm) discretization represents, for example, a sphere placed at the center of a box-shaped domain discretized with a regular mesh (DMDA).

The implementation uses both MPI and OpenMP to accelerate the computations. The domain decomposition of the sphere follows a PIC-like approach: each particle is assigned to a rank based on the mesh element it belongs to, and thus inherits the rank that “owns” that element. This is essentially the same strategy used in one of the PETSc examples.

However, I have observed that, for certain configurations, some ranks end up with very few or even no particles, which leads to load imbalance.

Has anyone experienced a similar issue?

Thanks,
Miguel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20260317/35eb434b/attachment.html>