Mailman 3 Re: yt-dev Digest, Vol 110, Issue 5 - yt-dev

Feb. 20, 2018


      Hi Nathan,
...
thanks! I won't be able to look at this in earnest until later next
week, so if you think you can look at it in the meantime, then great. I
have filed an issue. Sorry I wasn't able to get back to this last week.
Cheers,
  -Mark
On 12 February 2018 at 19:38, <yt-dev-request@python.org> wrote:
...
Send yt-dev mailing list submissions to
        yt-dev@python.org
To subscribe or unsubscribe via the World Wide Web, visit
        https://mail.python.org/mm3/mailman3/lists/yt-dev.python.org/
or, via email, send a message with subject or body 'help' to
        yt-dev-request@python.org
You can reach the person managing the list at
        yt-dev-owner@python.org
When replying, please edit your Subject line so it is more specific
than "Re: Contents of yt-dev digest..."
Today's Topics:
1. OctTree cut_region OffAxis Projections (Mark Richardson)
   2. Re: OctTree cut_region OffAxis Projections (Nathan Goldbaum)
----------------------------------------------------------------------
Date: Mon, 12 Feb 2018 17:49:21 -0500
From: Mark Richardson <mark.richardson.work@gmail.com>
Subject: [yt-dev] OctTree cut_region OffAxis Projections
To: yt-dev@python.org
Message-ID:
        <CA+qw_EzS2jGka-B+jqg0MJBqH5em2SAwi3GNuyb1J_xGtkgP1Q@mail.
gmail.com>
Content-Type: multipart/alternative;
        boundary="001a1148215464247f05650bb037"
--001a1148215464247f05650bb037
Content-Type: text/plain; charset="UTF-8"
Hey, I've previously been able to make off-axis projections of ramses
datasets with geometric data containers. However, I've recently tried to
do
this with a cut_region for dens gas, without success. Here is a script
that
should work with a ramses cosmo dataset:
import yt
import numpy as np
num = 1
Look = [1.,1.,0]
Up = [0.,1.,1.]
c = np.array([0.5,0.5,0.5])
rad = 0.2
width = 2*rad
extra = 1.05
left  = c - rad*extra
right = c + rad*extra
bb = [left,right]
fn = "../output_{0:05d}/info_{0:05d}.txt"
ds = yt.load(fn.format(num),bbox=bb)
region = ds.region(c,left,right)
d20 = region.cut_region(['obj["density"] > 1e-20'])
p1 = yt.OffAxisProjectionPlot(ds,Look,"density",
weight_field='density',center=c,width=width,data_source=regi
on,north_vector=Up)
p1.save("OAProjDens_Region_{0:03d}.png".format(nout))
p1 = yt.OffAxisProjectionPlot(ds,Look,"density",
weight_field='density',center=c,width=width,data_source=d20,
north_vector=Up)
p1.save("OAProjDens_D20_{0:03d}.png".format(nout))
The first OA Projection is generated correctly, but the second one crashes
with the message:
yt : [INFO     ] 2018-02-12 22:20:08,915 Making a fixed resolution buffer
of (('gas', 'density')) 800 by 800
yt : [INFO     ] 2018-02-12 22:22:05,737 Saving plot
OAProjDens_Region_001.png
yt : [INFO     ] 2018-02-12 22:22:06,641 xlim = -0.2 0.2
yt : [INFO     ] 2018-02-12 22:22:06,641 ylim = -0.2 0.2
yt : [INFO     ] 2018-02-12 22:22:06,641 zlim = -0.500000 0.500000
yt : [INFO     ] 2018-02-12 22:23:58,431 Making a fixed resolution buffer
of (('gas', 'density')) 800 by 800
Traceback (most recent call last):
  File "MakeOffAxisProjClouds.py", line 176, in <module>
    p1 = yt.OffAxisProjectionPlot(ds,Look,"density",
weight_field='density',center=c,width=width,data_source=d20,
north_vector=Up)
  File
"/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/plot_window.py",
line 1727, in __init__
    fontsize=fontsize)
  File
"/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/plot_window.py",
line 660, in __init__
    PlotWindow.__init__(self, *args, **kwargs)
  File
"/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/plot_window.py",
line 219, in __init__
    self._setup_plots()
  File
"/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/plot_window.py",
line 748, in _setup_plots
    self._recreate_frb()
  File
"/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/plot_window.py",
line 1733, in _recreate_frb
    super(OffAxisProjectionPlot, self)._recreate_frb()
  File
"/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/plot_window.py",
line 272, in _recreate_frb
    self._frb._get_data_source_fields()
  File
"/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/fixed_resolution.py",
line 155, in _get_data_source_fields
    self[f]
  File
"/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/fixed_resolution.py",
line 543, in __getitem__
    method=dd.method)
  File
"/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/volume_rendering/off_axis_projection.py",
line 189, in off_axis_projection
    for i, (grid, mask) in enumerate(data_source.blocks):
  File
"/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/data_
objects/selection_data_containers.py",
line 836, in blocks
    with obj._field_parameter_state(self.field_parameters):
AttributeError: 'OctreeSubsetBlockSlicePosition' object has no attribute
'_field_parameter_state'
I'm trying to find how to have d20 inherit _field_parameter_state from
region, but thought someone might be able to point me in the right
direction a bit faster.
Thanks!
   -Mark
--
Mark Richardson
MAT Postdoctoral Fellow
Department of Astrophysics
American Museum of Natural History
Mark.Richardson.Work@gmail.com
My Website <https://sites.google.com/site/marklarichardson/>
212 496 3432
--001a1148215464247f05650bb037
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable
...
</div><div><font size=3D"4">Thanks!<br></font><
/div><div>=C2=A0=C2=A0<font=
 size=3D"4"> -Mark</font><br></div><div><br>-- <br><div
class=3D"gmail_sign=
ature"><div dir=3D"ltr"><div><div dir=3D"ltr"><div dir=3D"ltr"><div
<div dir=3D"ltr"><div><font size=3D"4">Hey, I've previously been able
t=
o make off-axis projections of ramses datasets with geometric data
containe=
rs. However, I've recently tried to do this with a cut_region for
dens =
gas, without success. Here is a script that should work with a ramses
cosmo=
 dataset:<br></font><br><span style=3D"font-family:monospace
,monospace">imp=
ort yt<br>import numpy as np<br><br>num =3D 1<br>Look =3D [1.,1.,0]<br>Up
=
=3D [0.,1.,1.]<br><br></span><br><span style=3D"font-family:monospace
,monos=
pace"><span style=3D"font-family:monospace,monospace">c =3D
np.array([0.5,0=
.5,0.5])<br>rad =3D 0.2<br></span></span><span
style=3D"font-family:monospa=
ce,monospace"><span style=3D"font-family:monospace,monospace"><span
style=
=3D"font-family:monospace,monospace">width =3D 2*rad<br></span>extra =3D
1.=
05<br></span></span><span style=3D"font-family:monospace,monospace"><span
s=
tyle=3D"font-family:monospace,monospace">left =C2=A0=3D c -
rad*extra<br>ri=
ght =3D c + rad*extra<br>bb =3D [left,right]<br><br></span>fn =3D
"../=
output_{0:05d}/info_{0:05d}.txt"<br>ds =3D
yt.load(fn.format(num),bbox=
=3Dbb)<br><br>region =3D ds.region(c,left,right)<br><br>d20 =3D
region.cut_=
region(['obj["density"] > 1e-20'])<br><br>p1 =3D
yt.Of=
fAxisProjectionPlot(ds,Look,"density",
weight_field=3D'densit=
y',center=3Dc,width=3Dwidth,data_source=3Dregion,north_
vector=3DUp)<br>=
p1.save("OAProjDens_Region_{0:03d}.png".format(nout))<br><br>p1
=
=3D yt.OffAxisProjectionPlot(ds,Look,"density",
weight_field=3D&#=
39;density',center=3Dc,width=3Dwidth,data_source=3Dd20,
north_vector=3DU=
p)<br>p1.save("OAProjDens_D20_{0:03d}.png".format(
nout))</span><b=
r><br></div><font size=3D"4">The first OA Projection is generated
correctly=
, but the second one crashes with the message: </font><span
style=3D"font-f=
amily:monospace,monospace"><br><br>yt : [INFO =C2=A0 =C2=A0 ] 2018-02-12
22=
:20:08,915 Making a fixed resolution buffer of (('gas',
'densit=
y')) 800 by 800<br>yt : [INFO =C2=A0 =C2=A0 ] 2018-02-12 22:22:05,737
S=
aving plot OAProjDens_Region_001.png<br>yt : [INFO =C2=A0 =C2=A0 ]
2018-02-=
12 22:22:06,641 xlim =3D -0.2 0.2<br>yt : [INFO =C2=A0 =C2=A0 ]
2018-02-12 =
22:22:06,641 ylim =3D -0.2 0.2<br>yt : [INFO =C2=A0 =C2=A0 ] 2018-02-12
22:=
22:06,641 zlim =3D -0.500000 0.500000<br>yt : [INFO =C2=A0 =C2=A0 ]
2018-02=
-12 22:23:58,431 Making a fixed resolution buffer of (('gas',
'=
density')) 800 by 800<br>Traceback (most recent call last):<br>=C2=A0
F=
ile "MakeOffAxisProjClouds.py", line 176, in
<module><br>=
=C2=A0 =C2=A0 p1 =3D yt.OffAxisProjectionPlot(ds,Look,"density",
=
weight_field=3D'density',center=3Dc,width=3Dwidth,da
ta_source=3Dd20=
,north_vector=3DUp)<br>=C2=A0 File "/mnt/blackwhale/home/mri
chard/soft=
/yt/yt_cphyc/yt/visualization/plot_window.py", line 1727, in
__init__<=
br>=C2=A0 =C2=A0 fontsize=3Dfontsize)<br>=C2=A0 File
"/mnt/blackwhale/=
home/mrichard/soft/yt/yt_cphyc/yt/visualization/plot_window.py",
line =
660, in __init__<br>=C2=A0 =C2=A0 PlotWindow.__init__(self, *args,
**kwargs=
)<br>=C2=A0 File "/mnt/blackwhale/home/mri
chard/soft/yt/yt_cphyc/yt/vi=
sualization/plot_window.py", line 219, in __init__<br>=C2=A0 =C2=A0
se=
lf._setup_plots()<br>=C2=A0 File "/mnt/blackwhale/home/mri
chard/soft/y=
t/yt_cphyc/yt/visualization/plot_window.py", line 748, in
_setup_plots=
<br>=C2=A0 =C2=A0 self._recreate_frb()<br>=C2=A0 File
"/mnt/blackwhale=
/home/mrichard/soft/yt/yt_cphyc/yt/visualization/plot_window.py",
line=
 1733, in _recreate_frb<br>=C2=A0 =C2=A0 super(OffAxisProjectionPlot,
self)=
._recreate_frb()<br>=C2=A0 File "/mnt/blackwhale/home/mri
chard/soft/yt=
/yt_cphyc/yt/visualization/plot_window.py", line 272, in
_recreate_frb=
<br>=C2=A0 =C2=A0 self._frb._get_data_source_fields()<br>=C2=A0 File
"=
/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualizat
ion/fixed_resol=
ution.py", line 155, in _get_data_source_fields<br>=C2=A0 =C2=A0
self[=
f]<br>=C2=A0 File "/mnt/blackwhale/home/mri
chard/soft/yt/yt_cphyc/yt/v=
isualization/fixed_resolution.py", line 543, in
__getitem__<br>=C2=A0 =
=C2=A0 method=3Ddd.method)<br>=C2=A0 File "/mnt/blackwhale/home/mri
cha=
rd/soft/yt/yt_cphyc/yt/visualization/volume_rendering/off_
axis_projection.p=
y", line 189, in off_axis_projection<br>=C2=A0 =C2=A0 for i, (grid,
ma=
sk) in enumerate(data_source.blocks):<br>=C2=A0 File
"/mnt/blackwhale/=
home/mrichard/soft/yt/yt_cphyc/yt/data_objects/selection_
data_containers.py=
", line 836, in blocks<br>=C2=A0 =C2=A0 with
obj._field_parameter_stat=
e(self.field_parameters):<br>AttributeError:
'OctreeSubsetBlockSlicePos=
ition' object has no attribute '_field_parameter_state
9;</span><=
br><div><div><br></div><div><font size=3D"4">I'm trying to find how
to =
have d20 inherit _field_parameter_state from region, but thought someone
mi=
ght be able to point me in the right direction a bit faster.
<br><br></font=
dir=3D"=
ltr"><div dir=3D"ltr"><div dir=3D"ltr"><div dir=3D"ltr"><br><span
style=3D"=
font-size:12.8px">Mark Richardson</span><br style=3D"font-size:12.8px"><di
v=
 style=3D"font-size:small"><span style=3D"font-size:12.8px">MAT
Postdoctora=
l Fellow</span><div style=3D"font-size:12.8px">Department of
Astrophysics</=
div><div><span style=3D"font-size:12.8px">American Museum of Natural
Histor=
y</span><div style=3D"font-size:12.8px"><a href=3D"mailto:Mark.Richardson
.W=
ork@gmail.com" style=3D"color:rgb(17,85,204)"
target=3D"_blank">Mark.Richar=
dson.Work@gmail.com</a></div><div><a href=3D"https://sites.google.c
om/site/=
marklarichardson/ <https://sites.google.com/site/=marklarichardson/>"
style=3D"color:rgb(17,85,204)" target=3D"_blank">My Webs=
ite</a><br><div style=3D"font-size:12.8px">212 496 3432
<(212)%20496-3432></div></div></div></=
div></div></div></div></div></div></div></div></div></div>
</div></div></div>
--001a1148215464247f05650bb037--
------------------------------
Date: Mon, 12 Feb 2018 18:37:46 -0600
From: Nathan Goldbaum <nathan12343@gmail.com>
Subject: [yt-dev] Re: OctTree cut_region OffAxis Projections
To: yt-dev@python.org
Message-ID:
        <CAJXewO=eL-F5FtrpY83ddkN=_fh+MC1=Wqadx2bR0=cQrxmvZg@mail.gm
ail.com>
Content-Type: multipart/alternative;
        boundary="94eb2c0eeff24e015605650d35ca"
--94eb2c0eeff24e015605650d35ca
Content-Type: text/plain; charset="UTF-8"
Hi Mark,
This is caused by the fact that OctreeSubsetBlockSlicePosition is not a
subclass of YTSelectionContainer, which implements the appropriate API.
I don't remember offhand why it's designed that way. The fix would either
be to make OctreeSubsetBlockSlicePosition inherit from
YTSelectionContainer
(no idea how hard this would be or what it would break to do that) or
simply implement the _field_parameter_state context manager for the
OctreeSubsetBlockSlicePosition class.
For reference, here's the implementation of the field_parameter_state
context manager for the base data object class:
https://github.com/yt-project/yt/blob/master/yt/data_objects
/data_containers.py#L1084
If you want to take this on to fix the issue you're running into then
that's great. If not then I'd be happy to take a look at this this week.
For the latter option, it would be great if you could file an issue for
this, along with a script that triggers this problem and makes use of one
of the public RAMSES test datasets we have on yt-project.org/data. That
will make it easier for me to start on the fix.
-Nathan
On Mon, Feb 12, 2018 at 4:56 PM Mark Richardson
<mark.richardson.work@gmail.
com> wrote:
...
Hey, I've previously been able to make off-axis projections of ramses
datasets with geometric data containers. However, I've recently tried
to do
this with a cut_region for dens gas, without success. Here is a script
that
should work with a ramses cosmo dataset:
import yt
import numpy as np
num = 1
Look = [1.,1.,0]
Up = [0.,1.,1.]
c = np.array([0.5,0.5,0.5])
rad = 0.2
width = 2*rad
extra = 1.05
left  = c - rad*extra
right = c + rad*extra
bb = [left,right]
fn = "../output_{0:05d}/info_{0:05d}.txt"
ds = yt.load(fn.format(num),bbox=bb)
region = ds.region(c,left,right)
d20 = region.cut_region(['obj["density"] > 1e-20'])
p1 = yt.OffAxisProjectionPlot(ds,Look,"density",
weight_field='density',center=c,width=width,data_source=
region,north_vector=Up)
p1.save("OAProjDens_Region_{0:03d}.png".format(nout))
p1 = yt.OffAxisProjectionPlot(ds,Look,"density",
weight_field='density',center=c,width=width,data_source=d20,
north_vector=Up)
p1.save("OAProjDens_D20_{0:03d}.png".format(nout))
The first OA Projection is generated correctly, but the second one
crashes
with the message:
yt : [INFO     ] 2018-02-12 22:20:08,915 Making a fixed resolution
buffer
of (('gas', 'density')) 800 by 800
yt : [INFO     ] 2018-02-12 22:22:05,737 Saving plot
OAProjDens_Region_001.png
yt : [INFO     ] 2018-02-12 22:22:06,641 xlim = -0.2 0.2
yt : [INFO     ] 2018-02-12 22:22:06,641 ylim = -0.2 0.2
yt : [INFO     ] 2018-02-12 22:22:06,641 zlim = -0.500000 0.500000
yt : [INFO     ] 2018-02-12 22:23:58,431 Making a fixed resolution
buffer
of (('gas', 'density')) 800 by 800
Traceback (most recent call last):
  File "MakeOffAxisProjClouds.py", line 176, in <module>
    p1 = yt.OffAxisProjectionPlot(ds,Look,"density",
weight_field='density',center=c,width=width,data_source=d20,
north_vector=Up)
  File "/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/plot_window.py",
line 1727, in __init__
    fontsize=fontsize)
  File "/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/plot_window.py",
line 660, in __init__
    PlotWindow.__init__(self, *args, **kwargs)
  File "/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/plot_window.py",
line 219, in __init__
    self._setup_plots()
  File "/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/plot_window.py",
line 748, in _setup_plots
    self._recreate_frb()
  File "/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/plot_window.py",
line 1733, in _recreate_frb
    super(OffAxisProjectionPlot, self)._recreate_frb()
  File "/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/visualiza
tion/plot_window.py",
line 272, in _recreate_frb
    self._frb._get_data_source_fields()
  File "/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/
visualization/fixed_resolution.py", line 155, in
_get_data_source_fields
    self[f]
  File "/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/
visualization/fixed_resolution.py", line 543, in __getitem__
    method=dd.method)
  File "/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/
visualization/volume_rendering/off_axis_projection.py", line 189, in
off_axis_projection
    for i, (grid, mask) in enumerate(data_source.blocks):
  File "/mnt/blackwhale/home/mrichard/soft/yt/yt_cphyc/yt/
data_objects/selection_data_containers.py", line 836, in blocks
    with obj._field_parameter_state(self.field_parameters):
AttributeError: 'OctreeSubsetBlockSlicePosition' object has no
attribute
'_field_parameter_state'
I'm trying to find how to have d20 inherit _field_parameter_state from
region, but thought someone might be able to point me in the right
direction a bit faster.
Thanks!
   -Mark
--
Mark Richardson
MAT Postdoctoral Fellow
Department of Astrophysics
American Museum of Natural History
Mark.Richardson.Work@gmail.com
My Website <https://sites.google.com/site/marklarichardson/>
212 496 3432 <(212)%20496-3432>
_______________________________________________
yt-dev mailing list -- yt-dev@python.org
To unsubscribe send an email to yt-dev-leave@python.org
_______________________________________________
yt-dev mailing list -- yt-dev@python.org
To unsubscribe send an email to yt-dev-leave@python.org
--94eb2c0eeff24e015605650d35ca
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable
...
<div>I don't remember offhand why it's designed that way. The
fix =
would either be to make OctreeSubsetBlockSlicePosition inherit from
YTSelec=
tionContainer (no idea how hard this would be or what it would break to
do =
...
</div><div><br><div class=3D"gmail_quote"><div>On Mon, Feb 12, 2018 at
4:5=
6 PM Mark Richardson <<a href=3D"mailto:mark.richardson.work@gmail.com"
=
target=3D"_blank">mark.richardson.work@gmail.<wbr>com</a>>
wrote:<br></d=
iv><blockquote class=3D"gmail_quote" style=3D"margin:0px 0px 0px
0.8ex;bord=
er-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div><font
size=
=3D"4">Hey, I've previously been able to make off-axis projections of
r=
amses datasets with geometric data containers. However, I've recently
t=
ried to do this with a cut_region for dens gas, without success. Here is
a =
...
</span><span style=3D"font-family:monospace,monospace"><span
...
Look,"density", weight_field=3D'density
9;,center=3D<wbr>c=
,width=3Dwidth,data_source=3Dd20,<wbr>north_vector=3DUp)<br>=C2=A0 File
&qu=
ot;/mnt/blackwhale/home/<wbr>mrichard/soft/yt/yt_cphyc/yt/<w
br>visualizatio=
n/plot_window.py", line 1727, in __init__<br>=C2=A0 =C2=A0
fontsize=3D=
fontsize)<br>=C2=A0 File "/mnt/blackwhale/home/<wb
r>mrichard/soft/yt/y=
t_cphyc/yt/<wbr>visualization/plot_window.py", line 660, in
__init__<b=
r>=C2=A0 =C2=A0 PlotWindow.__init__(self, *args, **kwargs)<br>=C2=A0 File
&=
quot;/mnt/blackwhale/home/<wbr>mrichard/soft/yt/yt_cphyc/yt/
<wbr>visualizat=
ion/plot_window.py", line 219, in __init__<br>=C2=A0 =C2=A0
self._setu=
...
mrichard/soft/yt/yt_cphyc/yt/<wbr>visualization/plot_window.py",
...
resolution.py", line 543, in __getitem__<br>=C2=A0 =C2=A0
method=3Ddd=
.method)<br>=C2=A0 File "/mnt/blackwhale/home/<wb
r>mrichard/soft/yt/yt=
_cphyc/yt/<wbr>visualization/volume_<wbr>rendering/off_axis_
...
py", line 189, in off_axis_projection<br>=C2=A0 =C2=A0 for i,
(grid, =
mask) in enumerate(data_source.blocks):<br>=C2=A0 File
"/mnt/blackwhal=
e/home/<wbr>mrichard/soft/yt/yt_cphyc/yt/<wbr>data_objects/s
election_data_<=
wbr>containers.py", line 836, in blocks<br>=C2=A0 =C2=A0 with
obj._fie=
ld_parameter_state(<wbr>self.field_parameters):<br>AttributeError:
'<wb=
r>OctreeSubsetBlockSlicePosition<wbr>' object has no attribute
'_fi=
eld_parameter_state'</span><br><div><div><br></div><div><font
size=3D"4=
">I'm trying to find how to have d20 inherit _field_parameter_state
fro=
m region, but thought someone might be able to point me in the right
...
-- <br><div class=3D"gmail-m_5724714285761475125m_-6614609091351070164gm
ai=
l_signature"><div><div><div><div><div><div><div><div><br><span
...
<div><span style=3D"font-size:12.8px">American Museum of Natural
History</=
span><div style=3D"font-size:12.8px"><a href=3D"mailto:Mark.Richardson
.Work=
@gmail.com" style=3D"color:rgb(17,85,204)" target=3D"_blank">Mark.Richard
so=
n.Work@gmail.com</a></div><div><a href=3D"https://sites.google.c
om/site/mar=
klarichardson/ <https://sites.google.com/site/mar=klarichardson/>"
<div dir=3D"ltr"><div><div dir=3D"auto">Hi Mark,</div><br></div><div>This
i=
s caused by the fact that OctreeSubsetBlockSlicePosition is not a
subclass =
of YTSelectionContainer, which implements the appropriate
API.<br><br></div=
that) or simply implement the _field_parameter_state context manager for
th=
e OctreeSubsetBlockSlicePosition class.<br><br></div><div>For reference,
he=
re's the implementation of the field_parameter_state context manager
fo=
r the base data object class:<br><br><a href=3D"https://github.com/yt-
proje=
ct/yt/blob/master/yt/data_objects/data_containers.py#L1084
<https://github.com/yt-proje=ct/yt/blob/master/yt/data_objects/data_container...>
">https://github.=
com/yt-project/yt/blob/master/yt/data_objects/data_container
s.py#L1084</a><=
br><br></div><div>If you want to take this on to fix the issue you're
r=
unning into then that's great. If not then I'd be happy to take a
l=
ook at this this week. For the latter option, it would be great if you
coul=
d file an issue for this, along with a script that triggers this problem
an=
d makes use of one of the public RAMSES test datasets we have on <a
href=3D=
"http://yt-project.org/data">yt-project.org/data</a>. That will make it
eas=
ier for me to start on the fix.<br><br></div><div>-Nathan
<br></div><div><br=
script that should work with a ramses cosmo dataset:<br></font><br><span
st=
yle=3D"font-family:monospace,monospace">import yt<br>import numpy as
np<br>=
<br>num =3D 1<br>Look =3D [1.,1.,0]<br>Up =3D
[0.,1.,1.]<br><br></span><br>=
<span style=3D"font-family:monospace,monospace"><span
style=3D"font-family:=
monospace,monospace">c =3D np.array([0.5,0.5,0.5])<br>rad =3D
0.2<br></span=
style=3D"font=
-family:monospace,monospace"><span style=3D"font-family:monospace
,monospace=
">width =3D 2*rad<br></span>extra =3D 1.05<br></span></span><span
style=3D"=
font-family:monospace,monospace"><span style=3D"font-family:monospace
,monos=
pace">left =C2=A0=3D c - rad*extra<br>right =3D c + rad*extra<br>bb =3D
[le=
ft,right]<br><br></span>fn =3D "../output_{0:05d}/info_{
0:<wbr>05d}.tx=
t"<br>ds =3D yt.load(fn.format(num),bbox=3D<wbr>bb)<br><br>region
=3D =
ds.region(c,left,right)<br><br>d20 =3D region.cut_region(['obj[&q
uot;<w=
br>density"] > 1e-20'])<br><br>p1 =3D
yt.OffAxisProjectionPlot(=
ds,<wbr>Look,"density", weight_field=3D'density
9;,center=
=3D<wbr>c,width=3Dwidth,data_source=3D<wbr>region,north_vect
or=3DUp)<br>p1.=
save("OAProjDens_Region_{0:<wbr>03d}.png".format(n
out))<br><br>p1=
 =3D yt.OffAxisProjectionPlot(ds,<wbr>Look,"density",
weight_fiel=
d=3D'density',center=3D<wbr>c,width=3Dwidth,data_sou
rce=3Dd20,<wbr>=
north_vector=3DUp)<br>p1.save("OAProjDens_D20_{0:<wbr>0
3d}.png".f=
ormat(nout))</span><br><br></div><font size=3D"4">The first OA
Projection i=
s generated correctly, but the second one crashes with the message:
</font>=
<span style=3D"font-family:monospace,monospace"><br><br>yt : [INFO
=C2=A0 =
=C2=A0 ] 2018-02-12 22:20:08,915 Making a fixed resolution buffer of
(('=
;gas', 'density')) 800 by 800<br>yt : [INFO =C2=A0 =C2=A0 ]
201=
8-02-12 22:22:05,737 Saving plot OAProjDens_Region_001.png<br>yt : [INFO
=
=C2=A0 =C2=A0 ] 2018-02-12 22:22:06,641 xlim =3D -0.2 0.2<br>yt : [INFO
=C2=
=A0 =C2=A0 ] 2018-02-12 22:22:06,641 ylim =3D -0.2 0.2<br>yt : [INFO
=C2=A0=
 =C2=A0 ] 2018-02-12 22:22:06,641 zlim =3D -0.500000 0.500000<br>yt :
[INFO=
 =C2=A0 =C2=A0 ] 2018-02-12 22:23:58,431 Making a fixed resolution buffer
o=
f (('gas', 'density')) 800 by 800<br>Traceback (most
recent=
 call last):<br>=C2=A0 File "MakeOffAxisProjClouds.py", line
176,=
 in <module><br>=C2=A0 =C2=A0 p1 =3D yt.OffAxisProjectionPlot(ds,<w
br=
p_plots()<br>=C2=A0 File "/mnt/blackwhale/home/<wb
r>mrichard/soft/yt/y=
t_cphyc/yt/<wbr>visualization/plot_window.py", line 748, in
_setup_plo=
ts<br>=C2=A0 =C2=A0 self._recreate_frb()<br>=C2=A0 File
"/mnt/blackwha=
le/home/<wbr>mrichard/soft/yt/yt_cphyc/yt/<wbr>visualization
/plot_window.py=
", line 1733, in _recreate_frb<br>=C2=A0 =C2=A0
super(OffAxisProjectio=
nPlot, self)._recreate_frb()<br>=C2=A0 File
"/mnt/blackwhale/home/<wbr=
line=
 272, in _recreate_frb<br>=C2=A0 =C2=A0 self._frb._get_data_source_<wb
r>fie=
lds()<br>=C2=A0 File "/mnt/blackwhale/home/<wb
r>mrichard/soft/yt/yt_cp=
hyc/yt/<wbr>visualization/fixed_<wbr>resolution.py", line 155, in
_get=
_data_source_fields<br>=C2=A0 =C2=A0 self[f]<br>=C2=A0 File
"/mnt/blac=
kwhale/home/<wbr>mrichard/soft/yt/yt_cphyc/yt/<wbr>visualiza
tion/fixed_<wbr=
projection.<wbr=
direct=
ion a bit faster. <br><br></font></div><div><font
size=3D"4">Thanks!<br></f=
ont></div><div>=C2=A0=C2=A0<font size=3D"4">
-Mark</font><br></div><div><br=
style=3D"fon=
t-size:12.8px">Mark Richardson</span><br style=3D"font-size:12.8px"><div
st=
yle=3D"font-size:small"><span style=3D"font-size:12.8px">MAT Postdoctoral
F=
ellow</span><div style=3D"font-size:12.8px">Department of
Astrophysics</div=
style=3D"color:rgb(17,85,204)" target=3D"_blank">My Website=
</a><br><div style=3D"font-size:12.8px"><a href=3D"tel:(212)%20496-3432"
va=
lue=3D"+12124963432 <(212)%20496-3432>" target=3D"_blank">212 496 3432
<(212)%20496-3432></a></div></div></div></=
div></div></div></div></div></div></div></div></div></div>
</div></div></div>
______________________________<wbr>_________________<br>
yt-dev mailing list -- <a href=3D"mailto:yt-dev@python.org"
target=3D"_blan=
k">yt-dev@python.org</a><br>
To unsubscribe send an email to <a href=3D"mailto:yt-dev-leave@python.org"
=
target=3D"_blank">yt-dev-leave@python.org</a><br>
______________________________<wbr>_________________<br>
yt-dev mailing list -- <a href=3D"mailto:yt-dev@python.org"
target=3D"_blan=
k">yt-dev@python.org</a><br>
To unsubscribe send an email to <a href=3D"mailto:yt-dev-leave@python.org"
=
target=3D"_blank">yt-dev-leave@python.org</a><br>
</blockquote></div></div></div>
--94eb2c0eeff24e015605650d35ca--
------------------------------
Subject: Digest Footer
_______________________________________________
yt-dev mailing list -- yt-dev@python.org
To unsubscribe send an email to yt-dev-leave@python.org
------------------------------
End of yt-dev Digest, Vol 110, Issue 5
**************************************
--
Mark Richardson
MAT Postdoctoral Fellow
Department of Astrophysics
American Museum of Natural History
Mark.Richardson.Work@gmail.com
My Website <https://sites.google.com/site/marklarichardson/>
212 496 3432

Re: yt-dev Digest, Vol 110, Issue 5

Mark Richardson

tags

participants (1)