  Motivated by various problems in physics and applied mathematics, we look for
constraints and properties of real Fourier-positive functions, i.e. with
positive Fourier transforms. Properties of the "Dirac comb" distribution and of
its tensor products in higher dimensions lead to Poisson resummation, allowing
for a useful approximation formula of a Fourier transform in terms of a limited
number of terms. A connection with the Bochner theorem on positive definiteness
of Fourier-positive functions is discussed. As a practical application, we find
simple and rapid analytic algorithms for checking Fourier-positivity in 1- and
(radial) 2-dimensions among a large variety of real positive functions. This
may provide a step towards a classification of positive positive-definite
functions.
