The accuracy and computational cost of the Gaussian puff model largely depends on two parameters: the frequency at which the simulations are evaluated and the frequency at which puffs are created, denoted by ∆t and ∆p, respectively. The former affects the accuracy of each single puff simulation, with smaller ∆t increasing accuracy x ← computePuffPath(Ω, u)