In a previous blog post I described the construction of Confidence Bands for Smoothness. The idea is to express the uncertainty in the amount of smoothing to be applied when fitting functions to data using Gaussian process regression. Now you can do the same thing with INLA. The previous implementation is in STAN. This is quite fast but not as fast as INLA.