---

SparseSymMatrix.hh

Go to the documentation of this file.

#ifndef SPARSESYMMATRIX
#define SPARSESYMMATRIX

#include "SymMatrix.hh"
#include "IotrRefCount.hh"
#include "SimpleArrayData.hh"
#include "SimpleIntData.hh"

class SparseSymMatrix : public SymMatrix {
protected:
  SimpleArrayDataHandle mDataArray;
  SimpleIntDataHandle   mIrowArray;
  SimpleIntDataHandle   mKcolArray;
  double * mElements;
  int    * mIrow;
  int    * mKcol;
  int mLeadingDimension;
  /* NXN matrix with mz nonzeros */
  int mMz;
public:
  SparseSymMatrix( int n, int mz );
  SparseSymMatrix( double * v, int * irow, int * kcol,
                   int n, int mz );
  virtual ~SparseSymMatrix();
  bool operator==( const SparseSymMatrix & MM ) const;
  bool operator!=( const SparseSymMatrix & MM ) const
  {
    return !(*this == MM);
  }

  double * elements() { return mElements; };
  const double * elements() const { return mElements; };
  int * irow() { return mIrow; };
  const int * irow() const { return mIrow; };
  int * kcol() { return mKcol; };
  const int * kcol() const { return mKcol; };
  int leadingDimension() const { return mLeadingDimension; };

  virtual int rows() const { return mN; };
  virtual int columns() const { return mN; };
  int maxNonZeros() const { return mMz; };
  virtual void setToZero();
  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 symMaskedMult( double alpha, const IotrMask & mask,
                              const IotrVector & x,
                              double beta,  IotrVector & y ) const;
  virtual void symCopyFromSparseCols( const int irow[], const int kcol[],
                                      const double M[], int & has_upper );
  virtual void symCopyToSparseCols( int irow[], int len, 
                                    int kcol[],
                                    double M[], int & morei ) const;
  
  virtual void copyFromCuteSymmetric( double * M, int * irow,
                                      int * jcol, int nz );
  virtual void copyFromDenseSymmetric( const double *v, int vlength );
//   virtual void atPutSpCol( int col, const double A[], int lenA,
//                         const int irowA[], int& info );
//   virtual void shiftCols( int col, int shift, int& info );
//   virtual void fromGetSpRow( int row, int col,
//                           double A[], int lenA,
//                           int jcolA[], int& nnz,
//                           int colExtent, int& info ) const;
  virtual void inspect(const char * name_in = 0,
                       const char * filename = 0 ) const;
  virtual void inspectAsMatlab(const char * name_in = 0,
                               const char * filename = 0 ) const;
  virtual double lowerBoundEigenvalue() const;
};

typedef Handle<SparseSymMatrix> SparseSymMatrixHandle;

#endif

---