Bifurcation analysis is an essential tool in nonlinear dynamics. Bifurcation diagrams help to discover subtle features of investigated dynamics such as chaotic and periodic regimes, hidden attractors and fixed points. However, plotting high-resolution bifurcation diagrams can be a computationally challenging task, especially in multiparametric evaluation. It should be noted that the bifurcation analysis is a task with natural parallelism and thus can be efficiently solved using hybrid central-graphics processing architectures. In this paper, we propose an advanced algorithm and special software for plotting bifurcation diagrams using calculations accelerated by graphics processing unit in combination with a highly efficient semi-implicit ordinary differential equation solver. Time series processing is based on the extraction of amplitude and phase features for the density-based spatial clustering of applications with noise to determine oscillation periodicity. We showcase the features of the application of proposed solutions on a set of test chaotic systems. The performance of the analysis algorithms is investigated in comparison with conventional solutions based on central processing unit and several approaches known from the literature. We explicitly show that the proposed algorithm outperforms known solutions in both calculation speed and precision.