Commit b16102aa authored by Peter Jentsch's avatar Peter Jentsch
Browse files

time dep pi_base

parent a7519605
......@@ -3,9 +3,9 @@ function get_parameters()#(0.0000,0.00048,0.0005,0.16,-1.30,-1.24,-0.8,0.35,0.35
sim_length = 400,
num_households = 5000,
I_0_fraction = 0.005,
β_y = 0.00046,
β_m = 0.00044,
β_o = 0.6,
β_y = 0.00047,
β_m = 0.00041,
β_o = 0.75,
α_y = 0.4,
α_m = 0.4,
α_o = 0.4,
......@@ -15,7 +15,7 @@ function get_parameters()#(0.0000,0.00048,0.0005,0.16,-1.30,-1.24,-0.8,0.35,0.35
π_base_o = -0.95,
η = 0.0,
κ = 0.0,
ω = 0.00,
ω = 0.0005,
ω_en = 0.00,
γ = 0.0,
ξ = 5.0,
......
......@@ -98,11 +98,13 @@ end
Base.@propagate_inbounds @views function update_vaccination_opinion_state!(t,modelsol,total_infections)
@unpack π_base_y,π_base_m,π_base_o, η,γ, κ, ω, ω_en,γ,ξ = modelsol.params
@unpack infection_introduction_day, π_base_y,π_base_m,π_base_o, η,γ, κ, ω, ω_en,γ,ξ = modelsol.params
@unpack demographics,time_of_last_alert, nodes, soc_network,u_vac,u_next_vac,app_user,app_user_list = modelsol
app_user_pointer = 0
for i in 1:nodes
π_base = @SVector [π_base_y,π_base_m,π_base_o]
π_base = t<infection_introduction_day ?
(π_base_y,π_base_m,π_base_o) :
(π_base_y*1.25,π_base_m*1.25,π_base_o)
random_soc_network = sample(Random.default_rng(Threads.threadid()), soc_network.graph_list[t])
if !isempty(neighbors(random_soc_network,i))
......
......@@ -82,5 +82,5 @@ end
@testset "degree distribution generation" for _ in samples
end
\ No newline at end of file
No preview for this file type
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment