<div dir="ltr"><div><div>Hi Nathan,<br><br></div>Thanks for your suggestion. BTW, I have misinterpreted by situation. I'm running in 1/8 part of 3D cartesian coordinates with symmetry. I have looked at the fixed_resolution.py and realized that it can be resolved by modifying FixedResoutionBuffer class to work with mirror symmetry as it handles periodicity. How do you thing about this idea? Do you think it is easy to implement within the yt philosophy?<br>
<br></div>Chang-Goo<br><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Jul 10, 2014 at 9:38 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 Chang-Goo,<div><br></div><div>yt's philosophy for presenting data is to do it exactly as it is represented on disk. Since your data only fills up a fraction of the domain, that is how you see it.</div>
<div><br></div><div>That said, I definitely see why you would want to see the duplicated dataset covering the full domain. To do that, I think you'll need to modify the pixelizer or modify the pixelized representation produced by the pixelizer. You can take a look at how the pixelizers are set up in spherical coordinates in yt/geometry/spherical_coordinates.py, near line 76. Depending on the dimension you are slicing along, yt uses an orthographic pixelizer or a pixelizer optimized for cylindrical coordinates.</div>
<div><br></div><div>You will have to somehow detect that your data only fills half the domain and then mirror it before passing it to the pixelizer. Alternatively you could manipulate the pixelization its self before reutning it from the `pixelize` function..</div>
<div><br></div><div>I think this would be a generally useful addition to yt its self. If you do go down this path, we would very much like to see a pull request with these modifications so that others will be able to generate the pixelations you want.</div>
<div><br></div><div>If you don't want to modify yt, you'll need to do the appropriate data mirroring on the `data` array in your attached script.</div><div><br></div><div>Hope that helps,</div><div><br></div><div>
Nathan</div></div><div class="gmail_extra"><br><br><div class="gmail_quote"><div><div class="h5">On Thu, Jul 10, 2014 at 11:10 AM, Chang-Goo Kim <span dir="ltr"><<a href="mailto:cgkim@astro.princeton.edu" target="_blank">cgkim@astro.princeton.edu</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 class="h5"><div dir="ltr"><div><div>Hi all,<br><br></div>I have recently run a part of sphere by assuming mirror symmetry. So, the simulated domain only spans from 0 to L/2 in every direction, which should be -L/2 ~ L/2 for full domain. I then tried to make frb using slice and to_frb function. What I get is a slice in -L/4 ~L/4 domain with periodicity. How can I get correct slice? Here's a relevant part of script.<br>
<br> for iax in range(3):<br> slc = pf.slice(iax,0.0,center=(0.0,0.0,0.0))<br> slc_frb = slc.to_frb(pf.domain_width[0],512)<br> frb.append(slc_frb)<br><br> for i,field in enumerate(fields):<br>
ax = grid[iax*nrow+i]<br> data=np.array(slc_frb[field])<br> im = ax.imshow(data, origin='lower', norm=LogNorm())<br> im.set_extent(np.array(slc_frb.bounds))<br>
<br></div><div>Thanks,<br><br></div><div>Chang-Goo<br></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" 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">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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>