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

coarsen trapz grid

Browse source
This commit is contained in:
markus 2021-02-19 19:01:29 +01:00
parent 83aaaefab9
commit c7546cec40
1 changed files with 1 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
cara/montecarlo.py
View file

@ -358,7 +358,7 @@ class MCExposureModel:
exposure = np.zeros(self.concentration_model.infected.samples) exposure = np.zeros(self.concentration_model.infected.samples)
for start, stop in self.exposed.presence.boundaries(): for start, stop in self.exposed.presence.boundaries():
times = np.arange(start, stop, 0.01) times = np.arange(start, stop, 0.05)
concentrations = np.asarray([self.concentration_model.concentration(t) for t in times]) concentrations = np.asarray([self.concentration_model.concentration(t) for t in times])
integrals = np.trapz(concentrations, times, axis=0) integrals = np.trapz(concentrations, times, axis=0)
exposure += integrals exposure += integrals