%matplotlib inline
%config InlineBackend.figure_format = 'retina'
import matplotlib.pyplot as plt
import numpy as np

single neutrino mass

import sys
sys.path.append('../')
from CEmulator.Emulator import CBaseEmulator
csstemu = CBaseEmulator(verbose=True) # default neutrino_mass_split='single'

Loading the PkcbLin emulator...
Using 513 training samples.
Loading the PknnLin emulator...
Using 512 training samples [remove c0001 (no massive neutrino)].
The neutrino mass is treated as a single massive component.

\(\sigma_{8, \mathrm{cb}}\)

\(\sigma_{8, \mathrm{cb}}\) means the \(\sigma_8\) of the cb components (without the massive neutrino). When mnu = 0.0 eV, the \(\sigma_{8, \mathrm{cb}}\) should be the same as \(\sigma_{8}\)

Ob = 0.05
Oc = 0.25
h = 0.7
ns = 0.96
As = 2.1e-9
w0 = -1
wa = 0
mnu = 0.00
csstemu.set_cosmos(Omegab=Ob, Omegac=Oc, As=As, mnu=mnu,
                   H0=h*100, ns=ns, w=w0, wa=wa)

%time sigma8_emulator = csstemu.get_sigma8_cb(type='Emulator')
print('Emulator: ', sigma8_emulator)
%time sigma8_camb     = csstemu.get_sigma8_cb(type='CAMB')
print('CAMB:     ', sigma8_camb)
%time sigma8_class    = csstemu.get_sigma8_cb(type='CLASS')
print('CLASS:    ', sigma8_class)

CPU times: user 3.84 ms, sys: 727 µs, total: 4.56 ms
Wall time: 4.2 ms
Emulator:  0.8256041401171753
CPU times: user 2.93 s, sys: 27.9 ms, total: 2.96 s
Wall time: 3.62 s
CAMB:      0.8251244044446654
CPU times: user 991 ms, sys: 27.9 ms, total: 1.02 s
Wall time: 1.1 s
CLASS:     0.8250704348558016

print('Emulator/CAMB ',  sigma8_emulator/sigma8_camb)
print('Emulator/CLASS',  sigma8_emulator/sigma8_class)
print('CAMB/CLASS    ',  sigma8_camb/sigma8_class)

Emulator/CAMB  1.0005814101121306
Emulator/CLASS 1.0006468602422616
CAMB/CLASS     1.0000654120988752

\(\sigma_{8}\)

Commonly used \(\sigma_{8}\) include both the cb matter and massive neutrino component.

Ob = 0.05
Oc = 0.25
h = 0.7
ns = 0.96
As = 2.1e-9
w0 = -1
wa = 0
mnu = 0.06
csstemu.set_cosmos(Omegab=Ob, Omegac=Oc, As=As, mnu=mnu,
                   H0=h*100, ns=ns, w=w0, wa=wa)

%time sigma8_emulator = csstemu.get_sigma8(type='Emulator')
print('Emulator: ', sigma8_emulator)
%time sigma8_camb     = csstemu.get_sigma8(type='CAMB')
print('CAMB:     ', sigma8_camb)
%time sigma8_class    = csstemu.get_sigma8(type='CLASS')
print('CLASS:    ', sigma8_class)

CPU times: user 7.14 ms, sys: 2.01 ms, total: 9.14 ms
Wall time: 8.82 ms
Emulator:  0.8133836766724952
CPU times: user 4.3 s, sys: 13 ms, total: 4.31 s
Wall time: 4.32 s
CAMB:      0.8138252232299025
CPU times: user 2.89 s, sys: 14 ms, total: 2.9 s
Wall time: 2.91 s
CLASS:     0.8135117574801464

print('Emulator/CAMB ',  sigma8_emulator/sigma8_camb)
print('Emulator/CLASS',  sigma8_emulator/sigma8_class)
print('CAMB/CLASS    ',  sigma8_camb/sigma8_class)

Emulator/CAMB  0.9994574430174886
Emulator/CLASS 0.999842558135794
CAMB/CLASS     1.0003853241786291

degenerate neutrino mass

import sys
sys.path.append('../')
from CEmulator.Emulator import CBaseEmulator
csstemu = CBaseEmulator(verbose=True, neutrino_mass_split='degenerate')

Loading the PkcbLin emulator...
Using 513 training samples.
Loading the PknnLin emulator...
Using 512 training samples [remove c0001 (no massive neutrino)].
Loading the Tkcb_lin_N3_N1 emulator...
Using 513 training samples.
Loading the Tkmm_lin_N3_N1 emulator...
Using 513 training samples.

\(\sigma_{8, \mathrm{cb}}\)

Ob = 0.05
Oc = 0.25
h = 0.7
ns = 0.96
As = 2.1e-9
w0 = -1
wa = 0
mnu = 0.00
csstemu.set_cosmos(Omegab=Ob, Omegac=Oc, As=As, mnu=mnu,
                   H0=h*100, ns=ns, w=w0, wa=wa)

%time sigma8_emulator = csstemu.get_sigma8(type='Emulator')
print('Emulator: ', sigma8_emulator)
%time sigma8_camb     = csstemu.get_sigma8(type='CAMB')
print('CAMB:     ', sigma8_camb)
%time sigma8_class    = csstemu.get_sigma8(type='CLASS')
print('CLASS:    ', sigma8_class)

CPU times: user 7.43 ms, sys: 1.97 ms, total: 9.4 ms
Wall time: 9.05 ms
Emulator:  0.8251897903850405
CPU times: user 2.77 s, sys: 13 ms, total: 2.79 s
Wall time: 2.79 s
CAMB:      0.8251244044446654
CPU times: user 984 ms, sys: 7.99 ms, total: 992 ms
Wall time: 1.14 s
CLASS:     0.8250704348558016

print('Emulator/CAMB ',  sigma8_emulator/sigma8_camb)
print('Emulator/CLASS',  sigma8_emulator/sigma8_class)
print('CAMB/CLASS    ',  sigma8_camb/sigma8_class)

Emulator/CAMB  1.0000792437358814
Emulator/CLASS 1.0001446610182556
CAMB/CLASS     1.0000654120988752

\(\sigma_{8}\)

Ob = 0.05
Oc = 0.25
h = 0.7
ns = 0.96
As = 2.1e-9
w0 = -1
wa = 0
mnu = 0.06
csstemu.set_cosmos(Omegab=Ob, Omegac=Oc, As=As, mnu=mnu,
                   H0=h*100, ns=ns, w=w0, wa=wa)
print(csstemu.neutrino_mass_split)

degenerate

%time sigma8_emulator = csstemu.get_sigma8(type='Emulator')
print('Emulator: ', sigma8_emulator)
%time sigma8_camb     = csstemu.get_sigma8(type='CAMB')
print('CAMB:     ', sigma8_camb)
%time sigma8_class    = csstemu.get_sigma8(type='CLASS')
print('CLASS:    ', sigma8_class)

CPU times: user 7.72 ms, sys: 2.98 ms, total: 10.7 ms
Wall time: 10.4 ms
Emulator:  0.8148882730802531
CPU times: user 4.87 s, sys: 3.99 ms, total: 4.88 s
Wall time: 4.89 s
CAMB:      0.815649107231685
CPU times: user 7.36 s, sys: 9.06 ms, total: 7.37 s
Wall time: 7.39 s
CLASS:     0.8149908494556329

print('Emulator/CAMB ',  sigma8_emulator/sigma8_camb)
print('Emulator/CLASS',  sigma8_emulator/sigma8_class)
print('CAMB/CLASS    ',  sigma8_camb/sigma8_class)

Emulator/CAMB  0.9990672040897413
Emulator/CLASS 0.9998741380035759
CAMB/CLASS     1.0008076873212648