Add Planck18 Neff parameter

hersle · Oct 15, 2024 · e065509 · e065509
1 parent e060d30
commit e065509
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/docs/src/comparison.md b/docs/src/comparison.md
@@ -24,7 +24,6 @@ using Printf
 lmax = 6
 M = SymBoltz.ΛCDM(; lmax, h = nothing, Λanalytical = true)
 pars = SymBoltz.parameters_Planck18(M)
-pars[M.ν.Neff] = 3.044 # TODO
 
 function run_class(in::Dict{String, Any}, exec, inpath, outpath)
     merge!(in, Dict(
@@ -113,7 +112,7 @@ end
 
 k = 1e0 / u"Mpc" # 1/Mpc # disagreement on smaller scales
 sol1 = solve_class(pars, k)
-sol2 = solve(M, pars, k; solver = Rodas5P()) # TODO: KenCarp4 and Kvaerno5 "emulate" radiation streaming, while Rodas5P continues in an exact way
+sol2 = solve(M, pars, k; solver = SymBoltz.Rodas5P()) # TODO: KenCarp4 and Kvaerno5 "emulate" radiation streaming, while Rodas5P continues in an exact way
 
 # map results from both codes to common convention
 h = pars[M.g.h]

diff --git a/src/parameters.jl b/src/parameters.jl
@@ -5,8 +5,9 @@ function parameters_Planck18(M::CosmologyModel)
         M.g.h => h,
         M.γ.T0 => 2.7255,
         M.c.Ω0 => 0.1200 / h^2,
-        M.b.Ω0 => 0.02237 / h^2,
-        M.b.rec.Yp => 0.2454
+        M.b.Ω0 => 0.0224 / h^2,
+        M.b.rec.Yp => 0.2454,
+        M.ν.Neff => 2.99 # TODO: vs 3.046 or 3.044?
 
         # TODO: how to handle
         # 1) backwards with default for Λ.Ω0? --> skip E = 1