One of the biggest challenges in tokamak devices is to maintain plasma shape which is affected by constant disturbances. One of the possible solutions for this challenge is usage of model predictive control methodology. Optimization problem which represents used methodology was solved with dual fast gradient method. However, described implementation was too slow. The goal of this thesis is to re-implement the existing implementation using parallel technologies. New implementation must be able to provide solutions with the same level of accuracy as current serial implementation. We attempt to achieve faster execution with parallel program running on general purpose graphics processing units. Chosen technology is OpenCL framework.