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

move calculate_qr into emission_rate_when_present

Browse source
This commit is contained in:
markus 2021-02-03 14:56:36 +01:00
parent 6ba9d0c595
commit 388d3c6132
1 changed files with 2 additions and 11 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

13
cara/montecarlo.py
View file

@ -154,18 +154,9 @@ class MCInfectedPopulation(MCPopulation):
breathing_rates = self._generate_breathing_rates()
return qr_func(viral_loads, emission_concentration, mask_efficiency, self.qid, breathing_rates)
viral_loads = 10 ** viral_loads
@staticmethod
def _calculate_qr(viral_load: float, emission_concentration: float, mask_efficiency: float,
copies_per_quantum: float, breathing_rate: float) -> float:
"""
Calculates the quantum generation rate given a set of parameters.
"""
# Unit conversions
viral_load = 10 ** viral_load
return viral_load * emission_concentration * (1 - mask_efficiency) * breathing_rate / copies_per_quantum
return viral_loads * emission_concentration * (1 - mask_efficiency) * breathing_rates / self.qid
def individual_emission_rate(self, time) -> np.ndarray:
"""