Simple PGLS implementation

Probably not actually very useful for least squares fitting now, as
little support for testing different models.  But could be useful
for Bayesian fitting.  No transformations yet.