Covariance bug

I found this bug after getting some nonsense posterior covariances.  The covariance matrix `Kxs` was being overwritten during calculation.

I added a relevant test.  It passes for the basic solver but the last part (positive variance away from training points) fails for HODLR.  I am not sure why.