Support Particle Functionalities

[cindytsai]

Support Particle Functionalities

Initially, particle plots may generate false figure, if the MPI size is too large. After testing on different machine, this doesn't seems to be an issue in libyt. This is more or less related to memory space of one machine has. But we still cannot find where actually is this issue, we move to another issue.

• 41

:information_source: After testing on calab912 and eureka, it doesn't seem to be an issue. Instead, it is more or less related to memory on the machine. Though we still don't know why and where is the bug.

:heavy_check_mark: Particle plot (ParticlePlot, ParticleProjectionPlot) can successfully run in parallel in inline-analysis. Though some memory related issues may occurred.

Test Run on My Laptop with RAM = 16

All of the images are results from inline-analysis. MPI=1 is the correct and expected outcome.

• Test Problem: gamer Plummer
• Machine: My Laptop

  ParticlePlot

  import yt
  yt.enable_parallelism()
  def yt_inline():
  ds = yt.frontends.libyt.libytDataset()
  par = yt.ParticlePlot( ds, 'particle_position_x', 'particle_position_y', "particle_mass", center='c' )
  if yt.is_root():
  par.save()

• MPI=1
• MPI=2

ParticleProjectionPlot

import yt
yt.enable_parallelism()
def yt_inline():
    ds = yt.frontends.libyt.libytDataset()
    par = yt.ParticleProjectionPlot( ds, "z")
    if yt.is_root():
        par.save()

• MPI=1
• MPI=2
  • The upper-right cluster is different.

[cindytsai]

ParticlePlot

• Test Problem: gamer Plummer
• Machine: calab912
• Inline Script

  import yt
  yt.enable_parallelism()
  def yt_inline():
  ds = yt.frontends.libyt.libytDataset()
  par = yt.ParticlePlot( ds, 'particle_position_x', 'particle_position_y', "particle_mass", center='c' )
  if yt.is_root():
      par.save()

• Output: For MPI = 1, 2, 3, 4, they all generate the same figure. I think this is cause by memory... The previous tests are done on my laptop.

[cindytsai]

ParticlePlot

• Test Problem: gamer Plummer
• Machine: eureka
• MPI Version: 3.1.5-intel-oneapi
• Inline Script:

  import yt
  yt.enable_parallelism()
  def yt_inline():
    ds = yt.frontends.libyt.libytDataset()
    par = yt.ParticlePlot( ds, 'particle_position_x', 'particle_position_y', "particle_mass", center='c' )
    if yt.is_root():
        par.save()

• Run Command:

  OMPI_MCA_osc=sm,pt2pt mpirun -np 4 --output-filename log ./gamer

  • Should put OMPI_MCA_osc=sm.pt2pt at the front, otherwise will reach an assertion error when encountered RMA operations. This only happens in eureka.
• Output: MPI = 1 and MPI = 16 are the same.
  • MPI = 1
  • MPI = 16

[cindytsai]

ParticleProjectionPlot

• Test Problem: gamer Plummer
• Machine: eureka
• MPI Version: 3.1.5-intel-oneapi
• Inline Script:

  import yt
  yt.enable_parallelism()
  def yt_inline():
    ds = yt.frontends.libyt.libytDataset()
    par_prj = yt.ParticleProjectionPlot( ds, "z" )
    if yt.is_root():
        par_prj.save()

• Run Command:

  OMPI_MCA_osc=sm,pt2pt mpirun -np 32 --output-filename log ./gamer

• Output: MPI = 1 and MPI = 32 are the same.

  • MPI = 1

  • MPI = 32