Holography is a promising technology for three-dimensional images, but hologram computation is very time-consuming. To accelerate it, we propose a parallel polynomial approximation-based method with a graphical processing unit (GPU) and compare the computation time and quality of reconstruction. It achieved over 15 times faster computations without significant degradations.