/*
Context for bounded quasi-Newton-Krylov type optimization algorithms
*/

#pragma once

#include <../src/tao/bound/impls/bnk/bnk.h>
#include <../src/ksp/ksp/utils/lmvm/lmvm.h>
#include <../src/ksp/ksp/utils/lmvm/symbrdn/symbrdn.h>

typedef struct {
  PetscErrorCode (*solve)(Tao);
  Mat       B;
  PC        pc;
  PetscBool is_spd;
} TAO_BQNK;

#define BQNK_INIT_CONSTANT  0
#define BQNK_INIT_DIRECTION 1
#define BQNK_INIT_TYPES     2

PETSC_INTERN PetscErrorCode TaoSolve_BQNK(Tao);
PETSC_INTERN PetscErrorCode TaoSetUp_BQNK(Tao);
PETSC_INTERN PetscErrorCode TaoCreate_BQNK(Tao);