[yt-users] Enzo cic_particle_* error

[Casey W. Stark]

So I figured out a way to load particle fields onto AMRCoveringGridBase. I
also had to tack on NumberOfParticles inside the get_data method. cic_*
expects this attribute, and I couldn't figure out a better way to check for
it.

My solution is a serious hack right now and performs a lot of extra reads,
but it would be nice to support this eventually. I will stash it to a
branch for now. Would anyone be interested in coming back to this and
getting it into a good enough state to push back to yt/yt?

Best,
Casey


On Thu, Feb 23, 2012 at 7:47 PM, Casey W. Stark <caseywstark at gmail.com>wrote:

> Okay, that explains it. Thanks Matt.
>
> I will try the shadow region solution tomorrow.
>
> Best,
> Casey
>
>
> On Thu, Feb 23, 2012 at 7:26 PM, Matthew Turk <matthewturk at gmail.com>wrote:
>
>> Hi Casey,
>>
>> On Thu, Feb 23, 2012 at 10:23 PM, Casey W. Stark <caseywstark at gmail.com>
>> wrote:
>> > Hi all.
>> >
>> > I'm trying to make a covering grid of the full domain of an enzo
>> > (cosmological) simulation, getting the dark matter density and velocity
>> > fields. This works fine with the gas:
>> >
>> > In [1]: cg = pf.h.covering_grid(0, [0,0,0], [256,256,256],
>> > fields=["Density", "x-velocity", "y-velocity", "z-velocity"])
>> >
>> > But then I try the DM fields:
>> >
>> > In [2]: cg = pf.h.covering_grid(0, [0,0,0], [256,256,256],
>> > fields=["Dark_Matter_Density", "cic_particle_velocity_x",
>> > "cic_particle_velocity_y", "cic_particle_velocity_z"])
>> >
>> > ... ipython cruft ...
>> >
>> >
>> /project/projectdirs/nyx/yt-installs/yt-hopper/src/yt-hg/yt/data_objects/data_containers.pyc
>> > in __init__(self, level, left_edge, dims, fields, pf, num_ghost_zones,
>> > use_pbar, **kwargs)
>> >    3342         self.domain_width = na.rint((self.pf.domain_right_edge -
>> >    3343
>> > self.pf.domain_left_edge)/self.dds).astype('int64')
>> > -> 3344         self._refresh_data()
>> >    3345
>> >    3346     def _get_list_of_grids(self, buffer = 0.0):
>> >
>> >
>> /project/projectdirs/nyx/yt-installs/yt-hopper/src/yt-hg/yt/data_objects/data_containers.pyc
>> > in _refresh_data(self)
>> >    3359
>> >    3360     def _refresh_data(self):
>> > -> 3361         AMR3DData._refresh_data(self)
>> >    3362         self['dx'] = self.dds[0] *
>> na.ones(self.ActiveDimensions,
>> > dtype='float64')
>> >    3363         self['dy'] = self.dds[1] *
>> na.ones(self.ActiveDimensions,
>> > dtype='float64')
>> >
>> >
>> /project/projectdirs/nyx/yt-installs/yt-hopper/src/yt-hg/yt/data_objects/data_containers.pyc
>> > in _refresh_data(self)
>> >     307         """
>> >     308         self.clear_data()
>> > --> 309         self.get_data()
>> >     310
>> >     311     def keys(self):
>> >
>> >
>> /project/projectdirs/nyx/yt-installs/yt-hopper/src/yt-hg/yt/data_objects/data_containers.pyc
>> > in get_data(self, fields)
>> >    3377                 try:
>> >    3378                     #print "Generating", field
>> >
>> > -> 3379                     self._generate_field(field)
>> >    3380                     continue
>> >    3381                 except NeedsOriginalGrid, ngt_exception:
>> >
>> >
>> /project/projectdirs/nyx/yt-installs/yt-hopper/src/yt-hg/yt/data_objects/data_containers.pyc
>> > in _generate_field(self, field)
>> >    3404             # First we check the validator; this might even
>> raise!
>> >
>> >    3405             self.pf.field_info[field].check_available(self)
>> > -> 3406             self[field] = self.pf.field_info[field](self)
>> >    3407         else: # Can't find the field, try as it might
>> >    3408             raise KeyError(field)
>> >
>> >
>> /project/projectdirs/nyx/yt-installs/yt-hopper/src/yt-hg/yt/data_objects/field_info_container.pyc
>> > in __call__(self, data)
>> >     314         ii = self.check_available(data)
>> >     315         original_fields = data.keys()  # Copy
>> > --> 316         dd = self._function(self, data)
>> >     317         dd *= self._convert_function(data)
>> >     318         for field_name in data.keys():
>> >
>> >
>> /project/projectdirs/nyx/yt-installs/yt-hopper/src/yt-hg/yt/frontends/enzo/fields.pyc
>> > in _cic_particle_field(field, data)
>> >     332     particle_field = field.name[4:]
>> >     333     top = na.zeros(data.ActiveDimensions, dtype='float32')
>> > --> 334     if data.NumberOfParticles == 0: return top
>> >     335     particle_field_data = data[particle_field] *
>> > data['particle_mass']
>> >     336
>> > amr_utils.CICDeposit_3(data["particle_position_x"].astype(na.float64),
>> >
>> > AttributeError: 'AMRCoveringGrid' object has no attribute
>> > 'NumberOfParticles'
>> >
>> > Are you not supposed to use a covering grid with these fields or am I
>> > creating the covering grid incorrectly? I get reasonable looking plots
>> if I
>> > slice/project cic_particle_velocity_x, so I believe it's working
>> otherwise.
>>
>> As it stands, covering grids do not yet have particle support.
>>
>> Adding support should be relatively straightforward; inside get_data,
>> if it's a particle field, set up a shadow region from the left_edge to
>> the right_edge, and then retrieve the field from that.
>>
>> -Matt
>>
>> >
>> > Best,
>> > Casey
>> >
>> > _______________________________________________
>> > yt-users mailing list
>> > yt-users at lists.spacepope.org
>> > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>> >
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20120224/86d93afa/attachment.html>