<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Nov 10, 2015 at 5:16 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"><br><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Tue, Nov 10, 2015 at 4:16 PM, Dave Semeraro <span dir="ltr"><<a href="mailto:semeraro@tacc.utexas.edu" target="_blank">semeraro@tacc.utexas.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div lang="EN-US" link="#0563C1" vlink="#954F72">
<div>
<p class="MsoNormal">Hi YT,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I have an enormous Enzo file that I want to volume render. It takes way too long to render it at full resolution. I have extracted a covering grid and dumped the raw data to a file. Then I read that data into another rendering engine (Boo)
 and rendered it. </p></div></div></blockquote><div><br></div></span><div>You could also do this in yt by reloading the data using load_uniform_grid or by volume rendering a subvolume.</div></div></div></div></blockquote><div><br></div><div>And I should also point out that yt's volume renderer can be both OpenMP and MPI-parallelized and has been used to volume render really huge simulations in the past.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class=""><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" link="#0563C1" vlink="#954F72"><div><p class="MsoNormal">What I want to do is volume render the covering grid. I generate the covering grid like this:<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">max_level=1<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">lowerleft = [0.,0.,0.]<u></u><u></u></p>
<p class="MsoNormal">datasetname = "/work/01197/semeraro/Cooldata"<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">ds = yt.load(datasetname)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">cg = ds.covering_grid(max_level,lowerleft,ds.domain_dimensions*2**max_level)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Then I try to do something like this:<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">min,max = cg.quantities.extrema(‘Density’) <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">and I get a stack trace ending in  <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">  File "/work/01197/semeraro/maverick/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 578, in selector<u></u><u></u></p>
<p class="MsoNormal">    raise YTDataSelectorNotImplemented(self._type_name)<u></u><u></u></p>
<p class="MsoNormal">yt.utilities.exceptions.YTDataSelectorNotImplemented: Data selector 'covering_grid' not implemented.</p></div></div></blockquote><div><br></div></span><div>Can you just do cg['Density'].min(), cg['Density'].max()? covering_grid is a construction data container compared to a selection data container, so I think it's actually an inconsistency in the yt data object API that we're exposing the derived quantities like this.</div><span class=""><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" link="#0563C1" vlink="#954F72"><div><p class="MsoNormal"><u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Obviously I cant do that. So what am I doing wrong here? <u></u>
<u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Thanks, Dave<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Oh and I am making a movie for SC15 and if I get this figured out it will run in the TACC booth.</p></div></div></blockquote><div><br></div></span><div>Cool, I will be there as well. I'll have to check it out :)</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" link="#0563C1" vlink="#954F72"><div><p class="MsoNormal">
<u></u><u></u></p>
</div>
</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>
</blockquote></div><br></div></div>