Large-Scale Quasi-Newton Trust-Region Methods: High-Accuracy Solvers, Dense Initializations, And Extensions