pingu_98/cara
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

fix quanta_exposure

Browse source
This commit is contained in:
markus 2021-02-18 12:25:21 +01:00
parent b0646619ad
commit 8526489e96
1 changed files with 3 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

5
cara/montecarlo.py
View file

@ -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