Add QPT inversion routines with variances.

--BRJ and MPS