On efficient computation of multidimensional oscillatory integrals with local Fourier bases