Mathematical models are of great value in epidemiology to help understand the dynamics of the various infectious diseases, as well as in the conception of effective control strategies. The classical approach is to use differential equations to describe, in a quantitative manner, the spread of diseases within a particular population. An alternative approach is to represent each individual in the population as a string or vector of characteristic data and simulate the contagion and recovery processes by computational means. This type of model, referred in the literature as MBI (models based on individuals), has the advantage of being flexible as the characteristics of each individual can be quite complex, involving, for instance, age, sex, pre‐existing health conditions, environmental factors, social habits, etc. However, when it comes to simulations involving large populations, MBI may require a large computational effort in terms of memory storage and processing time. In order to cope with the problem of heavy computational effort, this paper proposes a parallel implementation of MBI using a graphics processor unit compatible with CUDA. It was found that, even in the case of a simple susceptible–infected–recovered model, the computational gains in terms of processing time are significant. Copyright © 2015 John Wiley & Sons, Ltd.