Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

find_clump not working on octree-AMR datasets #4322

Open
PaoloSuin opened this issue Feb 2, 2023 · 1 comment
Open

find_clump not working on octree-AMR datasets #4322

PaoloSuin opened this issue Feb 2, 2023 · 1 comment

Comments

@PaoloSuin
Copy link

This is perhaps related to #2610, since other issues with volume rendering (#1212 (comment) )

    ad = gas.all_data()
    master_clump = Clump(ad, ("gas", "density"))
    c_min = ad["gas", "density"].min()
    c_max = ad["gas", "density"].max()
    step = 2.0
    find_clumps(master_clump, c_min, c_max, step)

leads to

Traceback (most recent call last):
  File "/home/psuin/.local/lib/python3.10/site-packages/IPython/core/interactiveshell.py", line 3378, in run_code
    exec(code_obj, self.user_global_ns, self.user_ns)
  File "<ipython-input-1-e086ea99c2eb>", line 1, in <module>
    find_clumps(master_clump, c_min, c_max, step)
  File "/home/psuin/.local/lib/python3.10/site-packages/yt/data_objects/level_sets/clump_handling.py", line 416, in find_clumps
    clump.find_children(min_val, max_val=max_val)
  File "/home/psuin/.local/lib/python3.10/site-packages/yt/data_objects/level_sets/clump_handling.py", line 152, in find_children
    nj, cids = identify_contours(self.data, self.field, min_val, max_val)
  File "/home/psuin/.local/lib/python3.10/site-packages/yt/data_objects/level_sets/contour_finder.py", line 23, in identify_contours
    for (g, node, (sl, dims, gi)) in data_source.tiles.slice_traverse():
  File "/home/psuin/.local/lib/python3.10/site-packages/yt/data_objects/data_containers.py", line 1508, in tiles
    self._tree = AMRKDTree(self.ds, data_source=self)
  File "/home/psuin/.local/lib/python3.10/site-packages/yt/utilities/amr_kdtree/amr_kdtree.py", line 176, in __init__
    raise RuntimeError(
RuntimeError: AMRKDTree does not support particle or octree-based data.

Is there any solution or workaround?

@welcome
Copy link

welcome bot commented Feb 2, 2023

Hi, and welcome to yt! Thanks for opening your first issue. We have an issue template that helps us to gather relevant information to help diagnosing and fixing the issue.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant