I discussed with Perez-Giz and with Foreman-Mackey the creation of quasi-periodic oscillator Gaussian Process models for stars. We want to start by fitting with a damped simple harmonic oscillator kicked by a white-noise source (this has an exact solution as a Gaussian Process, worked out by Goodman and I am sure many others before him). We then want to evolve to non-harmonic oscillators that are better at modeling pulsating stars, but still with tunable incoherence. Applications include: Making the study of quasi-periodic oscillations in compact objects more probabilistic, and more faithful and complete searches for RR Lyrae stars. One problem is that you can't arbitrarily modify your covariance function (kernel function) and obey the rule that it must construct only positive definite variance tensors. I don't really see how to deal with that problem in a simple way, since there is no simple test one can apply to a kernel function that tells you whether or not it is permitted.