Coupling algorithms for calculating sensitivities of Smoluchowski's coagulation equation

From MaRDI portal
Publication:3084397

DOI10.1137/09075679XzbMATH Open1210.82046arXiv0910.3916OpenAlexW2051091969MaRDI QIDQ3084397FDOQ3084397


Authors: Peter L. W. Man, Markus Kraft, James R. Norris, I. Bailleul Edit this on Wikidata


Publication date: 16 March 2011

Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)

Abstract: In this paper, two new stochastic algorithms for calculating parametric derivatives of the solution to the Smoluchowski coagulation equation are presented. It is assumed that the coagulation kernel is dependent on these parameters. The new algorithms (called `Single' and `Double') work by coupling two Marcus-Lushnikov processes in such a way as to reduce the difference between their trajectories, thereby significantly reducing the variance of central difference estimators of the parametric derivatives. In the numerical results, the algorithms are shown have have a O(1/N) order of convergence as expected, where N is the initial number of particles. It was also found that the Single and Double algorithms provide much smaller variances. Furthermore, a method for establishing `efficiency' is considered, which takes into account the variances as well as CPU run times, and the `Double' is significantly more `efficient' compared to the `Independent' algorithm in most cases.


Full work available at URL: https://arxiv.org/abs/0910.3916




Recommendations





Cited In (5)





This page was built for publication: Coupling algorithms for calculating sensitivities of Smoluchowski's coagulation equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q3084397)