Is it safe to run a MCMC by interpolating in tabulated data of a model? For background, I have output of a model that involves a set of coupled non-linear differential equations. Calculating models on the fly is slow, so I've tried using a MCMC simulation where the function call is an interpolation in a grid of pre-calculated output, rather than the actual output for a given set of parameters.
My question is, basically, has this kind of application been studied before? I've been Googling for scholarly material for a while and looking in some textbooks for clues but I can't find anything that discusses this. My gut instinct is that the results of the MCMC calculation are as reliable as the interpolation scheme but I'm not sure if there are any formal breakdowns in the subtle details of the chains.
I have a reasonable background in maths and physics but I'm not a master statistician. I apologize if I've overlooked an obvious answer!