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)