<div dir="ltr">I'm getting a strange error when executing the command<div><br></div><div>starmass = ell.quantities['TotalQuantity']('StarMassMsun')[0]<br></div><div><br></div><div>This is on Stampede, and on my laptop running OSX I did not get the same error. Has anyone else seen something similar?</div>
<div><br></div><div>From</div><div>G.S.<br><div><br></div><div><div>In [100]: ell.quantities['TotalQuantity']('StarMassMsun')</div><div>---------------------------------------------------------------------------</div>
<div>ValueError Traceback (most recent call last)</div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/scripts/iyt in <module>()</div><div>----> 1 ell.quantities['TotalQuantity']('StarMassMsun')</div>
<div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/derived_quantities.pyc in __call__(self, *args, **kwargs)</div><div> 80 self._<a href="http://data_source.pf.h.io">data_source.pf.h.io</a>)</div>
<div> 81 if lazy_reader and not self.force_unlazy:</div><div>---> 82 return self._call_func_lazy(args, kwargs)</div><div> 83 else:</div><div> 84 return self._call_func_unlazy(args, kwargs)</div>
<div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/derived_quantities.pyc in _call_func_lazy(self, args, kwargs)</div><div> 87 self.retvals = [ [] for i in range(self.n_ret)]</div><div>
88 for gi,g in enumerate(self._get_grids()):</div><div>---> 89 rv = self.func(GridChildMaskWrapper(g, self._data_source), *args, **kwargs)</div><div> 90 if not iterable(rv): rv = (rv,)</div>
<div> 91 for i in range(self.n_ret): self.retvals[i].append(rv[i])</div><div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/derived_quantities.pyc in _TotalQuantity(data, fields)</div>
<div> 703 totals = []</div><div> 704 for field in fields:</div><div>--> 705 if data[field].size < 1:</div><div> 706 totals.append(0.0)</div><div> 707 continue</div>
<div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/derived_quantities.pyc in __getitem__(self, item)</div><div> 48 def __getitem__(self, item):</div><div> 49 if item not in self.local_cache:</div>
<div>---> 50 data = self.data_source._get_data_from_grid(self.grid, item)</div><div> 51 self.local_cache[item] = data</div><div> 52 return self.local_cache[item]</div><div><br></div>
<div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/data_containers.pyc in save_state(self, grid, field, *args, **kwargs)</div><div> 79 old_keys = grid.field_data.keys()</div><div> 80 grid.field_parameters = self.field_parameters</div>
<div>---> 81 tr = func(self, grid, field, *args, **kwargs)</div><div> 82 grid.field_parameters = old_params</div><div> 83 grid.field_data = YTFieldData( [(k, grid.field_data[k]) for k in old_keys] )</div>
<div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/data_containers.pyc in _get_data_from_grid(self, grid, field)</div><div> 2644 return np.array([f[i,:][pointI] for i in range(3)])</div>
<div> 2645 else:</div><div>-> 2646 tr = grid[field]</div><div> 2647 if tr.size == 1: # dx, dy, dz, cellvolume</div><div> 2648 tr = tr * np.ones(grid.ActiveDimensions, dtype='float64')</div>
<div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in __getitem__(self, key)</div><div> 145 """</div><div> 146 if key not in self.field_data:</div>
<div>--> 147 self.get_data(key)</div><div> 148 return self.field_data[key]</div><div> 149 </div><div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in get_data(self, field, convert)</div>
<div> 188 else: raise</div><div> 189 else:</div><div>--> 190 self._generate_field(field)</div><div> 191 return self.field_data[field]</div><div> 192 </div>
<div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in _generate_field(self, field)</div><div> 133 self[field] = temp_array[sl]</div><div> 134 else:</div>
<div>--> 135 self[field] = self.pf.field_info[field](self)</div><div> 136 else: # Can't find the field, try as it might</div><div> 137 raise exceptions.KeyError(field)</div>
<div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/field_info_container.pyc in __call__(self, data)</div><div> 383 ii = self.check_available(data)</div><div> 384 original_fields = data.keys() # Copy</div>
<div>--> 385 dd = self._function(self, data)</div><div> 386 dd *= self._convert_function(data)</div><div> 387 for field_name in data.keys():</div><div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/universal_fields.pyc in _StarMass(field, data)</div>
<div> 439 </div><div> 440 def _StarMass(field,data):</div><div>--> 441 return data["star_density"] * data["CellVolume"]</div><div> 442 add_field("StarMassMsun", units=r"M_{\odot}",</div>
<div> 443 function=_StarMass,</div><div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in __getitem__(self, key)</div><div> 145 """</div><div> 146 if key not in self.field_data:</div>
<div>--> 147 self.get_data(key)</div><div> 148 return self.field_data[key]</div><div> 149 </div><div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in get_data(self, field, convert)</div>
<div> 188 else: raise</div><div> 189 else:</div><div>--> 190 self._generate_field(field)</div><div> 191 return self.field_data[field]</div><div> 192 </div>
<div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in _generate_field(self, field)</div><div> 133 self[field] = temp_array[sl]</div><div> 134 else:</div>
<div>--> 135 self[field] = self.pf.field_info[field](self)</div><div> 136 else: # Can't find the field, try as it might</div><div> 137 raise exceptions.KeyError(field)</div>
<div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/field_info_container.pyc in __call__(self, data)</div><div> 383 ii = self.check_available(data)</div><div> 384 original_fields = data.keys() # Copy</div>
<div>--> 385 dd = self._function(self, data)</div><div> 386 dd *= self._convert_function(data)</div><div> 387 for field_name in data.keys():</div><div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/frontends/enzo/fields.pyc in _spdensity(field, data)</div>
<div> 360 blank, np.array(data.LeftEdge).astype(np.float64),</div><div> 361 np.array(data.ActiveDimensions).astype(np.int32),</div><div>--> 362 np.float64(data['dx']))</div>
<div> 363 return blank</div><div> 364 </div><div><br></div><div>/home1/01041/tg802953/dev-yt/src/yt-hg/yt/utilities/lib/CICDeposit.so in yt.utilities.lib.CICDeposit.CICDeposit_3 (yt/utilities/lib/CICDeposit.c:1643)()</div>
<div><br></div><div>ValueError: Buffer dtype mismatch, expected 'float64_t' but got 'float'</div></div></div></div>