<div dir="ltr"><div>I've reported this issue here: <a href="https://bitbucket.org/yt_analysis/yt/issues/1093/volumes-derived-from-chained-cut-regions">https://bitbucket.org/yt_analysis/yt/issues/1093/volumes-derived-from-chained-cut-regions</a></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 7, 2015 at 7:24 PM, Stephanie Tonnesen <span dir="ltr"><<a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Hi Nathan and all,<br><br></div>I have pasted a code that shows my problem using the isolated galaxy output. You will have to change the line of code that reads in the isolated galaxy output, but otherwise it should run for anyone. See the madness yourself!!! <br><br><a href="http://paste.yt-project.org/show/5943/" target="_blank">http://paste.yt-project.org/show/5943/</a><br><br></div>Thanks again,<br></div>Stephanie<br></div><div class="gmail_extra"><span class=""><br clear="all"><div><div><div dir="ltr"><div><div><div><font face="garamond,serif">--<br>Dr. Stephanie Tonnesen<br></font></div><font face="garamond,serif">Alvin E. Nashman Postdoctoral Fellow<br></font></div><font face="garamond,serif">Carnegie Observatories, Pasadena, CA<br></font></div><font face="garamond,serif"><a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a><br></font></div></div></div>
<br></span><div><div class="h5"><div class="gmail_quote">On Wed, Oct 7, 2015 at 4:33 PM, Stephanie Tonnesen <span dir="ltr"><<a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div>Hi Nathan,<br><br></div>Thanks for the quick response! I will see if I can reproduce this on a test dataset. <br><br></div>Different is the confusing part to me. These all should be selecting exactly the same region. If P3 and P4 are trying to select a region from ds.all_data() instead of ds.disk() then all three should be different. Another strange thing, if I select x > 65 kpc (only the top half of the cylinder, then the volume is 3.7384e-4, half of the values that don't seem right to me! <br><br></div>If I do <br></div>r4c = ds.cut_region(["obj['Metal_Density']/obj['density'] > 0.33"])<br></div>then the volume doesnt change, is 7.4768e-4<br><br></div>If I do:<br></div>r4 = ds.all_data()<br>r4c = r4.cut_region(["obj['Metal_Density']/obj['density'] > 0.33"])<br></div>then the volume is, as it should be, 0.3<br><br></div>*I think that something going wrong when I have more than one cut_region command.* If I do:<span><br><br>r3 = ds.disk([0.5,0.5,0.5],[-1,0,0],(4,"kpc"),(2,"kpc"))<br></span>r3c = r3.cut_region(["(obj['x'].in_units('kpc') > 60.0)"])<br></div><div> or<span><br>r2 = ds.disk([0.5,0.5,0.5],[-1,0,0],(4,"kpc"),(2,"kpc"))<br>r2c = r2.cut_region(["obj['Metal_Density']/obj['density'] > 0.33"])<br><br></span></div>then the volume is 9.346e-5, as it should be. <br><br></div>But when I put those together:<span><br><br>r3 = ds.disk([0.5,0.5,0.5],[-1,0,0],(4,"kpc"),(2,"kpc"))<br></span>r3r = r3.cut_region(["(obj['x'].in_units('kpc') > 60.0)"])<span><br>r3c = r3r.cut_region(["obj['Metal_Density']/obj['density'] > 0.33"])<br><br></span></div>the volume is 7.4768e-4 again!<br><br></div>And when I do this:<br>r4 = ds.all_data()<br>r4r = r4.cut_region(["obj['Metal_Density']/obj['density'] > 0.33"])<br></div>r4c = r4r.cut_region(["(obj['z'].in_units('kpc') > 0.0)"])<br><br></div>I get 2.4467, which is more than my box size. <br></div><div class="gmail_extra"><span><br clear="all"><div><div><div dir="ltr"><div><div><div><font face="garamond,serif">--<br>Dr. Stephanie Tonnesen<br></font></div><font face="garamond,serif">Alvin E. Nashman Postdoctoral Fellow<br></font></div><font face="garamond,serif">Carnegie Observatories, Pasadena, CA<br></font></div><font face="garamond,serif"><a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a><br></font></div></div></div>
<br></span><div><div><div class="gmail_quote">On Wed, Oct 7, 2015 at 3:30 PM, Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Stephanie,<div><br></div><div>I'm not sure off-hand why this is happening.</div><div><br></div><div>It may help one of us to be able to reproduce and understand what you're seeing if you can phrase your question in the form of an example that one of can run locally. This would make use of one of the test datasets on <a href="http://yt-project.org/data" target="_blank">yt-project.org/data</a>, and a full script that we can run locally pasted to <a href="http://paste.yt-project.org/data" target="_blank">paste.yt-project.org/data</a></div><div><br></div><div>Finally, just to clarify, the confusion is that the volumes of P3 and P4 are *bigger* than P2 and identical?</div><div><br></div><div>It might be the case that the derived quantity infrastructure is somehow selecting the data as if it were selecting from ds.all_data() instead of your disk selector. Can you try computing what the total volume of cells in a data object constructed by applying cut_regions on ds.all_data() instead of ds.disk()?</div><div><br></div><div>-Nathan</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Wed, Oct 7, 2015 at 5:21 PM, Stephanie Tonnesen <span dir="ltr"><<a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr"><div><div><div><div><div><div><div><div><div>Dear yt-users,<br><br></div>I am running<br><br>Version = 3.2-dev<br>Changeset = dc2467c4eae7 (yt) tip<br><br></div>I have been fiddling around with cut_region in a disk region and have gotten really confused. I was initially trying to see whether defining a disk region such that the cylinder height was on the x-axis changed how I need to make a cut in the cylinder height. Anyway, as I was fiddling with these cuts I found that if I made a "spatial" cut (rather than a different variable cut) I was getting different volumes. Let me show you:<br><br>small bit of background: the whole box is 130 kpc across, so 0.5 0.5 0.5 is 65 65 65 kpc. Also, I am cutting and pasting bits here so the indentation is off. <br><br>ds = yt.load("/data/stonnes/dwarf_nowind_lowmassgal_bigtracer/DD"+loop[i]+"/DD44sMBhn25now"+loop[i])<br><br>r2 = ds.disk([0.5,0.5,0.5],[-1,0,0],(4,"kpc"),(2,"kpc"))<br> r2c = r2.cut_region(["obj['Metal_Density']/obj['density'] > 0.33"])<br><br> r3 = ds.disk([0.5,0.5,0.5],[-1,0,0],(4,"kpc"),(2,"kpc"))<br> r3r = r3.cut_region(["(obj['x'].in_units('kpc') > 63.0)"])<br> r3r = r3r.cut_region(["(obj['x'].in_units('kpc') < 67.0)"])<br> r3c = r3r.cut_region(["obj['Metal_Density']/obj['density'] > 0.33"])<br> <br> r4 = ds.disk([0.5,0.5,0.5],[-1,0,0],(4,"kpc"),(2,"kpc"))<br> r4r = r4.cut_region(["(obj['x'].in_units('kpc') > 61.)"])<br> r4c = r4r.cut_region(["obj['Metal_Density']/obj['density'] > 0.33"])<br><br> P2 = np.append(P2,r2c.quantities.total_quantity([("index","cell_volume")]))<br> P3 = np.append(P3,r3c.quantities.total_quantity([("index","cell_volume")]))<br> P4 = np.append(P4,r4c.quantities.total_quantity([("index","cell_volume")]))<br><br></div>P2 volume is 9.346e-5 : within 2% of my calculation for the volume I should be getting<br></div>P3 volume is 7.4768e-4<br></div>P4 volume is 7.4768e-4<br><br></div>When I print out the x- y- and z- maxes and mins, they are the same for each of the selected regions. <br>(By doing this:<br>m2 = np.append(m2,r2c.quantities.extrema([("index","z")])[0])<br>rho2 = np.append(rho2,r2c.quantities.extrema([("index","z")])[1])<br><br>)<br><br></div>Can anyone help me figure out/tell me what is going on here?<br><br></div>Thanks!<br></div>Stephanie<br><br clear="all"><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div dir="ltr"><div><div><div><font face="garamond,serif">--<br>Dr. Stephanie Tonnesen<br></font></div><font face="garamond,serif">Alvin E. Nashman Postdoctoral Fellow<br></font></div><font face="garamond,serif">Carnegie Observatories, Pasadena, CA<br></font></div><font face="garamond,serif"><a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a><br></font></div></div></div>
</div></div></div></div></div></div></div></div></div></div></div></div></div>
<br></div></div>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>
<br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div></div></div>
</blockquote></div><br></div></div></div>
<br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>