Normal Form Linear Solve in Trust Region Methods

Currently we cannot use Normal Form equations in TR. This can be added similar to how Gauss Newton and Levenberg work. Allows us to use Iterative Solvers like GMRES.

In the long run, we want to support LSLQ and such solvers but that requires a minor rewrite of JacVec and VecJac to provide the correct adjoint operators