for n = 1, 2, . . . , N , and Ẑ is given in . The angular frequencies in ω log were taken to be logarithmically spaced with ten frequencies per decade between the absolute value of the minimum and maximum angular frequencies obtained by FFTing the time data t (9). The hyperparameters in were obtained by minimizing (19a) within the set (log 10 λ, log 10 α, . Having assumed Q to be positive definite, was solved using the MOSEK solver in CVXPY . All reported computations were performed on an Intel(R)