This is an implementation of the preconditioned conjugate gradient algorithm over radial basis functions that result in a positive definite interpolation matrix. It also generates an efficient preconditioner and applies it in a cheap way where appropriate.