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