1 #include <../src/tao/bound/impls/bqnk/bqnk.h> 2 3 static PetscErrorCode TaoSetUp_BQNKTL(Tao tao) 4 { 5 KSP ksp; 6 PetscBool valid; 7 8 PetscFunctionBegin; 9 PetscCall(TaoSetUp_BQNK(tao)); 10 PetscCall(TaoGetKSP(tao, &ksp)); 11 PetscCall(PetscObjectHasFunction((PetscObject)ksp, "KSPCGSetRadius_C", &valid)); 12 PetscCheck(valid, PetscObjectComm((PetscObject)tao), PETSC_ERR_SUP, "Not for KSP type %s. Must use a trust-region CG method for KSP (e.g. KSPNASH, KSPSTCG, KSPGLTR)", ((PetscObject)ksp)->type_name); 13 PetscFunctionReturn(PETSC_SUCCESS); 14 } 15 16 /*MC 17 TAOBQNKTL - Bounded Quasi-Newton-Krylov Trust-region with Line-search fallback, for nonlinear 18 minimization with bound constraints. This method approximates the Hessian-vector 19 product using a limited-memory quasi-Newton formula, and iteratively inverts the 20 Hessian with a Krylov solver. The quasi-Newton matrix and its settings can be 21 accessed via the prefix `-tao_bqnk_`. For options database, see `TAOBNK` 22 23 Level: beginner 24 25 .seealso: `Tao`, `TaoType`, `TAOBNK`, `TAOBQNKTR`, `TAOBQNKLS` 26 M*/ 27 PETSC_EXTERN PetscErrorCode TaoCreate_BQNKTL(Tao tao) 28 { 29 TAO_BNK *bnk; 30 TAO_BQNK *bqnk; 31 32 PetscFunctionBegin; 33 PetscCall(TaoCreate_BQNK(tao)); 34 tao->ops->setup = TaoSetUp_BQNKTL; 35 bnk = (TAO_BNK *)tao->data; 36 bqnk = (TAO_BQNK *)bnk->ctx; 37 bqnk->solve = TaoSolve_BNTL; 38 PetscFunctionReturn(PETSC_SUCCESS); 39 } 40