From 8526489e96870b3ccf91b916fedbf2ed181e0a46 Mon Sep 17 00:00:00 2001 From: markus Date: Thu, 18 Feb 2021 12:25:21 +0100 Subject: [PATCH] fix quanta_exposure --- cara/montecarlo.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/cara/montecarlo.py b/cara/montecarlo.py index d45acc4a..3e0c208e 100644 --- a/cara/montecarlo.py +++ b/cara/montecarlo.py @@ -344,8 +344,9 @@ class MCExposureModel: exposure = np.zeros(self.concentration_model.infected.samples) for start, stop in self.exposed.presence.boundaries(): - concentrations = np.asarray([self.concentration_model.concentration(t) for t in np.linspace(start, stop)]) - integrals = np.trapz(concentrations, axis=0) + times = np.arange(start, stop, 0.001) + concentrations = np.asarray([self.concentration_model.concentration(t) for t in times]) + integrals = np.trapz(concentrations, times, axis=0) exposure += integrals return exposure * self.repeats