1: #include <../src/tao/bound/impls/bqnk/bqnk.h>
2: #include <petscksp.h>
4: static PetscErrorCode TaoSetUp_BQNKTR(Tao tao)
5: {
6: KSP ksp;
7: PetscVoidFunction valid;
9: TaoSetUp_BQNK(tao);
10: TaoGetKSP(tao, &ksp);
11: PetscObjectQueryFunction((PetscObject)ksp, "KSPCGSetRadius_C", &valid);
13: return 0;
14: }
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
23: Level: beginner
24: .seealso `TAOBNK`, `TAOBQNKTR`, `TAOBQNKLS`
25: M*/
26: PETSC_EXTERN PetscErrorCode TaoCreate_BQNKTR(Tao tao)
27: {
28: TAO_BNK *bnk;
29: TAO_BQNK *bqnk;
31: 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: return 0;
37: }