Linear minimum-mean-squared error estimation of phase noise, which has a symmetric levy distribution and a possibly large magnitude, from observables at irregular instants.