1: #include <../src/tao/bound/impls/bqnk/bqnk.h> 2: #include <petscksp.h> 4: static PetscErrorCode TaoSetUp_BQNKTR(Tao tao) 5: { 6: KSP ksp; 7: PetscVoidFunction valid; 9: TaoSetUp_BQNK(tao); 10: TaoGetKSP(tao, &ksp); 11: PetscObjectQueryFunction((PetscObject)ksp, "KSPCGSetRadius_C", &valid); 13: return 0; 14: } 16: /*MC 17: TAOBQNKTR - Bounded Quasi-Newton-Krylov Trust Region method for nonlinear minimization with 18: bound constraints. This method approximates the Hessian-vector product using a 19: limited-memory quasi-Newton formula, and iteratively inverts the Hessian with a 20: Krylov solver. The quasi-Newton matrix and its settings can be accessed via the 21: prefix `-tao_bqnk_`. For options database, see TAOBNK 23: Level: beginner 25: .seealso: `Tao`, `TaoType`, `TAOBNK`, `TAOBQNKTR`, `TAOBQNKLS` 26: M*/ 27: PETSC_EXTERN PetscErrorCode TaoCreate_BQNKTR(Tao tao) 28: { 29: TAO_BNK *bnk; 30: TAO_BQNK *bqnk; 32: TaoCreate_BQNK(tao); 33: tao->ops->setup = TaoSetUp_BQNKTR; 34: bnk = (TAO_BNK *)tao->data; 35: bqnk = (TAO_BQNK *)bnk->ctx; 36: bqnk->solve = TaoSolve_BNTR; 37: return 0; 38: }