---

QpExampleData.h

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

#ifndef QPEXAMPLEDATA_H
#define QPEXAMPLEDATA_H

#include "Data.h"

#include "DoubleMatrixHandle.h"
#include "DenseSymMatrixHandle.h"
#include "OoqpVectorHandle.h"


#include "OoqpVector.h"
#include "DoubleMatrix.h"
class LinearAlgebraPackage;

class QpExampleData : public Data
{
protected:
  QpExampleData() {};

  LinearAlgebraPackage * la;

  int nx;

  int my;

  int mz;

  SymMatrix * Q;

  GenMatrix * A;

  GenMatrix * C;
 
  OoqpVector * c; 

  OoqpVector * b;

  OoqpVector * d;
public:

  QpExampleData(LinearAlgebraPackage * la, int nx_in, int my_in, int mz_in);
  virtual ~QpExampleData();

  virtual double datanorm();

  virtual void datainput() { };
  virtual void datarandom();
  virtual void print();

  void putQIntoAt( SymMatrix& M, int row, int col )
  {
    M.symAtPutSubmatrix( row, col, *Q, 0, 0, nx, nx );
  }

  void putAIntoAt( SymMatrix& M, int row, int col )
  {
    M.symAtPutSubmatrix( row, col, *A, 0, 0, my, nx );
  }

  void putCIntoAt( SymMatrix& M, int row, int col )
  {
    M.symAtPutSubmatrix( row, col, *C, 0, 0, mz, nx );
  }

  void Qmult( double beta,  OoqpVector& y,
              double alpha, OoqpVector& x )
  {
    Q->mult( beta, y, alpha, x );
  }

  void Amult( double beta,  OoqpVector& y,
                      double alpha, OoqpVector& x)
  {
    A->mult( beta, y, alpha, x );
  }

  void Cmult( double beta,  OoqpVector& y,
              double alpha, OoqpVector& x )
  {
    C->mult( beta, y, alpha, x );
  }
  
  void ATransmult( double beta,  OoqpVector& y,
                           double alpha, OoqpVector& x )
  {
    A->transMult( beta, y, alpha, x );
  }

  void CTransmult( double beta,  OoqpVector& y,
                   double alpha, OoqpVector& x )
  {
    C->transMult( beta, y, alpha, x);
  }

  void getg( OoqpVector& c_out ) { c_out.copyFrom( *c ); }

  void getb( OoqpVector& b_out ) { b_out.copyFrom( *b ); }

  void getd( OoqpVector& d_out ) { d_out.copyFrom( *d ); }

  GenMatrix * getC() { IotrAddRef(&C); return C; }
};

#endif

---