1e0ed867bSAlp Dener #include <../src/tao/bound/impls/bqnk/bqnk.h> 2e0ed867bSAlp Dener 3*3850be85SAlp Dener /*MC 4*3850be85SAlp Dener TAOBQNKTL - Bounded Quasi-Newton-Krylov Trust-region with Line-search fallback, for nonlinear 5*3850be85SAlp Dener minimization with bound constraints. This method approximates the Hessian-vector 6*3850be85SAlp Dener product using a limited-memory quasi-Newton formula, and iteratively inverts the 7*3850be85SAlp Dener Hessian with a Krylov solver. The quasi-Newton matrix and its settings can be 8*3850be85SAlp Dener accessed via the prefix `-tao_bqnk_` 9*3850be85SAlp Dener 10*3850be85SAlp Dener Options Database Keys: 11*3850be85SAlp Dener + -tao_bqnk_max_cg_its - maximum number of bounded conjugate-gradient iterations taken in each Newton loop 12*3850be85SAlp Dener . -tao_bqnk_init_type - trust radius initialization method ("constant", "direction", "interpolation") 13*3850be85SAlp Dener . -tao_bqnk_update_type - trust radius update method ("step", "direction", "interpolation") 14*3850be85SAlp Dener - -tao_bqnk_as_type - active-set estimation method ("none", "bertsekas") 15*3850be85SAlp Dener 16*3850be85SAlp Dener Level: beginner 17*3850be85SAlp Dener M*/ 18e0ed867bSAlp Dener PETSC_EXTERN PetscErrorCode TaoCreate_BQNKTL(Tao tao) 19e0ed867bSAlp Dener { 20e0ed867bSAlp Dener PetscErrorCode ierr; 21e0ed867bSAlp Dener 22e0ed867bSAlp Dener PetscFunctionBegin; 23e0ed867bSAlp Dener ierr = TaoCreate_BQNK(tao);CHKERRQ(ierr); 24e0ed867bSAlp Dener tao->ops->solve = TaoSolve_BNTL; 25e0ed867bSAlp Dener PetscFunctionReturn(0); 26e0ed867bSAlp Dener } 27