def raster(self, ftype, component, quick=True): """ ftype is amongst 'gas', 'star', those with has an sml """ field = self.F[ftype] result = zeros(dtype='f4', shape=(self.shape[0], self.shape[1])) ccode.rasterize(field, result, component, xrange=self.cut[0], yrange=self.cut[1], zrange=self.cut[2], quick=quick) return result
def wraster(self, ftype, component, weightcomponent, quick=True): field = self.F[ftype] old = field[component].copy() try: if len(field[component].shape) == 1: field[component][:] *= field[weightcomponent][:] else: field[component][:,0] *= field[weightcomponent][:] field[component][:,1] *= field[weightcomponent][:] field[component][:,2] *= field[weightcomponent][:] result = zeros(dtype='f4', shape = (self.shape[0], self.shape[1])) mass = zeros(dtype='f4', shape = (self.shape[0], self.shape[1])) ccode.rasterize(field, [mass, result], [weightcomponent, component], xrange=self.cut[0], yrange=self.cut[1], zrange=self.cut[2], quick=quick) return result, mass except Exception as e: raise e finally: field[component] = old
def wraster(self, ftype, component, weightcomponent, quick=True): field = self.F[ftype] old = field[component].copy() try: if len(field[component].shape) == 1: field[component][:] *= field[weightcomponent][:] else: field[component][:, 0] *= field[weightcomponent][:] field[component][:, 1] *= field[weightcomponent][:] field[component][:, 2] *= field[weightcomponent][:] result = zeros(dtype='f4', shape=(self.shape[0], self.shape[1])) mass = zeros(dtype='f4', shape=(self.shape[0], self.shape[1])) ccode.rasterize(field, [mass, result], [weightcomponent, component], xrange=self.cut[0], yrange=self.cut[1], zrange=self.cut[2], quick=quick) return result, mass except Exception as e: raise e finally: field[component] = old
def raster(self, ftype, component, quick=True): """ ftype is amongst 'gas', 'star', those with has an sml """ field = self.F[ftype] result = zeros(dtype='f4', shape = (self.shape[0], self.shape[1])) ccode.rasterize(field, result, component, xrange=self.cut[0], yrange=self.cut[1], zrange=self.cut[2], quick=quick) return result