Accurate Reaction-Diffusion Operator Splitting on Tetrahedral Meshes for Parallel Stochastic Molecular Simulations

44 mins 50 secs, 82.01 MB, MP3 44100 Hz, 249.74 kbits/sec

Share this media item:

Embed this media item:

<iframe width="_width_" height="_height_" src="https://sms.cam.ac.uk/media/2274084/embed" frameborder="0" scrolling="no" allowfullscreen></iframe>

Choose size:

About this item

Available Formats

About this item

Description:	De Schutter, E (Okinawa Institute of Science and Technology) Wednesday 22nd June 2016 - 11:45 to 12:30


Created:	2016-07-04 16:16
Collection:	Stochastic Dynamical Systems in Biology: Numerical Methods and Applications
Publisher:	Isaac Newton Institute
Copyright:	De Schutter, E
Language:	eng (English)
Distribution:	World (downloadable)
Explicit content:	No
Aspect Ratio:	16:9
Screencast:	No
Bumper:	UCS Default
Trailer:	UCS Default


Abstract:	Co-authors: Hepburn, Iain (OIST), Chen, Weiliang (OIST) Spatial stochastic molecular simulations in biology are limited by the intense computation required to track molecules in space either by particle tracking or voxel-based methods, meaning that the serial limit has already been reached in sub-cellular models. This calls for parallel simulations that can take advantage of the power of modern supercomputers. GPU parallel implementations have been described for particle tracking methods [1,2] and for voxel-based methods [3], where good parallel performance gain up to 2 order of magnitude have been demonstrated but this depends strongly on model specificity. MPI parallel implementations have gained less attention than GPU implementations to date but offer several advantages including a greater range of platform support from personal computers to advanced supercomputer clusters. An initial MPI implementation for irregular grids has been described and almost ideal speedup demonstrated but only up to 4 cores [4], which indicates the potential for good scalability of such implementations. We describe an operator splitting implementation for irregular grids with a novel method to improve accuracy over Lie-Trotter splitting that is somewhat comparable to tau-reduction but without the performance cost. We systematically investigate parallel performance for a range of models and mesh partitionings using the STEPS simulation platform [5]. Finally we introduce a whole cell parallel simulation of a published reaction-diffusion model [6] within a detailed, complete neuron morphology and demonstrate a speedup of 3 orders of magnitude over serial computations. [1] L Dematte 2012. IEEE/ACM Trans. Comput. Biol. Bioinf. 9: 655-667 [2] DV Gladkov et al. 2011. Proc. 19th High Perf. Comp. Symp. 151-158 [3] E Roberts, JE Stone, Z Luthey-Schulten 2013. J. Comp. Chem. 34: 245–255 [4] A Hellander et al. 2014. J. Comput. Phys. 266: 89-100 [5] I Hepburn et al. 2012. BMC Syst. Biol. 6:36 [6] H Anwar et al. 2013. J. Neurosci. 33: 15848-15867 Related Links http://steps.sourceforge.net/ - Software site

Abstract:

Co-authors: Hepburn, Iain (OIST), Chen, Weiliang (OIST)

Spatial stochastic molecular simulations in biology are limited by the intense computation required to track molecules in space either by particle tracking or voxel-based methods, meaning that the serial limit has already been reached in sub-cellular models. This calls for parallel simulations that can take advantage of the power of modern supercomputers. GPU parallel implementations have been described for particle tracking methods [1,2] and for voxel-based methods [3], where good parallel performance gain up to 2 order of magnitude have been demonstrated but this depends strongly on model specificity. MPI parallel implementations have gained less attention than GPU implementations to date but offer several advantages including a greater range of platform support from personal computers to advanced supercomputer clusters. An initial MPI implementation for irregular grids has been described and almost ideal speedup demonstrated but only up to 4 cores [4], which indicates the potential for good scalability of such implementations.

We describe an operator splitting implementation for irregular grids with a novel method to improve accuracy over Lie-Trotter splitting that is somewhat comparable to tau-reduction but without the performance cost. We systematically investigate parallel performance for a range of models and mesh partitionings using the STEPS simulation platform [5]. Finally we introduce a whole cell parallel simulation of a published reaction-diffusion model [6] within a detailed, complete neuron morphology and demonstrate a speedup of 3 orders of magnitude over serial computations.

[1] L Dematte 2012. IEEE/ACM Trans. Comput. Biol. Bioinf. 9: 655-667 [2] DV Gladkov et al. 2011. Proc. 19th High Perf. Comp. Symp. 151-158 [3] E Roberts, JE Stone, Z Luthey-Schulten 2013. J. Comp. Chem. 34: 245–255 [4] A Hellander et al. 2014. J. Comput. Phys. 266: 89-100 [5] I Hepburn et al. 2012. BMC Syst. Biol. 6:36 [6] H Anwar et al. 2013. J. Neurosci. 33: 15848-15867

Related Links

http://steps.sourceforge.net/ - Software site

Available Formats

Format	Quality	Bitrate	Size
MPEG-4 Video	640x360	1.94 Mbits/sec	651.66 MB	View	Download
iPod Video	480x270	521.51 kbits/sec	171.06 MB	View	Download
MP3 *	44100 Hz	249.74 kbits/sec	82.01 MB	Listen	Download
Auto	(Allows browser to choose a format it supports)

Streaming Media Service Upload

Accurate Reaction-Diffusion Operator Splitting on Tetrahedral Meshes for Parallel Stochastic Molecular Simulations