Want to draw spatially-correlated values wihtout all the fuss?

I think I’ve figured out a simple way to sample from a spatially-correlated normal distribution without needing to invert the large precision matrix directly. 


gist here. 

imported from: yetanothergeographer