1: #include <../src/tao/bound/impls/bqnk/bqnk.h>

  3: static PetscErrorCode TaoSetUp_BQNKTL(Tao tao)
  4: {
  5:   KSP       ksp;
  6:   PetscBool valid;

  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: }

 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`

 23:   Level: beginner

 25:   Notes:
 26:   The base class for this method is `TAOBNK`

 28:   The various algorithmic factors can only be supplied via the options database

 30: .seealso: `Tao`, `TAOBNK`, `TAONLS`, `TAONTL`, `TAONM`, `TaoType`, `TaoCreate()`, `TAOBQNKTR`, `TAOBQNKLS`
 31: M*/
 32: PETSC_EXTERN PetscErrorCode TaoCreate_BQNKTL(Tao tao)
 33: {
 34:   TAO_BNK  *bnk;
 35:   TAO_BQNK *bqnk;

 37:   PetscFunctionBegin;
 38:   PetscCall(TaoCreate_BQNK(tao));
 39:   tao->ops->setup = TaoSetUp_BQNKTL;
 40:   bnk             = (TAO_BNK *)tao->data;
 41:   bqnk            = (TAO_BQNK *)bnk->ctx;
 42:   bqnk->solve     = TaoSolve_BNTL;
 43:   PetscFunctionReturn(PETSC_SUCCESS);
 44: }