IRWLS with H as the weighting filter