I ran your code with the latest version of PyMC3 and the error you described doesn’t occur, but there is an issue with your covariance matrix D. I think it isn’t positive semidefinite - you should make sure it has all positive eigenvalues. You can see my modifications to your code here: gp-radon-model-broken · GitHub