We investigate the impact of the role of pairing correlation on the energy per particles of Wigner-Seitz cells in the inner crust of a neutron star. In particular, we compare some common approximations done to treat pairing effects and we estimate the possible error. To reduce the computational cost of the calculations required to determine the chemical composition of the crust, we present a new numerical method based on Gaussian Emulator Process.