<div dir="ltr">The new change worked flawlessly, thanks John.<div><br></div><div>From</div><div>G.S.</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Feb 17, 2014 at 12:38 PM, Geoffrey So <span dir="ltr"><<a href="mailto:gsiisg@gmail.com" target="_blank">gsiisg@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Thanks, will give it a try and let you know.<br>
<br>
From<br>
G.S.<br>
<br>
Sent from my iPhone<br>
<div class="HOEnZb"><div class="h5"><br>
> On Feb 17, 2014, at 5:01 AM, John Wise <<a href="mailto:jwise@physics.gatech.edu">jwise@physics.gatech.edu</a>> wrote:<br>
><br>
> Hi Geoffrey,<br>
><br>
> This pull request has been merged, so all you have to do is update your yt installation on Stampede now.<br>
><br>
> Thanks,<br>
> John<br>
><br>
>> On 02/16/2014 08:00 PM, John Wise wrote:<br>
>> Hi Geoffrey,<br>
>><br>
>> This will happen if you're running with 32-bit precision in Enzo, and<br>
>> the particle fields except for position will be in 32-bit. However, the<br>
>> CICDeposit3 routine expects 64-bit floats as input. I've created a pull<br>
>> request that fixes this bug.<br>
>><br>
>> <a href="https://bitbucket.org/yt_analysis/yt/pull-request/706/input-enzo-deposited-fields-should-be/diff" target="_blank">https://bitbucket.org/yt_analysis/yt/pull-request/706/input-enzo-deposited-fields-should-be/diff</a><br>
>><br>
>><br>
>> You can go ahead and test it by pulling my changeset and recompiling.<br>
>><br>
>> hg pull -r f9127d6 <a href="https://bitbucket.org/jwise77/yt" target="_blank">https://bitbucket.org/jwise77/yt</a><br>
>> hg up<br>
>><br>
>> Please see if this works for you and comment on the PR page.<br>
>><br>
>> Cheers,<br>
>> John<br>
>><br>
>>> On 02/14/2014 05:12 PM, Geoffrey So wrote:<br>
>>> I'm getting a strange error when executing the command<br>
>>><br>
>>> starmass = ell.quantities['TotalQuantity']('StarMassMsun')[0]<br>
>>><br>
>>> This is on Stampede, and on my laptop running OSX I did not get the same<br>
>>> error. Has anyone else seen something similar?<br>
>>><br>
>>> From<br>
>>> G.S.<br>
>>><br>
>>> In [100]: ell.quantities['TotalQuantity']('StarMassMsun')<br>
>>> ---------------------------------------------------------------------------<br>
>>><br>
>>> ValueError Traceback (most recent call<br>
>>> last)<br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/scripts/iyt in <module>()<br>
>>> ----> 1 ell.quantities['TotalQuantity']('StarMassMsun')<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/derived_quantities.pyc<br>
>>><br>
>>> in __call__(self, *args, **kwargs)<br>
>>> 80 self._<a href="http://data_source.pf.h.io" target="_blank">data_source.pf.h.io</a><br>
>>> <<a href="http://data_source.pf.h.io" target="_blank">http://data_source.pf.h.io</a>>)<br>
>>> 81 if lazy_reader and not self.force_unlazy:<br>
>>> ---> 82 return self._call_func_lazy(args, kwargs)<br>
>>> 83 else:<br>
>>> 84 return self._call_func_unlazy(args, kwargs)<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/derived_quantities.pyc<br>
>>><br>
>>> in _call_func_lazy(self, args, kwargs)<br>
>>> 87 self.retvals = [ [] for i in range(self.n_ret)]<br>
>>> 88 for gi,g in enumerate(self._get_grids()):<br>
>>> ---> 89 rv = self.func(GridChildMaskWrapper(g,<br>
>>> self._data_source), *args, **kwargs)<br>
>>> 90 if not iterable(rv): rv = (rv,)<br>
>>> 91 for i in range(self.n_ret):<br>
>>> self.retvals[i].append(rv[i])<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/derived_quantities.pyc<br>
>>><br>
>>> in _TotalQuantity(data, fields)<br>
>>> 703 totals = []<br>
>>> 704 for field in fields:<br>
>>> --> 705 if data[field].size < 1:<br>
>>> 706 totals.append(0.0)<br>
>>> 707 continue<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/derived_quantities.pyc<br>
>>><br>
>>> in __getitem__(self, item)<br>
>>> 48 def __getitem__(self, item):<br>
>>> 49 if item not in self.local_cache:<br>
>>> ---> 50 data =<br>
>>> self.data_source._get_data_from_grid(self.grid, item)<br>
>>> 51 self.local_cache[item] = data<br>
>>> 52 return self.local_cache[item]<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/data_containers.pyc<br>
>>><br>
>>> in save_state(self, grid, field, *args, **kwargs)<br>
>>> 79 old_keys = grid.field_data.keys()<br>
>>> 80 grid.field_parameters = self.field_parameters<br>
>>> ---> 81 tr = func(self, grid, field, *args, **kwargs)<br>
>>> 82 grid.field_parameters = old_params<br>
>>> 83 grid.field_data = YTFieldData( [(k, grid.field_data[k])<br>
>>> for k in old_keys] )<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/data_containers.pyc<br>
>>><br>
>>> in _get_data_from_grid(self, grid, field)<br>
>>> 2644 return np.array([f[i,:][pointI] for i in range(3)])<br>
>>> 2645 else:<br>
>>> -> 2646 tr = grid[field]<br>
>>> 2647 if tr.size == 1: # dx, dy, dz, cellvolume<br>
>>> 2648 tr = tr * np.ones(grid.ActiveDimensions,<br>
>>> dtype='float64')<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in<br>
>>> __getitem__(self, key)<br>
>>> 145 """<br>
>>> 146 if key not in self.field_data:<br>
>>> --> 147 self.get_data(key)<br>
>>> 148 return self.field_data[key]<br>
>>> 149<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in<br>
>>> get_data(self, field, convert)<br>
>>> 188 else: raise<br>
>>> 189 else:<br>
>>> --> 190 self._generate_field(field)<br>
>>> 191 return self.field_data[field]<br>
>>> 192<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in<br>
>>> _generate_field(self, field)<br>
>>> 133 self[field] = temp_array[sl]<br>
>>> 134 else:<br>
>>> --> 135 self[field] = self.pf.field_info[field](self)<br>
>>> 136 else: # Can't find the field, try as it might<br>
>>> 137 raise exceptions.KeyError(field)<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/field_info_container.pyc<br>
>>><br>
>>> in __call__(self, data)<br>
>>> 383 ii = self.check_available(data)<br>
>>> 384 original_fields = data.keys() # Copy<br>
>>> --> 385 dd = self._function(self, data)<br>
>>> 386 dd *= self._convert_function(data)<br>
>>> 387 for field_name in data.keys():<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/universal_fields.pyc<br>
>>><br>
>>> in _StarMass(field, data)<br>
>>> 439<br>
>>> 440 def _StarMass(field,data):<br>
>>> --> 441 return data["star_density"] * data["CellVolume"]<br>
>>> 442 add_field("StarMassMsun", units=r"M_{\odot}",<br>
>>> 443 function=_StarMass,<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in<br>
>>> __getitem__(self, key)<br>
>>> 145 """<br>
>>> 146 if key not in self.field_data:<br>
>>> --> 147 self.get_data(key)<br>
>>> 148 return self.field_data[key]<br>
>>> 149<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in<br>
>>> get_data(self, field, convert)<br>
>>> 188 else: raise<br>
>>> 189 else:<br>
>>> --> 190 self._generate_field(field)<br>
>>> 191 return self.field_data[field]<br>
>>> 192<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in<br>
>>> _generate_field(self, field)<br>
>>> 133 self[field] = temp_array[sl]<br>
>>> 134 else:<br>
>>> --> 135 self[field] = self.pf.field_info[field](self)<br>
>>> 136 else: # Can't find the field, try as it might<br>
>>> 137 raise exceptions.KeyError(field)<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/field_info_container.pyc<br>
>>><br>
>>> in __call__(self, data)<br>
>>> 383 ii = self.check_available(data)<br>
>>> 384 original_fields = data.keys() # Copy<br>
>>> --> 385 dd = self._function(self, data)<br>
>>> 386 dd *= self._convert_function(data)<br>
>>> 387 for field_name in data.keys():<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/frontends/enzo/fields.pyc in<br>
>>> _spdensity(field, data)<br>
>>> 360 blank,<br>
>>> np.array(data.LeftEdge).astype(np.float64),<br>
>>> 361<br>
>>> np.array(data.ActiveDimensions).astype(np.int32),<br>
>>> --> 362 np.float64(data['dx']))<br>
>>> 363 return blank<br>
>>> 364<br>
>>><br>
>>> /home1/01041/tg802953/dev-yt/src/yt-hg/yt/utilities/lib/CICDeposit.so in<br>
>>> yt.utilities.lib.CICDeposit.CICDeposit_3<br>
>>> (yt/utilities/lib/CICDeposit.c:1643)()<br>
>>><br>
>>> ValueError: Buffer dtype mismatch, expected 'float64_t' but got 'float'<br>
>>><br>
>>><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>
> --<br>
> John Wise<br>
> Assistant Professor of Physics<br>
> Center for Relativistic Astrophysics, Georgia Tech<br>
> <a href="http://cosmo.gatech.edu" target="_blank">http://cosmo.gatech.edu</a><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>
</div></div></blockquote></div><br></div>