---

DoubleMatrix.h

Go to the documentation of this file.

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

#ifndef DOUBLEMATRIX_H
#define DOUBLEMATRIX_H

#include <string.h>
#include <memory.h>
#include <iostream.h>
#include "OoqpVectorHandle.h"
#include "IotrRefCount.h"

class DoubleLinearSolver;

class DoubleStorage : public IotrRefCount {
public:
  DoubleStorage() { };

  virtual void atPutDense( int row, int col, double * A, int lda,
                           int rowExtent, int colExtent ) = 0;
  virtual void fromGetDense( int row, int col, double * A, int lda,
                             int rowExtent, int colExtent ) = 0;
  virtual void atPutSpRow( int row, double A[], int lenA, int jcolA[],
                           int& info ) = 0;

  virtual void fromGetSpRow( int row, int col,
                             double A[], int lenA, int jcolA[], int& nnz,
                             int colExtent, int& info ) = 0;

  virtual void getSize( int& m, int& n )  = 0;

  virtual void getDiagonal( OoqpVector& vec ) = 0;
  virtual void setToDiagonal( OoqpVector& vec ) = 0;

  virtual void atPutDiagonal( int idiag, OoqpVector& x ) = 0;
  virtual void fromGetDiagonal( int idiag, OoqpVector& x ) = 0;
  virtual void SymmetricScale ( OoqpVector& vec ) = 0;
  virtual void ColumnScale ( OoqpVector& vec ) = 0;
  virtual void RowScale ( OoqpVector& vec ) = 0;
  virtual void scalarMult( double num) = 0;
  virtual ~DoubleStorage() {};
};

class DoubleMatrix : public IotrRefCount {
public:
  DoubleMatrix() {}

  virtual int isKindOf( int matrixType ) = 0;

  virtual void fromGetDense( int row, int col, double * A, int lda,
                             int rowExtent, int colExtent ) = 0;

  virtual void fromGetSpRow( int row, int col,
                             double A[], int lenA, int jcolA[], int& nnz,
                             int colExtent, int& info ) = 0;

  virtual void putSparseTriple( int irow[], int len, int jcol[], double A[], 
                                int& info ) = 0;

  virtual void mult ( double beta,  OoqpVector& y,
                      double alpha, OoqpVector& x ) = 0;

  virtual void transMult ( double beta,   OoqpVector& y,
                           double alpha,  OoqpVector& x ) = 0;

  virtual double abmaxnorm() = 0;

  virtual void writeToStream(ostream& out) const = 0;

  virtual void getDiagonal( OoqpVector& vec ) = 0;
  virtual void setToDiagonal( OoqpVector& vec ) = 0;

  virtual void atPutDiagonal( int idiag, OoqpVector& x ) = 0;
  virtual void fromGetDiagonal( int idiag, OoqpVector& x )= 0;

  virtual void SymmetricScale ( OoqpVector& vec ) = 0;
  virtual void ColumnScale ( OoqpVector& vec ) = 0;
  virtual void RowScale ( OoqpVector& vec ) = 0;
  virtual void scalarMult( double num) = 0;

  virtual void getSize( int& m, int& n )  = 0;

  virtual ~DoubleMatrix() {};
};

class SymMatrix    : public DoubleMatrix {
public:

  virtual void symAtPutSubmatrix( int destRow, int destCol,
                                  DoubleMatrix& M,
                                  int srcRow, int srcCol,
                                  int rowExtent, int colExtent ) = 0;
  
  virtual void randomizePSD(double * seed) = 0;

  virtual void symAtPutSpRow( int col, double A[], int lenA, int irowA[],
                              int& info ) = 0;
  virtual int size() = 0;
};

class GenMatrix    : public DoubleMatrix {
public:
  virtual void atPutSubmatrix( int destRow, int destCol, 
                               DoubleMatrix& M,
                               int srcRow, int srcCol,
                               int rowExtent, int colExtent ) = 0;

  virtual void atPutDense( int row, int col, double * A, int lda,
                           int rowExtent, int colExtent ) = 0;

  virtual void atPutSpRow( int col, double A[], int lenA, int jcolA[],
                           int& info ) = 0;

  virtual void randomize(double alpha, double beta, double * seed) = 0;
};

#endif

---