The performance of solvers for integrate-and-fire models with exact spike timing
-
1
RIKEN Brain Science Institute, Japan
-
2
BCCN Freiburg, Germany
Discrete time neuronal network simulation strategies typically constrain spike times to a grid determined by the computational step size. This approach can have the effect of introducing artificial synchrony [1]. However, continuous time (i.e. event-driven) approaches can be computationally demanding, both with respect to calculating future spike times and to event management, particularly for large network sizes. To address this problem, Morrison et al. [2] presented a general method of handling off-grid spiking in combination with exact subthreshold integration in globally time-driven simulations [3,4]. Within each time step the neuron model emulates an event-driven environment to process incoming spikes, whereas the timing of outgoing spikes is based on interpolation. Therefore, the computation step size is a decisive factor for both integration error and simulation time.
An alternative approach for calculating the exact spike times of integrate-and-fire neurons with exponential currents was recently published by Brette [5]. The problem of accurate detection of the first threshold crossing of the membrane potential is converted into finding the largest root of a polynomial. Common numerical means like Descartes' rule and Sturm's theorem are applicable. Although this approach was developed in the context of event-driven simulations, we take advantage of its ability to predict future threshold crossings in the time-driven environment of NEST [3].
We compare the accuracy of the two approaches in single-neuron simulations and the efficiency in a balanced random network of 10,000 neurons [6]. We show that the network simulation time when using the polynomial method depends only weakly on the computational step size, and the single neuron integration error is independent of it. Although the polynomial method attains the maximum precision expected from double numerics for all input rates and computation step sizes, the interpolation method is more efficient for input rates above a critical value.
These results suggest that the cost of processing incoming spikes, rather than the calculation of outgoing spikes, is the dominant factor. Therefore, we extend the model of Morrison et al. [2] by replacing the interpolation of threshold crossings with the computationally more expensive, but numerically more exact, Newton-Raphson technique. The resulting implementation achieves the maximum precision at all computation step sizes and is more efficient for all input and output rates than either the interpolation or the polynomial implementation.
Acknowledgements:
Partially funded by DIP F1.2, BMBF Grant 01GQ0420 to the Bernstein Center for Computational Neuroscience Freiburg, and EU Grant 15879 (FACETS).
References
1. Hansel D et al. (1998) Neural Comput, 10(2):467-483
2. Morrison A et al. (2007) Neural Comput, 19(1):47-79
3. Gewaltig M-O & Diesmann M (2007) Scholarpedia, 2(4):1430
4. Rotter S & Diesmann M (1999) Biol Cybern, 81(5-6):381-402
5. Brette R (2007) Neural Comput, 19(10):2604-2609
6. Brunel N (2000) J Comput Neurosci, 8(3):183-208
Conference:
Neuroinformatics 2008, Stockholm, Sweden, 7 Sep - 9 Sep, 2008.
Presentation Type:
Poster Presentation
Topic:
Computational Neuroscience
Citation:
Diesmann
M,
Hanuschkin
A,
Kunkel
S,
Helias
M and
Morrison
A
(2008). The performance of solvers for integrate-and-fire models with exact spike timing.
Front. Neuroinform.
Conference Abstract:
Neuroinformatics 2008.
doi: 10.3389/conf.neuro.11.2008.01.014
Copyright:
The abstracts in this collection have not been subject to any Frontiers peer review or checks, and are not endorsed by Frontiers.
They are made available through the Frontiers publishing platform as a service to conference organizers and presenters.
The copyright in the individual abstracts is owned by the author of each abstract or his/her employer unless otherwise stated.
Each abstract, as well as the collection of abstracts, are published under a Creative Commons CC-BY 4.0 (attribution) licence (https://creativecommons.org/licenses/by/4.0/) and may thus be reproduced, translated, adapted and be the subject of derivative works provided the authors and Frontiers are attributed.
For Frontiers’ terms and conditions please see https://www.frontiersin.org/legal/terms-and-conditions.
Received:
25 Jul 2008;
Published Online:
25 Jul 2008.
*
Correspondence:
Markus Diesmann, RIKEN Brain Science Institute, Wakoshi, Saitama, Japan, diesmann_frontiers@fz-juelich.de