<meta http-equiv="Content-Type" content="text/html; charset=utf-8">Ah yes, the "ds" attribute of a data object is a reference to the dataset itself. Sorry that wasn't clear initially.<span></span><div><br>On Saturday, May 2, 2015, Aaron Smith <<a href="mailto:asmith@astro.as.utexas.edu">asmith@astro.as.utexas.edu</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
Hi all,<br>
<br>
Joe decided (offlist) to directly access the region object itself.
Because re.ds references the dataset that the region came from,
instead of using
<br>
<ol style="margin:0px;padding:0px 0px 0px 48px;color:rgb(172,172,172);font-family:Consolas,Menlo,Monaco,'Lucida Console','Liberation Mono','DejaVu Sans Mono','Bitstream Vera Sans Mono',monospace,serif;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:21px;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(248,248,248)">
<li>
<div>ds<span> </span><span style="color:rgb(102,204,102)">=</span><span> </span>load<span style="color:black">(</span><span style="color:rgb(72,61,139)">'myfile.hdf5'</span><span style="color:rgb(102,204,102)">,</span>over_refine_factor<span style="color:rgb(102,204,102)">=</span><span style="color:rgb(255,69,0)">1</span><span style="color:black">)</span></div>
</li>
<li>
<div><span style="color:rgb(220,20,60)">re</span><span> </span><span style="color:rgb(102,204,102)">=</span><span> </span>ds.<span style="color:black">region</span><span style="color:black">(</span>center<span style="color:rgb(102,204,102)">=</span><span style="color:black">[</span><span style="color:rgb(255,69,0)">750</span><span style="color:black">]</span>*<span style="color:rgb(255,69,0)">3</span><span style="color:rgb(102,204,102)">,</span><span> </span>left_edge<span style="color:rgb(102,204,102)">=</span><span style="color:black">[</span><span style="color:rgb(255,69,0)">700</span><span style="color:black">]</span>*<span style="color:rgb(255,69,0)">3</span><span style="color:rgb(102,204,102)">,</span><span> </span>right_edge<span style="color:rgb(102,204,102)">=</span><span style="color:black">[</span><span style="color:rgb(255,69,0)">800</span><span style="color:black">]</span>*<span style="color:rgb(255,69,0)">3</span><span style="color:black">)</span></div>
</li>
<li>
<div>ad<span> </span><span style="color:rgb(102,204,102)">=</span><span> </span><span style="color:rgb(220,20,60)">re</span>.<span style="color:black">ds</span>.<span style="color:black">all_data</span><span style="color:black">(</span><span style="color:black">)</span></div>
</li>
<li>
<div>density<span> </span><span style="color:rgb(102,204,102)">=</span><span> </span>ad<span style="color:black">[</span><span style="color:black">(</span><span style="color:rgb(72,61,139)">"PartType0"</span><span style="color:rgb(102,204,102)">,</span><span style="color:rgb(72,61,139)">"density"</span><span style="color:black">)</span><span style="color:black">]</span></div>
</li>
<li>
<div>wdens<span> </span><span style="color:rgb(102,204,102)">=</span><span> </span>np.<span style="color:black">where</span><span style="color:black">(</span>density<span> </span><span style="color:rgb(102,204,102)">==</span><span> </span>np.<span style="color:rgb(0,128,0)">max</span><span style="color:black">(</span>density<span style="color:black">)</span><span style="color:black">)</span></div>
</li>
<li>
<div>coordinates<span> </span><span style="color:rgb(102,204,102)">=</span><span> </span>ad<span style="color:black">[</span><span style="color:black">(</span><span style="color:rgb(72,61,139)">"PartType0"</span><span style="color:rgb(102,204,102)">,</span><span style="color:rgb(72,61,139)">"Coordinates"</span><span style="color:black">)</span><span style="color:black">]</span></div>
</li>
<li>
<div>cen<span> </span><span style="color:rgb(102,204,102)">=</span><span> </span>coordinates<span style="color:black">[</span>wdens<span style="color:black">]</span><span style="color:black">[</span><span style="color:rgb(255,69,0)">0</span><span style="color:black">]</span></div>
</li>
</ol>
the question can be resolved with<br>
<ol style="margin:0px;padding:0px 0px 0px 48px;color:rgb(172,172,172);font-family:Consolas,Menlo,Monaco,'Lucida Console','Liberation Mono','DejaVu Sans Mono','Bitstream Vera Sans Mono',monospace,serif;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:21px;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(248,248,248)">
<li>
<div>ds<span> </span><span style="color:rgb(102,204,102)">=</span><span> </span>load<span style="color:black">(</span><span style="color:rgb(72,61,139)">'myfile.hdf5'</span><span style="color:rgb(102,204,102)">,</span>over_refine_factor<span style="color:rgb(102,204,102)">=</span><span style="color:rgb(255,69,0)">1</span><span style="color:black">)</span></div>
</li>
<li>
<div><span style="color:rgb(220,20,60)">re</span><span> </span><span style="color:rgb(102,204,102)">=</span><span> </span>ds.<span style="color:black">region</span><span style="color:black">(</span>center<span style="color:rgb(102,204,102)">=</span><span style="color:black">[</span><span style="color:rgb(255,69,0)">750</span><span style="color:black">]</span>*<span style="color:rgb(255,69,0)">3</span><span style="color:rgb(102,204,102)">,</span><span> </span>left_edge<span style="color:rgb(102,204,102)">=</span><span style="color:black">[</span><span style="color:rgb(255,69,0)">700</span><span style="color:black">]</span>*<span style="color:rgb(255,69,0)">3</span><span style="color:rgb(102,204,102)">,</span><span> </span>right_edge<span style="color:rgb(102,204,102)">=</span><span style="color:black">[</span><span style="color:rgb(255,69,0)">800</span><span style="color:black">]</span>*<span style="color:rgb(255,69,0)">3</span><span style="color:black">)</span></div>
</li>
<li>
<div>density<span> </span><span style="color:rgb(102,204,102)">=</span><span> </span>re<span style="color:black">[</span><span style="color:black">(</span><span style="color:rgb(72,61,139)">"PartType0"</span><span style="color:rgb(102,204,102)">,</span><span style="color:rgb(72,61,139)">"density"</span><span style="color:black">)</span><span style="color:black">]</span></div>
</li>
<li>
<div>wdens<span> </span><span style="color:rgb(102,204,102)">=</span><span> </span>np.<span style="color:black">where</span><span style="color:black">(</span>density<span> </span><span style="color:rgb(102,204,102)">==</span><span> </span>np.<span style="color:rgb(0,128,0)">max</span><span style="color:black">(</span>density<span style="color:black">)</span><span style="color:black">)</span></div>
</li>
<li>
<div>coordinates<span> </span><span style="color:rgb(102,204,102)">=</span><span> </span>re<span style="color:black">[</span><span style="color:black">(</span><span style="color:rgb(72,61,139)">"PartType0"</span><span style="color:rgb(102,204,102)">,</span><span style="color:rgb(72,61,139)">"Coordinates"</span><span style="color:black">)</span><span style="color:black">]</span></div>
</li>
<li>
<div>cen<span> </span><span style="color:rgb(102,204,102)">=</span><span> </span>coordinates<span style="color:black">[</span>wdens<span style="color:black">]</span><span style="color:black">[</span><span style="color:rgb(255,69,0)">0</span><span style="color:black">]</span></div>
</li>
</ol>
<br>
Cheers,<br>
Aaron<br>
<pre cols="72">--
Aaron Smith
NSF Graduate Research Fellow
Department of Astronomy
University of Texas at Austin
<a href="http://www.as.utexas.edu/~asmith" target="_blank">www.as.utexas.edu/~asmith</a></pre>
<div><br>
On 5/2/15 5:27 PM, Joseph Smidt wrote:<br>
</div>
<blockquote type="cite">
<pre>Hey everyone,
I am trying to find the densest point in Gadget data in the
subvolume with left edge [700]*3 and right edge [800]*3 where 700 and
800 are in code units.
This is what I do [1], which was inspired by this page [2],
except I try and take the all_data confined to a region. The output
[3] says however where you will note the center is not inside the
region. It found the densest point in the whole data set, not the
region.
So how would I do this instead for a subvolume of Gadget data? Thanks.
[1] <a href="http://pastebin.com/kyUgEi3a" target="_blank">http://pastebin.com/kyUgEi3a</a>
[2] <a href="http://yt-project.org/doc/cookbook/gadget_notebook.html" target="_blank">http://yt-project.org/doc/cookbook/gadget_notebook.html</a>
[3] <a href="http://pastebin.com/E4aLb5VC" target="_blank">http://pastebin.com/E4aLb5VC</a>
</pre>
</blockquote>
<br>
</div>
</blockquote></div>