#pragma once

#include <petsc/private/taoimpl.h>
#include <petscmath.h>

#define BMRM_INFTY 1e30 /* single precision: ~\pm 10^{38.53}; PetscReal precision: ~\pm 10^{308.25} */
#define ALPHA_MIN  1e-10
#define ALPHA_MAX  1e10
#define EPS_SV     1e-15
#define EPS        1e-20
#define TOL_LAM    1e-15
#define TOL_R      1e-10
#define INCRE_DIM  1000

/* Context for BMRM solver */
typedef struct {
  VecScatter scatter; /* Scatter context  */
  Vec        local_w;
  PetscReal  lambda;
} TAO_BMRM;

typedef struct Vec_Chain {
  Vec               V;
  struct Vec_Chain *next;
} Vec_Chain;

/* Context for Dai-Fletcher solver */
typedef struct {
  PetscInt   maxProjIter;
  PetscInt   maxPGMIter;
  PetscInt  *ipt, *ipt2, *uv;
  PetscReal *g, *y, *tempv, *d, *Qd, *t, *xplus, *tplus, *sk, *yk;

  PetscInt dim;

  PetscInt cur_num_cp;

  /* Variables (i.e. Lagrangian multipliers) */
  PetscReal *x;

  /* Linear part of the objective function  */
  PetscReal *f;

  /* Hessian of the QP */
  PetscReal **Q;

  /* Constraint matrix  */
  PetscReal *a;

  /* RHS of the equality constraint */
  PetscReal b;

  /* Lower bound vector for the variables */
  PetscReal *l;

  /* Upper bound vector for the variables */
  PetscReal *u;

  /* Tolerance for optimization error */
  PetscReal tol;
} TAO_DF;