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.