sandwich_from_whitened_resid function

GLS standard errors from whitened residuals