else:
y_old = y_new
if j<0 or j>= len(A):
print(j, len(A))
res.append(int(j))
rep[j] = rep[j] + 1
return(res, rep)
#___________________________________
# Load settings
conf = cv19.parser()
conf.check_file(argv)
conf.read_config_file()
conf.load_filenames()
conf.load_parameters()
#___________________________________
# parameters
R_0 = 2.2
beta = 0.7
sigma = 0.05
gamma = beta / R_0
population = 100000
# IDEA: # r, rep = random_gen(pdf.cumsum(), population) #}}} #___________________________________ # Load settings config = cv19.parser() config.check_file(argv) config.read_config_file() config.load_filenames() config.load_parameters() p = config.p IDM = cv19.IDModel() # Infectious desease model #___________________________________ # Experiment configuration # number of age bins age_bins = 5