Steihaug–Toint Subproblem Solver

A classical algorithm for solving large-scale trust-region subproblems. It uses conjugate gradient iterations while monitoring negative curvature and boundary conditions, requiring only one matrix-vector product per iteration, making it suitable for large-scale or operator-form Hessians.