Sampling a function periodically replicates its spectrum. As a bilinear function of the signal, the associated Wigner distribution function contains cross terms between the replicas. Often neglected, these cross terms affect numerical simulations of paraxial optical systems. We develop expressions for these cross terms and show their effect on an example calculation.