<div dir="ltr">Hi Wolfram,<div><br></div><div>The error you're seeing happens when the covering grid isn't properly filled with data sampled from the original AMR dataset. Unfortunately, It looks like the covering grid isn't being generated correctly in the current version of yt. I'm no expert in this part of the codebase so I don't think I can help out much in debugging the issue.</div>
<div><br></div><div>Can you go ahead and open an issue about this? The link to open a new issue is here: <a href="https://bitbucket.org/yt_analysis/yt/issues/new">https://bitbucket.org/yt_analysis/yt/issues/new</a></div>
<div>
<br></div><div>There's no need to register a bitbucket account to report an issue, but if you do you will get e-mails when there are updates and when the issue is fixed.</div><div><br></div><div>It would ease debugging and increase the probability that your issue will be fixed if you share a script that runs successfully in the old version of yt and unsuccessfully in the current develoment version, as well as the mercurial changeset hash (obtained by running `yt instinfo` from the command line) for the old yt version. I believe Matt has some Nyx test datasets, but if your script acts on a smallish (< 1 GB) dataset that you can share (via e.g. dropbox, google drive, or a personal website), that would also ease debugging.</div>
<div><br></div><div>Again, sorry to hear about the regression, hopefully we'll be able to figure out the root cause of the issue here.</div><div><br></div><div>Cheers,</div><div><br></div><div>Nathan</div></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Wed, Aug 7, 2013 at 3:17 AM, Wolfram Schmidt <span dir="ltr"><<a href="mailto:schmidt@astro.physik.uni-goettingen.de" target="_blank">schmidt@astro.physik.uni-goettingen.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Nathan,<br>
<br>
thanks a lot for your suggestions! I updated yt to the current development version, which didn't fix the problem. So I tried what you suggested:<br>
<br>
On 08/06/2013 07:52 PM, <a href="mailto:yt-users-request@lists.spacepope.org" target="_blank">yt-users-request@lists.<u></u>spacepope.org</a> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
and then try your script again. If the issue isn't fixed by updating, I<br>
would then try applying the following diff to the nyx frontend:<br>
<br>
diff -r f936432ed45d yt/frontends/nyx/data_<u></u>structures.py<br>
--- a/yt/frontends/nyx/data_<u></u>structures.py Tue Jul 30 10:08:24 2013<br>
-0700<br>
+++ b/yt/frontends/nyx/data_<u></u>structures.py Tue Aug 06 08:10:25 2013<br>
-0700<br>
@@ -611,7 +611,7 @@<br>
# aliases we need<br>
self.parameters["TopGridRank"] =<br>
len(self.parameters["<u></u>TopGridDimensions"]<br>
self.dimensionality = self.parameters["TopGridRank"]<br>
- self.domain_dimensions = self.parameters["<u></u>TopGridDimensions"]<br>
+ self.domain_dimensions =<br>
np.array(self.parameters["<u></u>TopGridDimensions"])<br>
self.refine_by = self.parameters.get("RefineBy"<u></u>, 2) # 2 is silent<br>
defaul<br>
<br>
# Nyx is always cosmological.<br>
<br>
</blockquote>
I replaced the definition of self.domain_dimensions with the np.array thing. Now, the script actually reads in the data, but then another problem is encountered when an FFT is applied to the data:<br>
<br>
Searching grids for values 100% ||||||||||||||||||||||||||||||<u></u>||||||||||||||||||||||||||||||<u></u>||||||||||||||||||||||||||||||<u></u>||||||||||||||||||||||||| Time: 00:03:04<br>
yt : [ERROR ] 2013-08-07 11:56:48,958 Covering problem: 16777216 cells are uncovered<br>
Traceback (most recent call last):<br>
File "/usr/users/wschmid2/scripts/<u></u>spect_sb.py", line 72, in <module><br>
velx_hat = na.fft.fftn( cg['xmom']/cg['density'] )/cg['xmom'].size<br>
File "/home/uni05/wschmid2/yt-i686/<u></u>src/yt-hg/yt/data_objects/<u></u>data_containers.py", line 332, in __getitem__<br>
self.get_data(key)<br>
File "/home/uni05/wschmid2/yt-i686/<u></u>src/yt-hg/yt/data_objects/<u></u>data_containers.py", line 3776, in get_data<br>
raise KeyError(n_bad)<br>
KeyError: 16777216<br>
<br>
The line with the fft actually used to work and it still works just fine if I use an old yt version (which, fortunately, I was able to dig out). Do you understand what is the problem here?<br>
<br>
Cheers,<br>
Wolfram<br>
<br>
______________________________<u></u>_________________<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" target="_blank">http://lists.spacepope.org/<u></u>listinfo.cgi/yt-users-<u></u>spacepope.org</a><br>
</blockquote></div><br></div>