(color online) (a) Plasma current density distribution of typical EAST upper-single-null configuration and (b) comparison between magnetic flux surface from Green’s function (blue lines) and that from GPU parallel G-S solver (red lines) contributed by plasma current shown in Fig. 5(a), and their difference is so small that blue and red lines coincide with each other.