See Improve NUT sampling rate in Gaussian Process when using custom Covariance - #2 by bwengals, which I meant to post here