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