Exploring the performance of spatial stochastic simulation algorithms

  • Authors:
  • Matthias Jeschke;Roland Ewald;Adelinde M. Uhrmacher

  • Affiliations:
  • Institute of Computer Science, University of Rostock, Joachim-Jungius-Str. 10, 18059 Rostock, Germany;Institute of Computer Science, University of Rostock, Joachim-Jungius-Str. 10, 18059 Rostock, Germany;Institute of Computer Science, University of Rostock, Joachim-Jungius-Str. 10, 18059 Rostock, Germany

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2011

Quantified Score

Hi-index 31.45

Visualization

Abstract

Since the publication of Gillespie's direct method, diverse methods have been developed to improve the performance of stochastic simulation methods and to enter the spatial realm. In this paper we discuss a spatial @t-leaping variant (S@t) that extends the basic leap method. S@t takes reaction and both outgoing and incoming diffusion events into account when calculating a leap candidate. A performance analysis shall reveal details on the achieved success in balancing speed and accuracy in comparison to other methods. However, performance analysis of spatial stochastic algorithms requires significant effort - it is crucial to choose suitable (benchmark) models and to carefully define model and simulation setups that take problem and simulation design spaces into account.