Partial implementation of fitting-style TRSHG

This is with LinearOperator ideas and analytic psuedoinverse of the Phi
matrix from "variable projection" methods