---

QpGenLinsys.h

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

#ifndef QPGENLINSYS
#define QPGENLINSYS

#include "LinearSystem.h"
#include "DoubleMatrixHandle.h"
#include "OoqpVectorHandle.h"
#include "OoqpVector.h"

class Data;
class QpGenData;
class QpGen;
class Variables;
class Residuals;
class DoubleLinearSolver;
class LinearAlgebraPackage;


class QpGenLinsys : public LinearSystem {
protected:

  OoqpVectorHandle nomegaInv;

  OoqpVectorHandle rhs;

  QpGen * factory;

  QpGenLinsys() {};

  int nx, my, mz;

  OoqpVectorHandle dd, dq;

  OoqpVectorHandle ixupp, icupp, ixlow, iclow;

  int nxupp, nxlow, mcupp, mclow;

public:
  QpGenLinsys(  QpGen * factory,
                QpGenData * data,
                LinearAlgebraPackage * la );

  virtual ~QpGenLinsys() {}


  virtual void factor(Data *prob, Variables *vars);

  virtual void solve(Data *prob, Variables *vars, Residuals *res,
                     Variables *step);

  virtual void joinRHS( OoqpVector& rhs,  OoqpVector& rhs1,
                        OoqpVector& rhs2, OoqpVector& rhs3 );

  virtual void separateVars( OoqpVector& vars1, OoqpVector& vars2,
                             OoqpVector& vars3, OoqpVector& vars );

  virtual void solveXYZS( OoqpVector& stepx, OoqpVector& stepy,
                          OoqpVector& stepz, OoqpVector& steps,
                          OoqpVector& ztemp, QpGenData * data );

  virtual void solveCompressed( OoqpVector& rhs ) = 0;

  virtual void putXDiagonal( OoqpVector& xdiag ) = 0;

  virtual void putZDiagonal( OoqpVector& zdiag ) = 0;

  virtual void computeDiagonals( OoqpVector& dd, OoqpVector& omega,
                                 OoqpVector& t,  OoqpVector& lambda,
                                 OoqpVector& u,  OoqpVector& pi,
                                 OoqpVector& v,  OoqpVector& gamma,
                                 OoqpVector& w,  OoqpVector& phi );
};

#endif

---