“…Among the works done on GPU acceleration in option pricing, most have focused on numerical implementation on GPU based on lattice methods such as in Dai, et al [10] , Ganesan, et al [14] , Solomon, et al [15] , Zhang, et al [16] , Suo, et al [17] ; or based on Monte Carlo simulation such as in Michael and Françoise [18] , Podlozhnyuk and Mark [19] , Abbas-Turki and Bernard [20] , Trainor and Crookes [21] , Grauer-Gray, et al [22] and Yu, et al [23] ; or Fourier methods and partial differential approximations such as in Dang [24] , Dang, Christara and Jackson [25] , Surkov [4] and COS approach on GPU such as in Zhang and Oosterlee [26] . These implementations were done for a particular option pricing model and option type such as American, Asian, Basket, European or multi-asset options.…”