def __init__(self, initial_parameters=None, density_contrasts=None, num_filtering_scales=50, path=None, window_function="top hat", number_of_processors=10): """ Instantiates :class:`Inertia` given the Args: """ if initial_parameters is None: initial_parameters = parameters.InitialConditionsParameters( path=path) self.initial_parameters = initial_parameters self.delta_class = density.DensityContrasts( initial_parameters=initial_parameters, num_filtering_scales=num_filtering_scales, window_function=window_function) self.shear_class = shear.Shear( initial_parameters=initial_parameters, num_filtering_scales=num_filtering_scales, number_of_processors=number_of_processors) DensityPeaks.__init__(self, initial_parameters=initial_parameters, density_contrasts=density_contrasts, num_filtering_scales=num_filtering_scales, path=path, window_function=window_function)
def __init__(self, initial_parameters=None, density_contrasts=None, num_filtering_scales=50, path=None, window_function="top hat"): """ Instantiates :class:`DensityPeaks` given the Args: """ if initial_parameters is None: initial_parameters = parameters.InitialConditionsParameters( path=path) self.initial_parameters = initial_parameters self.delta_class = density.DensityContrasts( initial_parameters=initial_parameters, num_filtering_scales=num_filtering_scales, window_function=window_function) if density_contrasts is None: density_contrasts = self.delta_class.density_contrasts self.density_contrasts = density_contrasts assert self.density_contrasts.shape == (num_filtering_scales, initial_parameters.shape**3)
import numpy as np from mlhalos import parameters from mlhalos import window from mlhalos import density import pynbody from regression.adaboost import gbm_04_only as gbm_fun import matplotlib.pyplot as plt from sklearn.externals import joblib # path="/Users/lls/Documents/mlhalos_files/" saving_path_traj = "/share/data2/lls/regression/gradboost/randomly_sampled_training/ic_traj_smoothed_above_1e15Msol/" ic = parameters.InitialConditionsParameters() w = window.WindowParameters(initial_parameters=ic) d = density.DensityContrasts(initial_parameters=ic) m = np.linspace(np.log10(3e10), np.log10(1e15), 50) width = np.append(np.diff(m), np.diff(m)[-1]) m_all_p = ic.initial_conditions["mass"].sum() m1 = np.arange(np.log10(1e15), np.log10(m_all_p*5000), step=width[-1])[1:] M = pynbody.array.SimArray(10**m1) M.units = "Msol" r_smoothing = w.get_smoothing_radius_corresponding_to_filtering_mass(ic, M) den = d.get_smooth_density_for_radii_list(ic, r_smoothing) traj_high_smoothing = den/ic.mean_density np.save(saving_path_traj + "density_contrasts_smoothed_above_1e15Msol.npy", traj_high_smoothing)
min_halo_number=0, max_halo_number=400, min_mass_scale=3e10, max_mass_scale=1e15) sp = shear.ShearProperties(initial_parameters=initial_params, num_filtering_scales=50, shear_scale="all", number_of_processors=60, path=path) ell = sp.density_subtracted_ellipticity prol = sp.density_subtracted_prolateness np.save( "/share/data1/lls/reseed50/features/density_subtracted_ellipticity.npy", ell) np.save( "/share/data1/lls/reseed50/features/density_subtracted_prolateness.npy", prol) Dc = density.DensityContrasts(initial_parameters=initial_params, num_filtering_scales=50, path=path) density_cont = Dc.density_contrasts np.save("/share/data1/lls/reseed50/features/density_contrasts.npy", density_cont) del density_cont del Dc
# path = "/home/lls/stored_files" # # path = "/Users/lls/Documents/mlhalos_files" # # ic = parameters.InitialConditionsParameters(path=path, load_final=True) # z_01 = pynbody.load(path + "/Nina-Simulations/double/snapshot_099") # z_01.physical_units() # # d = density.DensityContrasts(initial_parameters=ic, # snapshot=z_01, # num_filtering_scales=50, window_function="top hat", volume="sphere", path="/home/lls") # # dencon_01 = d.density_contrasts # np.save("/share/data2/lls/features_w_periodicity_fix/z01_density_contrasts.npy", dencon_01) path = "/home/lls/stored_files" # path = "/Users/lls/Documents/mlhalos_files" ic = parameters.InitialConditionsParameters(path=path, load_final=True) z_0 = pynbody.load(path + "/Nina-Simulations/double/snapshot_104") z_0.physical_units() d = density.DensityContrasts(initial_parameters=ic, snapshot=z_0, num_filtering_scales=50, window_function="top hat", volume="sphere", path="/home/lls") dencon_01 = d.density_contrasts np.save("/share/data2/lls/features_w_periodicity_fix/z0_density_contrasts.npy", dencon_01)