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