The FFT script I posted yesterday does it
http://paste.yt-project.org/show/4528/
I do:
cube = pf.covering_grid(max_level, left_edge=pf.domain_left_edge,
dims=pf.domain_dimensions,
fields=[irho, iu, iv, iw])
then if I don't hack:
pf.periodicity = (True, True, True)
before doing:
rho = cube[irho].d
I get the error. I imagine any dataset would show this error if you hack
one of the periodicity components to False before accessing the cube.
On Fri, Apr 18, 2014 at 12:17 PM, Matthew Turk
Hi Mike,
If I create a covering_grid for a full domain that is not
On Fri, Apr 18, 2014 at 12:02 PM, Michael Zingale
wrote: triply-periodic, the code will abort when I try to access the data in the covering grid with an error generated by selection_routines.pyx that is looking to make sure that the boundary conditions are supported. I understand that this support is not yet in place, but I don't think that it should be an issue for a covering_grid.
I agree, it should not. Can you tell me how you're generating the covering grid? If it is set with the origin at domain_left_edge and dimensions equal to the dimensions of the domain (at whatever refinement level you've specified) this should not error.
The covering grid procedure itself should not ever need boundary
conditions,
if I understand what it intends to do. If you are just sampling the coarse data (or averaging the overlying fine data) into your zone, there is no need to do any interpolation that would require BCs. For the SmoothedCoveringGrid stuff, I understand that BCs should come into play. Am I missing something? or should this check on the BCs be considered a bug in this case?
This sounds like a bug to me!
Mike
-- Michael Zingale Associate Professor
Dept. of Physics & Astronomy * Stony Brook University * Stony Brook, NY 11794-3800 phone: 631-632-8225 e-mail: Michael.Zingale@stonybrook.edu web: http://www.astro.sunysb.edu/mzingale
_______________________________________________ yt-dev mailing list yt-dev@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
_______________________________________________ yt-dev mailing list yt-dev@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
-- Michael Zingale Associate Professor Dept. of Physics & Astronomy * Stony Brook University * Stony Brook, NY 11794-3800 *phone*: 631-632-8225 *e-mail*: Michael.Zingale@stonybrook.edu *web*: http://www.astro.sunysb.edu/mzingale