---

QpGenVars.h

/* OOQP                                                               *
 * Authors: E. Michael Gertz, Stephen J. Wright                       *
 * (C) 2001 University of Chicago. See Copyright Notification in OOQP */

#ifndef QPGENVARS
#define QPGENVARS

#include "Variables.h"
#include "OoqpVectorHandle.h"

class QpGen;
class QpGenData;
class LinearAlgebraPackage;
class MpsReader;

#ifdef TESTING
class QpGenVarsTester;
#endif

class QpGenVars : public Variables {
#ifdef TESTING
  friend QpGenVarsTester;
#endif
protected:
  int nx, nxupp, nxlow;
  int my;
  int mz, mcupp, mclow;

  OoqpVectorHandle ixlow;
  OoqpVectorHandle ixupp;
  OoqpVectorHandle icupp;
  OoqpVectorHandle iclow;

public:  

  OoqpVectorHandle x;
  OoqpVectorHandle s;
  OoqpVectorHandle y;
  OoqpVectorHandle z;

  OoqpVectorHandle v;
  OoqpVectorHandle phi;

  OoqpVectorHandle w;
  OoqpVectorHandle gamma;

  OoqpVectorHandle t;
  OoqpVectorHandle lambda;
  
  OoqpVectorHandle u;
  OoqpVectorHandle pi;

  QpGenVars( OoqpVector * x_in, OoqpVector * s_in,
             OoqpVector * y_in, OoqpVector * z_in,
             OoqpVector * v_in, OoqpVector * gamma_in,
             OoqpVector * w_in, OoqpVector * phi_in,
             OoqpVector * t_in, OoqpVector * lambda_in,
             OoqpVector * u_in, OoqpVector * pi_in,
             OoqpVector * ixlow_in, OoqpVector * ixupp_in,
             OoqpVector * iclow_in, OoqpVector * icupp_in );

  QpGenVars( LinearAlgebraPackage * la,
             int nx_, int my_, int mz_,
             OoqpVector * ixlow, OoqpVector * ixupp,
             OoqpVector * iclow, OoqpVector * icupp );

  virtual ~QpGenVars();
  
  virtual double mu();

  virtual double mustep(Variables *step_in, double alpha);

  virtual void saxpy( Variables *b, double alpha );
  virtual void negate();
  
  virtual double stepbound( Variables *b );

  virtual double findBlocking( Variables * step, 
                               double & primalValue,
                               double & primalStep,
                               double & dualValue,
                               double & dualStep,
                               int& firstOrSecond );

  virtual void interiorPoint( double alpha, double beta );

  virtual void shiftBoundVariables( double alpha, double beta );

  virtual int isInteriorPoint();

  virtual double violation();

  virtual void print();
  virtual void printSolution( MpsReader * reader, QpGenData * prob,
                              int& iErr );

  virtual void unscaleSolution( QpGenData * data);
  virtual void unscaleBounds  ( QpGenData * data);

  virtual int  validNonZeroPattern();
  
  virtual void copy(Variables *b);

  virtual double onenorm();
  virtual double infnorm();
};

enum { no_block = 0,
       t_block = 1,
       lambda_block = 2,
       u_block = 3,
       pi_block = 4,
       v_block = 5,
       gamma_block = 6,
       w_block = 7,
       phi_block = 8
};  

#endif

---