On a model of three-dimensional bursting and its parallel implementation