We present a modification of the Rose-Machta algorithm [N. Rose and J. Machta, Phys. Rev. E 100, 063304 (2019)2470-004510.1103/PhysRevE.100.063304] and estimate the density of states for a two-dimensional Blume-Capel model, simulating 10^{5} replicas in parallel for each set of parameters. We perform a finite-size analysis of the specific heat and Binder cumulant, determine the critical temperature along the critical line, and evaluate the critical exponents. The obtained results are in good agreement with those previously obtained using various methods-Markov chain Monte Carlo simulation, Wang-Landau simulation, transfer matrix, and series expansion. The simulation results clearly illustrate the typical behavior of specific heat along the critical lines and through the tricritical point.