can feed a Q matrix to fit.spatautocorel
sample_u tested against samplexuv_with_prior_u

It is faster, as expected, and results seem resonable. Tests are in
spatcontrol.R (but need to define a Q).