1 #include <../src/tao/bound/impls/bqnk/bqnk.h> 2 #include <petscksp.h> 3 4 static PetscErrorCode TaoSetUp_BQNKTR(Tao tao) { 5 KSP ksp; 6 PetscVoidFunction valid; 7 8 PetscFunctionBegin; 9 PetscCall(TaoSetUp_BQNK(tao)); 10 PetscCall(TaoGetKSP(tao, &ksp)); 11 PetscCall(PetscObjectQueryFunction((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(0); 14 } 15 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 22 23 Level: beginner 24 .seealso `TAOBNK`, `TAOBQNKTR`, `TAOBQNKLS` 25 M*/ 26 PETSC_EXTERN PetscErrorCode TaoCreate_BQNKTR(Tao tao) { 27 TAO_BNK *bnk; 28 TAO_BQNK *bqnk; 29 30 PetscFunctionBegin; 31 PetscCall(TaoCreate_BQNK(tao)); 32 tao->ops->setup = TaoSetUp_BQNKTR; 33 bnk = (TAO_BNK *)tao->data; 34 bqnk = (TAO_BQNK *)bnk->ctx; 35 bqnk->solve = TaoSolve_BNTR; 36 PetscFunctionReturn(0); 37 } 38