In this paper, we study the linear Schrödinger equation over the d-dimensional torus, with small values of the perturbing potential. We consider numerical approximations of the associated solutions obtained by a symplectic splitting method (to discretize the time variable) in combination with the Fast Fourier Transform algorithm (to discretize the space variable). In this fully discrete setting, we prove that the regularity of the initial datum is preserved over long times, i.e. times that are...