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

clarify formula

Browse source
This commit is contained in:
markus 2021-01-19 13:10:43 +01:00
parent ee527631f4
commit 3308ffc872
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

@ -22,7 +22,7 @@ def calculate_qr(viral_load: float, emission: float, diameter: float, mask_effic
viral_load = 10 ** viral_load viral_load = 10 ** viral_load
emission = (emission * 3600) if breathing_rate is None else (emission * 1e6) emission = (emission * 3600) if breathing_rate is None else (emission * 1e6)
volume = 4 * np.pi * (diameter / 2) ** 3 / 3 volume = (4 * np.pi * (diameter / 2)**3) / 3
if breathing_rate is None: if breathing_rate is None:
breathing_rate = 1 breathing_rate = 1
return viral_load * emission * volume * (1 - mask_efficiency) * breathing_rate / copies_per_quantum return viral_load * emission * volume * (1 - mask_efficiency) * breathing_rate / copies_per_quantum