def sigma_gen(self, r=None, sigma=None, CDF=None): """ A Wrapper for make_sigma.sigma_gen See make_sigma.sigma_gen for documentation Upon executing, generates sigma, pdf, and cdf_inv and saves to ICobj USAGE: # Generate sigma object from r, sigma arrays sigma = ICobj.maker.sigma_gen(r, sigma) # Use pre-calculated CDF array sigma = ICobj.maker.sigma_gen(r, sigma, CDF) # Generate using a profile defined in sigma_profile.py ICobj.settings.sigma.kind = 'powerlaw' sigma = ICobj.maker.sigma_gen() r and sigma should be 1-D SimArrays. sigma is the surface density evaluated at r """ # Generate sigma if r is None: r, sigma = sigma_profile.make_profile(self._parent) sigma = make_sigma.sigma_gen(r, sigma, CDF) # Copy sigma to the parent (IC) object self._parent.sigma = sigma print 'Sigma stored in <IC instance>.sigma'
def sigma(self, sigma_input): """ Generates a sigma object and stores it in ICobj.rho IF ICobj.settings.sigma.kind = 'file': sigma_input should be a dictionary containing 'sigma' and 'r' Example: sigma_dict = pickle.load(open('sigmafile.p','rb')) IC.settings.sigma.kind = 'file' IC.add.sigma(sigma_dict) """ sigma_out = make_sigma.sigma_gen(self._parent.settings, sigma_input) # Assign to parent ICobj: self._parent.sigma = sigma_out print 'Sigma stored in <IC instance>.sigma'
def sigma_gen(self): """ A Wrapper for make_sigma.sigma_gen See make_sigma.sigma_gen for documentation Upon executing, generates sigma, pdf, and cdf_inv according to settings.sigma.kind and saves to ICobj USAGE: ICobj.maker.sigma_gen() """ # Generate sigma sigma = make_sigma.sigma_gen(self._parent.settings) # Copy sigma to the parent (IC) object self._parent.sigma = sigma print 'Sigma stored in <IC instance>.sigma'
def sigma_gen(self, r, sigma): """ A Wrapper for make_sigma.sigma_gen See make_sigma.sigma_gen for documentation Upon executing, generates sigma, pdf, and cdf_inv and saves to ICobj USAGE: ICobj.maker.sigma_gen(r, sigma) r and sigma should be 1-D SimArrays. sigma is the surface density evaluated at r """ # Generate sigma sigma = make_sigma.sigma_gen(r, sigma) # Copy sigma to the parent (IC) object self._parent.sigma = sigma print 'Sigma stored in <IC instance>.sigma'