1e0ed867bSAlp Dener #include <../src/tao/bound/impls/bqnk/bqnk.h> 2e0ed867bSAlp Dener 3*3850be85SAlp Dener /*MC 4*3850be85SAlp Dener TAOBQNKLS - Bounded Quasi-Newton-Krylov Line Search method for nonlinear minimization with 5*3850be85SAlp Dener bound constraints. This method approximates the Hessian-vector product using a 6*3850be85SAlp Dener limited-memory quasi-Newton formula, and iteratively inverts the Hessian with a 7*3850be85SAlp Dener Krylov solver. The quasi-Newton matrix and its settings can be accessed via the 8*3850be85SAlp Dener 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_BQNKLS(Tao tao) 19e0ed867bSAlp Dener { 20e0ed867bSAlp Dener TAO_BNK *bnk; 21e0ed867bSAlp Dener PetscErrorCode ierr; 22e0ed867bSAlp Dener 23e0ed867bSAlp Dener PetscFunctionBegin; 24e0ed867bSAlp Dener ierr = TaoCreate_BQNK(tao);CHKERRQ(ierr); 25e0ed867bSAlp Dener tao->ops->solve = TaoSolve_BNLS; 26e0ed867bSAlp Dener bnk = (TAO_BNK*)tao->data; 27e0ed867bSAlp Dener bnk->update_type = BNK_UPDATE_STEP; 28e0ed867bSAlp Dener PetscFunctionReturn(0); 29e0ed867bSAlp Dener } 30