/*
  Private data structure used for bqpip method
*/

#pragma once

#include <petsc/private/taoimpl.h>

typedef struct {
  /* Parameters */
  PetscInt predcorr;

  /* Problem variables, vectors and index sets */
  PetscInt n; /* Dimension of the Problem */
  PetscInt m; /* Number of constraints */

  /* Problem statistics */
  PetscReal dinfeas;
  PetscReal pinfeas;
  PetscReal pobj;
  PetscReal dobj;
  PetscReal gap;
  PetscReal rgap;
  PetscReal mu;
  PetscReal pathnorm;
  PetscReal psteplength;
  PetscReal dsteplength;
  PetscReal rnorm;

  /* Variable Vectors */
  Vec G;
  Vec DG;
  Vec T;
  Vec DT;
  Vec Z;
  Vec DZ;
  Vec S;
  Vec DS;
  Vec GZwork;
  Vec TSwork;
  Vec XL, XU;

  /* Work Vectors */
  Vec R3;
  Vec R5;
  Vec HDiag;
  Vec Work;

  Vec DiagAxpy;
  Vec RHS;
  Vec RHS2;

  /* Data */
  Vec       C;
  PetscReal d;
} TAO_BQPIP;