This adds the new matrix. It has been integrated in the LBFG solver to validate the functionality is complete. Has a few minor cleanups and improvements too.
Closes RAP-183