A new, computationally- and statistically-efficient algorithm, the Fast
χ2 algorithm, can find a periodic signal with harmonic content in
irregularly-sampled data with non-uniform errors. The algorithm calculates the
minimized χ2 as a function of frequency at the desired number of
harmonics, using Fast Fourier Transforms to provide O(NlogN) performance.
The code for a reference implementation is provided.Comment: Source code for the reference implementation is available at
http://public.lanl.gov/palmer/fastchi.html . Accepted by ApJ. 24 pages, 4
figure