[yt-users] strange error with StarMassMsun

Geoffrey So gsiisg at gmail.com
Fri Feb 14 14:12:52 PST 2014


I'm getting a strange error when executing the command

starmass = ell.quantities['TotalQuantity']('StarMassMsun')[0]

This is on Stampede, and on my laptop running OSX I did not get the same
error.  Has anyone else seen something similar?

From
G.S.

In [100]: ell.quantities['TotalQuantity']('StarMassMsun')
---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
/home1/01041/tg802953/dev-yt/src/yt-hg/scripts/iyt in <module>()
----> 1 ell.quantities['TotalQuantity']('StarMassMsun')

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/derived_quantities.pyc
in __call__(self, *args, **kwargs)
     80                           self._data_source.pf.h.io)
     81         if lazy_reader and not self.force_unlazy:
---> 82             return self._call_func_lazy(args, kwargs)
     83         else:
     84             return self._call_func_unlazy(args, kwargs)

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/derived_quantities.pyc
in _call_func_lazy(self, args, kwargs)
     87         self.retvals = [ [] for i in range(self.n_ret)]
     88         for gi,g in enumerate(self._get_grids()):
---> 89             rv = self.func(GridChildMaskWrapper(g,
self._data_source), *args, **kwargs)
     90             if not iterable(rv): rv = (rv,)
     91             for i in range(self.n_ret):
self.retvals[i].append(rv[i])

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/derived_quantities.pyc
in _TotalQuantity(data, fields)
    703     totals = []
    704     for field in fields:
--> 705         if data[field].size < 1:
    706             totals.append(0.0)
    707             continue

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/derived_quantities.pyc
in __getitem__(self, item)
     48     def __getitem__(self, item):
     49         if item not in self.local_cache:
---> 50             data = self.data_source._get_data_from_grid(self.grid,
item)
     51             self.local_cache[item] = data
     52         return self.local_cache[item]

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/data_containers.pyc
in save_state(self, grid, field, *args, **kwargs)
     79         old_keys = grid.field_data.keys()
     80         grid.field_parameters = self.field_parameters
---> 81         tr = func(self, grid, field, *args, **kwargs)
     82         grid.field_parameters = old_params
     83         grid.field_data = YTFieldData( [(k, grid.field_data[k]) for
k in old_keys] )

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/data_containers.pyc
in _get_data_from_grid(self, grid, field)
   2644             return np.array([f[i,:][pointI] for i in range(3)])
   2645         else:
-> 2646             tr = grid[field]
   2647             if tr.size == 1: # dx, dy, dz, cellvolume
   2648                 tr = tr * np.ones(grid.ActiveDimensions,
dtype='float64')

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in
__getitem__(self, key)
    145         """
    146         if key not in self.field_data:
--> 147             self.get_data(key)
    148         return self.field_data[key]
    149

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in
get_data(self, field, convert)
    188                     else: raise
    189             else:
--> 190                 self._generate_field(field)
    191         return self.field_data[field]
    192

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in
_generate_field(self, field)
    133                 self[field] = temp_array[sl]
    134             else:
--> 135                 self[field] = self.pf.field_info[field](self)
    136         else: # Can't find the field, try as it might
    137             raise exceptions.KeyError(field)

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/field_info_container.pyc
in __call__(self, data)
    383         ii = self.check_available(data)
    384         original_fields = data.keys() # Copy
--> 385         dd = self._function(self, data)
    386         dd *= self._convert_function(data)
    387         for field_name in data.keys():

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/universal_fields.pyc
in _StarMass(field, data)
    439
    440 def _StarMass(field,data):
--> 441     return data["star_density"] * data["CellVolume"]
    442 add_field("StarMassMsun", units=r"M_{\odot}",
    443           function=_StarMass,

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in
__getitem__(self, key)
    145         """
    146         if key not in self.field_data:
--> 147             self.get_data(key)
    148         return self.field_data[key]
    149

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in
get_data(self, field, convert)
    188                     else: raise
    189             else:
--> 190                 self._generate_field(field)
    191         return self.field_data[field]
    192

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/grid_patch.pyc in
_generate_field(self, field)
    133                 self[field] = temp_array[sl]
    134             else:
--> 135                 self[field] = self.pf.field_info[field](self)
    136         else: # Can't find the field, try as it might
    137             raise exceptions.KeyError(field)

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/data_objects/field_info_container.pyc
in __call__(self, data)
    383         ii = self.check_available(data)
    384         original_fields = data.keys() # Copy
--> 385         dd = self._function(self, data)
    386         dd *= self._convert_function(data)
    387         for field_name in data.keys():

/home1/01041/tg802953/dev-yt/src/yt-hg/yt/frontends/enzo/fields.pyc in
_spdensity(field, data)
    360                            blank,
np.array(data.LeftEdge).astype(np.float64),
    361
 np.array(data.ActiveDimensions).astype(np.int32),
--> 362                            np.float64(data['dx']))
    363     return blank
    364

/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)()

ValueError: Buffer dtype mismatch, expected 'float64_t' but got 'float'
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140214/5117454c/attachment.htm>


More information about the yt-users mailing list