[yt-users] sub-region projections in parallel

[Matthew Turk]

Hi Stephen,

Right now, parallel projections are only of the entire domain.  The
reason behind this is that the current means of doing domain
decomposition are based on Regions, but Regions cannot themselves be
subdivided.  When I wrote the parallel projections the idea was that
sub-regions of projection would never be too prohibitive for
projecting on a single processor.  I guess this was myopic?

As Eric said in his reply, I think Britton's field_cuts might help you
out here.  The way that came to my mind was that of an extracted
hierarchy, which I'm working on (wherein an entire hierarchy object is
constructed from some subset of an existing hierarchy) but which
probably won't be ready for a couple weeks.

Britton?  Any ideas?

-Matt

On Tue, Feb 3, 2009 at 10:51 AM, Stephen Skory <stephenskory at yahoo.com> wrote:
> Hi all,
>
> I'd like to make some parallel projections for a sub-region. What is the appropriate way to do that? I've tried something like this, below, but it makes a projection for the whole box, not the sub-region. I've tried a few other ways, something like this <http://paste.enzotools.org/show/38/>, but it works in only serial but not in parallel, where it crashes, like this <http://paste.enzotools.org/show/39/> (that's abbreviated output).
>
> i = 1
> pc = PlotCollection(pf.h.region([0.5]*3,[i*.1,0.,0.],[i*.1 + .1,1.,1.]),center=[0.5]*3)
> pc.add_projection("Density", 0)
> pc.save("slice-0.%d" % i)
>
> Thanks!
>
> _______________________________________________________
> sskory at physics.ucsd.edu           o__  Stephen Skory
> http://physics.ucsd.edu/~sskory/ _.>/ _Graduate Student
> ________________________________(_)_\(_)_______________
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>