From 487c930f8f3f308786a58705352b75e2fa9b3cad Mon Sep 17 00:00:00 2001 From: Nicolas Mounet Date: Wed, 28 Apr 2021 14:20:13 +0200 Subject: [PATCH] Dealing with vectorised concentrations in exposure model --- cara/models.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/cara/models.py b/cara/models.py index bee4a0ff..1a71f1aa 100644 --- a/cara/models.py +++ b/cara/models.py @@ -701,7 +701,7 @@ class ExposureModel: def integrate(fn, start, stop): values = np.linspace(start, stop) - return np.trapz([fn(v) for v in values], values) + return np.trapz([fn(v) for v in values], values, axis=0) for start, stop in self.exposed.presence.boundaries(): exposure += integrate(self.concentration_model.concentration, start, stop)