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. 

Originally posted on yetanothergeographer.tumblr.com.