Improve CMB power spectrum #147
ContinuousIntegration.yml
on: push
Documentation
44m 26s
Matrix: test
Annotations
3 errors and 3 warnings
Documentation:
../../../.julia/packages/Documenter/Bs999/src/utilities/utilities.jl#L44
failed to run `@example` block in src/models.md:9-18
```@example RMΛ
using SymBoltz, Unitful, UnitfulAstro, Plots
M = RMΛ()
pars = Dict(M.r.Ω₀ => 5e-5, M.m.Ω₀ => 0.3, M.g.h => 1.0, M.r.T₀ => 0.0) # TODO: don't pass h and T₀ to avoid infinite loop
ks = [1e-3, 1e-2, 1e-1, 1e-0] / u"Mpc"
sol = solve(M, pars, ks)
p1 = plot(sol, log10(M.g.a), [M.r.ρ, M.m.ρ, M.Λ.ρ, M.G.ρ] ./ M.G.ρ; N = 10000)
p2 = plot(sol, ks, log10(M.g.a), M.g.Φ; N = 10000)
plot(p1, p2, layout = (2, 1), size = (600, 600))
```
exception =
type CosmologyModel has no field b
Stacktrace:
[1] getproperty(M::CosmologyModel, prop::Symbol)
@ SymBoltz ~/work/SymBoltz.jl/SymBoltz.jl/src/solve.jl:53
[2] solve(M::CosmologyModel, pars::Dict{Symbolics.Num, Float64}; aini::Float64, solver::OrdinaryDiffEqRosenbrock.Rodas5P{0, true, Nothing, typeof(OrdinaryDiffEqCore.DEFAULT_PRECS), Val{:forward}, true, nothing, typeof(OrdinaryDiffEqCore.trivial_limiter!), typeof(OrdinaryDiffEqCore.trivial_limiter!)}, reltol::Float64, backwards::Bool, thermo::Bool, debug_initialization::Bool, guesses::Dict{Any, Any}, kwargs::@kwargs{})
@ SymBoltz ~/work/SymBoltz.jl/SymBoltz.jl/src/solve.jl:148
[3] solve
@ ~/work/SymBoltz.jl/SymBoltz.jl/src/solve.jl:98 [inlined]
[4] solve(M::CosmologyModel, pars::Dict{Symbolics.Num, Float64}, ks::Vector{Unitful.Quantity{Float64, 𝐋^-1, Unitful.FreeUnits{(Mpc^-1,), 𝐋^-1, nothing}}}; aini::Float64, solver::OrdinaryDiffEqSDIRK.KenCarp4{0, true, Nothing, OrdinaryDiffEqNonlinearSolve.NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEqCore.DEFAULT_PRECS), Val{:forward}, true, nothing, typeof(OrdinaryDiffEqCore.trivial_limiter!)}, reltol::Float64, backwards::Bool, verbose::Bool, thread::Bool, kwargs::@kwargs{})
@ SymBoltz ~/work/SymBoltz.jl/SymBoltz.jl/src/solve.jl:164
[5] solve(M::CosmologyModel, pars::Dict{Symbolics.Num, Float64}, ks::Vector{Unitful.Quantity{Float64, 𝐋^-1, Unitful.FreeUnits{(Mpc^-1,), 𝐋^-1, nothing}}})
@ SymBoltz ~/work/SymBoltz.jl/SymBoltz.jl/src/solve.jl:159
[6] top-level scope
@ models.md:14
[7] eval
@ ./boot.jl:430 [inlined]
[8] #60
@ ~/.julia/packages/Documenter/Bs999/src/expander_pipeline.jl:803 [inlined]
[9] cd(f::Documenter.var"#60#62"{Module, Expr}, dir::String)
@ Base.Filesystem ./file.jl:112
[10] (::Documenter.var"#59#61"{Documenter.Page, Module, Expr})()
@ Documenter ~/.julia/packages/Documenter/Bs999/src/expander_pipeline.jl:802
[11] (::IOCapture.var"#5#9"{DataType, Documenter.var"#59#61"{Documenter.Page, Module, Expr}, IOContext{Base.PipeEndpoint}, IOContext{Base.PipeEndpoint}, IOContext{Base.PipeEndpoint}, IOContext{Base.PipeEndpoint}})()
@ IOCapture ~/.julia/packages/IOCapture/Y5rEA/src/IOCapture.jl:170
[12] with_logstate(f::IOCapture.var"#5#9"{DataType, Documenter.var"#59#61"{Documenter.Page, Module, Expr}, IOContext{Base.PipeEndpoint}, IOContext{Base.PipeEndpoint}, IOContext{Base.PipeEndpoint}, IOContext{Base.PipeEndpoint}}, logstate::Base.CoreLogging.LogState)
@ Base.CoreLogging ./logging/logging.jl:522
[13] with_logger(f::Function, logger::Base.CoreLogging.ConsoleLogger)
@ Base.CoreLogging ./logging/logging.jl:632
[14] capture(f::Documenter.var"#59#61"{Documenter.Page, Module, Expr}; rethrow::Type, color::Bool, passthrough::Bool, capture_buffer::IOBuffer, io_context::Vector{Any})
@ IOCapture ~/.julia/packages/IOCapture/Y5rEA/src/IOCapture.jl:167
[15] runner(::Type{Documenter.Expanders.ExampleBlocks}, node::MarkdownAST.Node{Nothing}, page::Documenter.Page, doc::Documenter.Document)
@ Documenter ~/.julia/packages/Documenter/Bs999/src/expander_pipeline.jl:801
|
Documentation:
../../../.julia/packages/Documenter/Bs999/src/utilities/utilities.jl#L44
failed to run `@example` block in src/models.md:87-102
```@example BDRMΛ
using SymBoltz, ModelingToolkit, Unitful, UnitfulAstro, Plots
M = SymBoltz.BDRMΛ()
D = Differential(M.t)
pars_fixed = Dict(M.r.Ω₀ => 5e-5, M.m.Ω₀ => 0.3, M.g.h => 1.0, M.r.T₀ => 0.0, M.G.ω => 10.0, D(M.G.ϕ) => 0.0) # unspecified: M.Λ.Ω₀, M.G.ϕ
pars_guess = Dict(M.G.ϕ => 0.95, M.Λ.Ω₀ => 0.7) # initial guesses for shooting method
pars_shoot = shoot(M, pars_fixed, pars_guess, [M.g.ℰ ~ 1, M.G.G ~ 1]; thermo = false, backwards = false) # exact solutions
pars = merge(pars_fixed, pars_shoot) # merge fixed and shooting parameters
k = 1e-0 / u"Mpc"
sol = solve(M, pars, k, backwards = false) # TODO: set background integration direction during model creation
p1 = plot(sol, log10(M.g.a), M.G.G; N = 10000)
p2 = plot(sol, k, log10(M.g.a), M.G.δϕ; N = 50000)
plot(p1, p2, layout = (2, 1))
```
exception =
type CosmologyModel has no field b
Stacktrace:
[1] getproperty(M::CosmologyModel, prop::Symbol)
@ SymBoltz ~/work/SymBoltz.jl/SymBoltz.jl/src/solve.jl:53
[2] solve(M::CosmologyModel, pars::Dict{Symbolics.Num, Float64}; aini::Float64, solver::OrdinaryDiffEqRosenbrock.Rodas5P{0, true, Nothing, typeof(OrdinaryDiffEqCore.DEFAULT_PRECS), Val{:forward}, true, nothing, typeof(OrdinaryDiffEqCore.trivial_limiter!), typeof(OrdinaryDiffEqCore.trivial_limiter!)}, reltol::Float64, backwards::Bool, thermo::Bool, debug_initialization::Bool, guesses::Dict{Any, Any}, kwargs::@kwargs{})
@ SymBoltz ~/work/SymBoltz.jl/SymBoltz.jl/src/solve.jl:148
[3] solve
@ ~/work/SymBoltz.jl/SymBoltz.jl/src/solve.jl:98 [inlined]
[4] solve(M::CosmologyModel, pars::Dict{Symbolics.Num, Float64}, ks::Vector{Unitful.Quantity{Float64, 𝐋^-1, Unitful.FreeUnits{(Mpc^-1,), 𝐋^-1, nothing}}}; aini::Float64, solver::OrdinaryDiffEqSDIRK.KenCarp4{0, true, Nothing, OrdinaryDiffEqNonlinearSolve.NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEqCore.DEFAULT_PRECS), Val{:forward}, true, nothing, typeof(OrdinaryDiffEqCore.trivial_limiter!)}, reltol::Float64, backwards::Bool, verbose::Bool, thread::Bool, kwargs::@kwargs{})
@ SymBoltz ~/work/SymBoltz.jl/SymBoltz.jl/src/solve.jl:164
[5] solve
@ ~/work/SymBoltz.jl/SymBoltz.jl/src/solve.jl:159 [inlined]
[6] #solve#128
@ ~/work/SymBoltz.jl/SymBoltz.jl/src/solve.jl:201 [inlined]
[7] top-level scope
@ models.md:98
[8] eval
@ ./boot.jl:430 [inlined]
[9] #60
@ ~/.julia/packages/Documenter/Bs999/src/expander_pipeline.jl:803 [inlined]
[10] cd(f::Documenter.var"#60#62"{Module, Expr}, dir::String)
@ Base.Filesystem ./file.jl:112
[11] (::Documenter.var"#59#61"{Documenter.Page, Module, Expr})()
@ Documenter ~/.julia/packages/Documenter/Bs999/src/expander_pipeline.jl:802
[12] (::IOCapture.var"#5#9"{DataType, Documenter.var"#59#61"{Documenter.Page, Module, Expr}, IOContext{Base.PipeEndpoint}, IOContext{Base.PipeEndpoint}, IOContext{Base.PipeEndpoint}, IOContext{Base.PipeEndpoint}})()
@ IOCapture ~/.julia/packages/IOCapture/Y5rEA/src/IOCapture.jl:170
[13] with_logstate(f::IOCapture.var"#5#9"{DataType, Documenter.var"#59#61"{Documenter.Page, Module, Expr}, IOContext{Base.PipeEndpoint}, IOContext{Base.PipeEndpoint}, IOContext{Base.PipeEndpoint}, IOContext{Base.PipeEndpoint}}, logstate::Base.CoreLogging.LogState)
@ Base.CoreLogging ./logging/logging.jl:522
[14] with_logger(f::Function, logger::Base.CoreLogging.ConsoleLogger)
@ Base.CoreLogging ./logging/logging.jl:632
[15] capture(f::Documenter.var"#59#61"{Documenter.Page, Module, Expr}; rethrow::Type, color::Bool, passthrough::Bool, capture_buffer::IOBuffer, io_context::Vector{Any})
@ IOCapture ~/.julia/packages/IOCapture/Y5rEA/src/IOCapture.jl:167
[16] runner(::Type{Documenter.Expanders.ExampleBlocks}, node::MarkdownAST.Node{Nothing}, page::Documenter.Page, doc::Documenter.Document)
@ Documenter ~/.julia/packages/Documenter/Bs999/src/expander_pipeline.jl:801
|
Documentation
Process completed with exit code 1.
|
Documentation:
src/solve.jl#L175
Multi-threading was requested, but Julia is running with 1 thread.
|
Documentation:
src/solve.jl#L175
Multi-threading was requested, but Julia is running with 1 thread.
|
Documentation:
src/solve.jl#L175
Multi-threading was requested, but Julia is running with 1 thread.
|