Error variance computation for ordinary krigging?

[Xunius]

Hi,

I'm studying your code for ordinary krigging, and I notice that the variance of krigging error is computed differently from wiki and this reference, and they also differ from each other, so I'm bit confused.

Can you give a look at the relevant code and see if there is any error?