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