PrMatSparse.h

//===========================================================================
// GoTools - SINTEF Geometry Tools version 1.1
//
// GoTools module: parametrization
//
// Copyright (C) 2000-2005 SINTEF ICT, Applied Mathematics, Norway.
//
// This program is free software; you can redistribute it and/or          
// modify it under the terms of the GNU General Public License            
// as published by the Free Software Foundation version 2 of the License. 
//
// This program is distributed in the hope that it will be useful,        
// but WITHOUT ANY WARRANTY; without even the implied warranty of         
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the          
// GNU General Public License for more details.                           
//
// You should have received a copy of the GNU General Public License      
// along with this program; if not, write to the Free Software            
// Foundation, Inc.,                                                      
// 59 Temple Place - Suite 330,                                           
// Boston, MA  02111-1307, USA.                                           
//
// Contact information: E-mail: tor.dokken@sintef.no                      
// SINTEF ICT, Department of Applied Mathematics,                         
// P.O. Box 124 Blindern,                                                 
// 0314 Oslo, Norway.                                                     
//
// Other licenses are also available for this software, notably licenses
// for:
// - Building commercial software.                                        
// - Building software whose source code you wish to keep private.        
//===========================================================================
#ifndef PRMATSPARSE_H
#define PRMATSPARSE_H

#include "PrMatrix.h"
#include "errormacros.h"
#include <vector>

/*<PrMatSparse-syntax: */

class PrMatSparse : public PrMatrix
{
private:
  int m_, n_; // matrix dimensions
  int p_;     // number of non-zeros
  std::vector<int> irow_; // the indexes in a_ and jcol_ of the
                     // first non-zeros of the m rows
  std::vector<int> jcol_; // the p indexes j of the non-zero elements
  std::vector<double> a_;   // the p non-zero elements
  

public:

  virtual int rows() const;
  virtual int colmns() const;
  virtual double operator () (int i, int j) const;
  virtual void prod(const PrVec& x, PrVec& y) const;
  virtual void print(std::ostream& os);
  virtual void read(std::istream& is);
  virtual ~PrMatSparse();




  void matProd(PrMatSparse& B, PrMatSparse& C) const;
  void scalProd(double d);
  void printFull(std::ostream& os);

  PrMatSparse() : m_(0), n_(0), p_(0) {}
  PrMatSparse(int m, int n, int num_nonzero);
  PrMatSparse(int m, int n, int num_nonzero,
              const int* irow, const int* jcol, const double* data);
  void redim(int m, int n, int num_nonzero);
  void setToMatrix(const PrMatrix& m, double tol = 0.0);

    int& irow(int k);        
    const int& irow(int k) const;  
    int& jcol(int k);        
    const int& jcol(int k) const; 
    double& operator () (int k);
    const double& operator () (int k) const;
};


/*>PrMatSparse-syntax: */

//-----------------------------------------------------------------------------
inline void PrMatSparse::scalProd(double d)
//-----------------------------------------------------------------------------
{
    for (std::vector<double>::iterator it = a_.begin(); it != a_.end(); ++it) {
        (*it) *= d;
    }
}


//-----------------------------------------------------------------------------
inline const int& PrMatSparse::irow(int k) const
//-----------------------------------------------------------------------------
{
  return irow_[k];
}

//-----------------------------------------------------------------------------
inline int& PrMatSparse::irow(int k)
//-----------------------------------------------------------------------------
{
  return irow_[k];
}

//-----------------------------------------------------------------------------
inline const int& PrMatSparse::jcol(int k) const
//-----------------------------------------------------------------------------
{
  return jcol_[k];
}

//-----------------------------------------------------------------------------
inline int& PrMatSparse::jcol(int k)
//-----------------------------------------------------------------------------
{
  return jcol_[k];
}

//-----------------------------------------------------------------------------
inline const double& PrMatSparse::operator() (int k) const
//-----------------------------------------------------------------------------
{
  return a_[k];
}

//-----------------------------------------------------------------------------
inline double& PrMatSparse::operator() (int k)
//-----------------------------------------------------------------------------
{
  return a_[k];
}

/*Class:PrMatSparse

Name:              PrMatSparse
Syntax:            @PrMatSparse-syntax
Keywords:
Description:       This class implements a matrix
Member functions:

Constructors:
Files:
Example:

See also:
Developed by:      SINTEF Applied Mathematics, Oslo, Norway
Author:            Michael Floater, SINTEF
Date:              Dec. 98
*/

#endif // PRMATSPARSE_H

---