Skip to content
Snippets Groups Projects
Commit 03b8007a authored by Peter Jentsch's avatar Peter Jentsch
Browse files

woops wrong parameter

parent 4d549016
No related branches found
No related tags found
No related merge requests found
......@@ -13,9 +13,9 @@ function get_parameters()#(0.0000,0.00048,0.0005,0.16,-1.30,-1.24,-0.8,0.35,0.35
π_base_y = -1.37,
π_base_m = -1.46,
π_base_o = -0.95,
η = 0.0,
η = 0.00,
κ = 0.0,
ω = 0.0,
ω = 0.0055,
ω_en = 0.0005,
Γ = 1/7,
ξ = 5.0,
......
......@@ -39,7 +39,7 @@ Base.@propagate_inbounds @views function update_alert_durations!(t,modelsol) # B
end
end
@views function update_infection_state!(t,modelsol; record_degrees = false)
Base.@propagate_inbounds @views function update_infection_state!(t,modelsol; record_degrees = false)
@unpack β_y,β_m,β_o,α_y,α_m,α_o,recovery_rate,immunizing,immunization_begin_day = modelsol.params
@unpack u_inf,u_vac,u_next_inf,demographics,inf_network,status_totals, immunization_countdown = modelsol
......@@ -175,11 +175,14 @@ function solve!(modelsol,recording::T) where T
modelsol.u_inf[inf_index] = Infected
modelsol.status_totals[Int(Infected)] += 1
end
end
if t>modelsol.params.immunization_begin_day
update_alert_durations!(t,modelsol)
end
update_vaccination_opinion_state!(t,modelsol,modelsol.status_totals[Int(Infected)])
update_infection_state!(t,modelsol; record_degrees = (T <: HeatmapRecorder))
update_vaccination_opinion_state!(t,modelsol,modelsol.status_totals[Int(Infected)])
#advance agent states based on the new network
modelsol.u_vac .= modelsol.u_next_vac
......
module CovidAlertVaccinationModel
using Intervals: Ending
using Base: Float64
using LightGraphs
using RandomNumbers.Xorshifts
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment