TaoCreate_XYZ() docstring updates

fixes for code style and single/double/quad builds

initial scaling term changed back to old code

fixing output files

Changes to the initial diagonal scaling vector in SymBrdn propagated across

removed diagonal Broyden matrix and buried it into SymBrdn as a scaling option

options prefixes fixed for derived algorithms

LMVM forward and inverse comparison in Rosenbrock improved

New BQNLS algorithm added based on the BNK kernels

Added diagonal Broyden scaling and preconditioning to BQNK algorithms

Significant improvements to MatMult and MatSolve efficiency for LMVM matrices so that QNK methods become competitive

Quasi-Newton-Krylov methods added that utilize MatMult on LMVM matrices inside a Krylov solver