Hi John,

Thanks!  I *think* it works--here is my only confusion.  If I set c_max to some value less than the max density (say, 1e-24 in the code I attached previously), it doesn't crash anymore but it doesn't actually plot a clump contour either.  If I set c_max to the max density setting (or I guess a bit above, 10**np.floor(np.log10(data_source[field]).max()+1) ), then I get a contour plotted.  Does this seem right to you?  

Also, they both give me this message "/mnt/xfs1/home/stonnesen/yt-conda/lib/python3.6/site-packages/matplotlib/contour.py:1180: UserWarning: No contour levels w
ere found within the data range.
  warnings.warn("No contour levels were found"
"
Even though I seem to have one central contour in the c_max = 10**np.floor(np.log10(data_source[field]).max()+1) case.

Thanks!
Stephanie

--
Dr. Stephanie Tonnesen
Associate Research Scientist
CCA, Flatiron Institute
New York, NY

stonnes@gmail.com


On Tue, Oct 2, 2018 at 7:38 AM Wise, John H <jwise@physics.gatech.edu> wrote:
Hi Stephanie,

I actually ran into this bug in the past week.  I fixed in for my case
by making this change in clump_handling.py

-            elif (child._validate()):
+            elif ((child.children is not None) and child._validate()):

I haven't had the time to submit a PR. If this works for you, let me
know, and I can I create one this week.

Thanks,
John


On 10/01/2018 09:48 PM, Stephanie Tonnesen wrote:
> Hi,
>
> Sorry about the crazy delay, here is a script that reproduces the error
> with the IsolatedGalaxy dataset!
>
> The script is attached.  I am happy to help someone look into it, or
> look more myself with a bit of extra guidance.
>
> Thanks!!
>
> Stephanie
>
> --
> Dr. Stephanie Tonnesen
> Associate Research Scientist
> CCA, Flatiron Institute
> New York, NY
>
> stonnes@gmail.com <mailto:stonnes@gmail.com>
>
>
> On Fri, Sep 14, 2018 at 6:51 AM Stephanie Tonnesen <stonnes@gmail.com
> <mailto:stonnes@gmail.com>> wrote:
>
>     Hi Britton!
>
>     No, I have not resolved the issue, unfortunately.  I will try and
>     get a script that reproduces the problem up soon.
>
>     Best,
>     Stephanie
>
>     On Thu, Sep 13, 2018 at 7:45 PM Britton Smith
>     <brittonsmith@gmail.com <mailto:brittonsmith@gmail.com>> wrote:
>
>         HI Stephanie,
>
>         Have you managed to resolve this issue? If so, would you be
>         willing to submit a pull request to get this fix in the main
>         repo? If not, could you post a script that could be used to
>         reproduce this error with some of the public data?
>
>         Britton
>
>         On Fri, Sep 7, 2018 at 1:05 PM Nathan Goldbaum
>         <nathan12343@gmail.com <mailto:nathan12343@gmail.com>> wrote:
>
>             Generally an error like this is due to a bug in yt. I don't
>             know enough about the code you're using to provide useful
>             advice, but perhaps we should be checking if `children` is
>             `None` and doing something else. The "workaround" is likely
>             to fix whatever bug you've run into.
>
>             On Fri, Sep 7, 2018 at 2:56 PM, Stephanie Tonnesen
>             <stonnes@gmail.com <mailto:stonnes@gmail.com>> wrote:
>
>                 Hi yt-users!
>
>                 I have been fiddling some more with yt's clump finder
>                 and have another question.  I want to use density to
>                 find clumps, but want my maximum threshold value to be
>                 below the fiducial setting--instead of above the maximum
>                 possible density in the code I just want the finder to
>                 stop searching at, say 3.2e-24.
>
>                 The code was working but running out of memory until I
>                 set the max threshold, and I now get this error:
>
>                 ....lots of pleasant looking output messages....
>                 Updating joins ... :   0%|          | 0/7 [00:00<?, ?it/s]
>                 Updating joins ... : 100%|██████████| 7/7 [00:00<00:00,
>                 48933.55it/s]
>                 yt : [INFO     ] 2018-09-07 14:08:55,007 Investigating 2
>                 children.
>                 yt : [INFO     ] 2018-09-07 14:08:55,007 Finding clumps:
>                 min: 3.200000e-24, max: 3.200000e-24, step: 2.000000
>                 yt : [INFO     ] 2018-09-07 14:08:55,041 Finding clumps:
>                 min: 3.200000e-24, max: 3.200000e-24, step: 2.000000
>                 yt : [INFO     ] 2018-09-07 14:08:55,042 Eliminating
>                 invalid, childless clump with 2 cells.
>                 yt : [INFO     ] 2018-09-07 14:08:55,042 1 of 2 children
>                 survived, linking its children to parent.
>                 Traceback (most recent call last):
>                    File "yt_clump_fingert.py", line 59, in <module>
>                      find_clumps(master_clump, c_min, c_max, step)
>                    File
>                 "/mnt/home/stonnesen/yt-conda/yt/yt/analysis_modules/level_sets/clump_handling.py",
>                 line 444, in find_clumps
>                      find_clumps(child, min_val*d_clump, max_val, d_clump)
>                    File
>                 "/mnt/home/stonnesen/yt-conda/yt/yt/analysis_modules/level_sets/clump_handling.py",
>                 line 438, in find_clumps
>                      find_clumps(clump, min_val*d_clump, max_val, d_clump)
>                    File
>                 "/mnt/home/stonnesen/yt-conda/yt/yt/analysis_modules/level_sets/clump_handling.py",
>                 line 438, in find_clumps
>                      find_clumps(clump, min_val*d_clump, max_val, d_clump)
>                    File
>                 "/mnt/home/stonnesen/yt-conda/yt/yt/analysis_modules/level_sets/clump_handling.py",
>                 line 438, in find_clumps
>                      find_clumps(clump, min_val*d_clump, max_val, d_clump)
>                    File
>                 "/mnt/home/stonnesen/yt-conda/yt/yt/analysis_modules/level_sets/clump_handling.py",
>                 line 461, in find_clumps
>                      for child in clump.children:
>                 TypeError: 'NoneType' object is not iterable
>
>                 I tried to circumvent this error by changing my max to
>                 2e-24 (not a multiple of 2 of 1e-25), and got the same
>                 error when the min was 3.2e-24 and the max was 2e-24 (!).
>
>                 Is there any way I can get around this error and have
>                 the clumpfinder stop looking once it has clump edges at
>                 my maximum density?
>
>                 Thanks!
>                 Stephanie
>                 --
>                 Dr. Stephanie Tonnesen
>                 Associate Research Scientist
>                 CCA, Flatiron Institute
>                 New York, NY
>
>                 stonnes@gmail.com <mailto:stonnes@gmail.com>
>
>                 _______________________________________________
>                 yt-users mailing list -- yt-users@python.org
>                 <mailto:yt-users@python.org>
>                 To unsubscribe send an email to
>                 yt-users-leave@python.org <mailto:yt-users-leave@python.org>
>
>
>             _______________________________________________
>             yt-users mailing list -- yt-users@python.org
>             <mailto:yt-users@python.org>
>             To unsubscribe send an email to yt-users-leave@python.org
>             <mailto:yt-users-leave@python.org>
>
>         _______________________________________________
>         yt-users mailing list -- yt-users@python.org
>         <mailto:yt-users@python.org>
>         To unsubscribe send an email to yt-users-leave@python.org
>         <mailto:yt-users-leave@python.org>
>
>     --
>     --
>     Dr. Stephanie Tonnesen
>     Associate Research Scientist
>     CCA, Flatiron Institute
>     New York, NY
>
>     stonnes@gmail.com <mailto:stonnes@gmail.com>
>
>
>
> _______________________________________________
> yt-users mailing list -- yt-users@python.org
> To unsubscribe send an email to yt-users-leave@python.org
>

--
John Wise
Associate Professor of Physics
Center for Relativistic Astrophysics, Georgia Tech
http://cosmo.gatech.edu
_______________________________________________
yt-users mailing list -- yt-users@python.org
To unsubscribe send an email to yt-users-leave@python.org