MatShift support added to MATSUBMAT. Mimics MATSHELL functionality. Needed for matrix-free Hessians in TAO Newton-Line-Search methods.

added missing test for recycling the QN approximation in LMVM algorithms

Test updates and bug fixes:
- Added new tests for previously untested LMVM scaling types.
- Slightly changed MATSHELL safeguard on the Hessian to avoid shifts for all virtual inactive submats.
- Typo

Test updates and bug fixes:
- Added new tests for previously untested LMVM scaling types.
- Slightly changed MATSHELL safeguard on the Hessian to avoid shifts for all virtual inactive submats.
- Typo fixed in test flags and associated outputs updated.

show more ...

added safeguard for MATSHELL Hessians

reverting back to using shared MatSolve interface

Custom solve function for LMVM matrices, disconnected from the common MatSolve interface

fixes for code style and single/double/quad builds

initial scaling term changed back to old code

fixing output files

curvature check updated

Changes to the initial diagonal scaling vector in SymBrdn propagated across

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

diagonal brdn scaling now replicates old LMVM code more accurately

options prefixes fixed for derived algorithms

LMVM forward and inverse comparison in Rosenbrock improved

New BQNLS algorithm added based on the BNK kernels

diagonal rescaling added to MATLMVMDIAGBRDN

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

BQNTR algorithm added

BFGS and DFP curvature condition check on updates

LMVM matrices can now define J0 with another LMVM matrix, and DiagBrdn scaling implemented for relevant TAO algorithms

Name refactoring for LMVM matrices

BNK algorithms fixed up for the new LMVM matrix