fitting-style TRSHG is done (untested)

This is using the LinearOperator ideas and analytic blockwise inversion
of the sparse matrices