This section shows some benchmarking.

using BenchmarkTools
using Epidemics

# an epidemic of 500 days
time_end = 500.0

# benchmark the default model with 3 age groups, intervention, and vaccination
@benchmark epidemic_default(β=[1.3/7], σ=[0.5], γ=[1/7], population = Population(), time_end=time_end, increment=1.0)

BenchmarkTools.Trial: 564 samples with 1 evaluation.
 Range (min … max):  7.983 ms … 49.443 ms  ┊ GC (min … max): 0.00% … 82.89%
 Time  (median):     8.151 ms              ┊ GC (median):    0.00%
 Time  (mean ± σ):   8.837 ms ±  2.972 ms  ┊ GC (mean ± σ):  4.46% ±  9.65%

  █▆      ▃▄                                                  
  ████▇▅▁▄██▇▅▁▄▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▄▁▄▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▁▄▅▅▄▄▁▁▅ ▆
  7.98 ms      Histogram: log(frequency) by time     17.4 ms <

 Memory estimate: 3.78 MiB, allocs estimate: 85939.