---

DenseGenMatrix.hh

Go to the documentation of this file.

#ifndef DENSEGENMATRIX
#define DENSEGENMATRIX

#include "GenMatrix.hh"
#include "IotrRefCount.hh" 
#include "SimpleArrayData.hh"
#include "ElementGatherSeq.hh"
#include "IotrHandle.hh"

class DenseGenMatrix : public GenMatrix {
protected:
  SimpleArrayDataHandle data;
  double ** mMat;
  int mLeadingDimension;
  int mRows;
  int mCols;
  /* my changes */

  // Disable C++ value semantics (is this necessary?)
  DenseGenMatrix() : data(), mMat(0),
                     mLeadingDimension(0), mRows(0), mCols(0) {};
public:
  DenseGenMatrix( int m, int n);
  DenseGenMatrix( double * v, int m, int n);
  virtual ~DenseGenMatrix();

  double * operator[]( int i ) { return mMat[i]; }
  const double * operator[]( int i ) const { return mMat[i]; }
 
  double * elements() { return mMat[0]; };
  const double * elements() const { return mMat[0]; };
  double **Mat() { return mMat; };
  const double * const *Mat () const { return mMat; };
  int leadingDimension() const { return mLeadingDimension; };

  virtual int rows() const;
  virtual int columns() const;
  virtual void clearRange( int istart, int jstart,
                           int iextent, int jextent );
  virtual void writefToStream( std::ostream& out,
                               const char format[] ) const;
  virtual void mult( double alpha, const IotrVector & x,
                     double beta,  IotrVector & y ) const;

  virtual void filteredMult( double alpha, const ElementGather & filter,
                             const IotrVector & x,
                             double beta,  IotrVector & y ) const;
  virtual void transposeMult( double alpha, const IotrVector & x,
                              double beta,  IotrVector & y ) const;
  virtual void filteredTransposeMult( double alpha, const IotrVector & x,
                                      double beta,  const ElementGather& filt,
                                      IotrVector & y ) const;
//   virtual void multSelect( double alpha, IotrVector & x,
//                         double beta, IotrVector & y,
//                         ElementGatherSeq & iequality);
  virtual void copyFromSparseCols( const int irow[], const int kcol[],
                                   const double M[] );
  virtual void copyToSparseCols( int irow[], int len, 
                                 int kcol[], double M[],
                                 int & morei ) const;
  virtual void inspect( const char * name = 0,
                        const char * filename = 0) const;
  virtual void inspectAsMatlab( const char * name = 0,
                                const char * filename = 0) const;
};

typedef Handle<DenseGenMatrix> DenseGenMatrixHandle;

#endif

---