csiborgtools/notebooks/plot_mass_function.ipynb at 4d7827006a5a772070a514839ba2b8010db6b097

mirror of https://github.com/Richard-Sti/csiborgtools.git synced 2024-12-22 23:18:01 +00:00

Richard Stiskalek 5dd8c668fa

Gaussian smoothing of density fields (#33 )

* Simplify smoothing support and looping over nonzero

* Simplify comments

* add now()

* add cat length

* add smoothed calculation

* add smoothing

* Add sorting

* Edit what is ignored

* Move notebooks

* Add nonsymmetric smoothed overlap

* Update NB

* Add support for reading in the smoothed overlap

* Switch to the true overlap definition

* Reader of the true overlap

* rem occups

* Import moved to a class

* Move definition

* Edit submission script

* Update to account for the new definition

* backup nb

* Switch back to properly initialising arrays

* Fix addition bug

* Update NB

* Fix little bug

* Update nb

2023-03-27 09:22:03 +01:00

104 KiB

Raw Blame History

In [1]:

import numpy as np
import matplotlib.pyplot as plt
from h5py import File
from scipy.stats import spearmanr

import csiborgtools

%matplotlib inline
%load_ext autoreload
%autoreload 2

In [2]:

paths = csiborgtools.read.Paths(**csiborgtools.paths_glamdring)

# d = np.load(paths.field_interpolated("SDSS", "csiborg2_main", 16817, "density", "SPH", 1024))

In [33]:

survey = csiborgtools.SDSS()(apply_selection=False)
# survey = csiborgtools.SDSSxALFALFA()(apply_selection=False)

In [35]:

for kind in ["main", "random"]:
    x, smooth = csiborgtools.summary.read_interpolated_field(survey, f"csiborg2_{kind}", "density", "SPH", 1024, paths)
    np .savez(f"../data/{survey.name}_{kind}_density_SPH_1024.npz", val=x, smooth_scales=smooth)

Reading fields:   0%|          | 0/20 [00:00<?, ?it/s]Reading fields: 100%|██████████| 20/20 [00:11<00:00,  1.80it/s]
Reading fields: 100%|██████████| 20/20 [00:10<00:00,  1.86it/s]

In [37]:

Out[37]:

(20, 641409, 5)

In [24]:

np.load("../data/SDSS_main_density_SPH_1024.npz")["val"]

Out[24]:

array([[[nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        ...,
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan]],

       [[nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        ...,
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan]],

       [[nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        ...,
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan]],

       ...,

       [[nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        ...,
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan]],

       [[nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        ...,
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan]],

       [[nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        ...,
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan],
        [nan, nan, nan, nan, nan]]], dtype=float32)

In [ ]:

104 KiB Raw Blame History

104 KiB

Raw Blame History