Annotate Clumps Not Working in Odd MPI Size

[cindytsai]

Annotate Clumps Not Working in Odd MPI Size

• Test Problem: gamer MHD Vortex
• Inline Script:

  import numpy as np
  import yt
  from yt.data_objects.level_sets.api import Clump, find_clumps
  yt.enable_parallelism()
  def yt_inline():
      ds = yt.frontends.libyt.libytDataset()
      data_source = ds.all_data()
  
      c_min = 10 ** np.floor(np.log10(data_source[("gas", "density")]).min())
      c_max = 10 ** np.floor(np.log10(data_source[("gas", "density")]).max() + 1)
  
      master_clump = Clump(data_source, ("gas", "density"))
      master_clump.add_validator("min_cells", 20)
  
      find_clumps(master_clump, c_min, c_max, 2.0)
      leaf_clumps = master_clump.leaves
  
      prj = yt.ProjectionPlot(ds, "z", ("gamer","Dens"), center="c")
      prj.annotate_clumps(leaf_clumps)
  
      # Either having this is clause or not, it will still failed when MPI = 3.
      if yt.is_root():
          prj.save()

• Description:
  • FAILED in MPI = 3, stuck at somewhere else, other than class IOHandlerlibyt.
  • SUCCEED in MPI = 1, and MPI = 4.
    • Just like the projection plot, because i didn't set a range that actually grab a clump.