csiborgtools/notebooks/flow/reconstruction_slice.ipynb
Richard Stiskalek 32e36afdc3
Add CLONES support (#145)
* Add submit script

* Add support for gadget2

* Add OMP

* Lower resolution

* Update resolution again

* Add selection on galactic latitude

* Update name

* Update nb

* Update nb

* Add CLONES params

* Add import

* Add CLONES path

* Add CLONES reader

* Quick fix

* Add CLONES

* Update slice nb

* Update nb

* Update nb

* Tiny updates

* Update nb

* Add name

* Add CLONES support

* Add CF4 field

* Add CLONES

* Update script

* Update nb

* Update nb
2024-09-12 17:04:25 +02:00

893 KiB

In [1]:
import sys
import numpy as np
import matplotlib.pyplot as plt
import joblib
import scienceplots
sys.path.append("../")
import csiborgtools

plt.style.use(["science", "notebook"])
%matplotlib widget
%load_ext autoreload
%autoreload 2
In [2]:
d = np.load("/mnt/extraspace/rstiskalek/csiborg/split/clumps_07444_00951.npz")
In [12]:
%timeit d["232"]
772 µs ± 4.01 µs per loop (mean ± std. dev. of 7 runs, 1,000 loops each)
In [ ]:

In [2]:
paths = csiborgtools.read.CSiBORGPaths(**csiborgtools.paths_glamdring)
# cat = csiborgtools.read.ClumpsCatalogue(7444, paths)
In [3]:
paths.knn_path(7444, "auto", "la")
Out[3]:
'/mnt/extraspace/rstiskalek/csiborg/knn/auto/knncdf_07444_la.npz'
In [8]:
np.savez("test.npz", a=np.random.rand(510, 510, 510), b=np.random.rand(510, 510, 510))
In [15]:
d = np.load("test.npz")
In [16]:
d.files
Out[16]:
['a', 'b']
In [7]:
d.files
Out[7]:
['a', 'b']
In [ ]:

In [26]:
np.save("test.npy", np.array([]))
In [27]:
np.load("test.npy")
Out[27]:
array([], dtype=float64)
In [ ]:

In [ ]:

In [21]:
# paths.split_path(123, 7444, 952)
Out[21]:
'/mnt/extraspace/rstiskalek/csiborg/split/ic_00952/out_07444_123.npz'
In [4]:
nsim = 7444
nsnap = max(paths.get_snapshots(7444))
reader = csiborgtools.read.ParticleReader(paths)
In [14]:
# clumpind = reader.read_clumps(nsnap, nsim, cols="index")["index"]
In [ ]:

In [ ]:

In [17]:
parts = np.load("/mnt/extraspace/rstiskalek/csiborg/initmatch/clump_7468_particles.npy", allow_pickle=True)
In [ ]:

In [ ]: