Covariance calculations are inefficient in LOD imputation

[ericbair-sciome]

Currently, the LOD imputation procedure calculates the full covariance matrix prior to each round of imputation, which wastes both computing time and memory. I think it should be possible to rewrite the function so that it computes only the subset of the covariance matrix that it needs, which should speed things up substantially.