a_mat.h

Go to the documentation of this file.

/*
         Copyright 2010-2016 Pierre SMARS (smars@yuntech.edu.tw)
         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, either version 2 of the License, or
         (at your option) any later version.
 
         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, see <http://www.gnu.org/licenses/>.
 */
#ifndef _VE_MAT_H_
#define _VE_MAT_H_
 
#include <valarray>
#include <iostream>
#include "a_mat_base.h"
 
typedef unsigned int UI;
 
class a_mat;
 
a_mat operator+(a_mat& a, a_mat& b);
a_mat operator-(a_mat& a, a_mat& b);
a_mat operator*(a_mat& a, a_mat& b);
a_mat operator*(a_mat& a, const double v);
a_mat operator*(const double v, a_mat& a);
a_mat operator/(a_mat& a, const double v);
 
class a_mat: public a_mat_base
{
        public:
                a_mat(UI  i,  UI j, double v = 0); //v : value 
                a_mat(int i, int j, double v = 0);
                a_mat(const a_mat& m);
                ~a_mat();
 
                inline UI maxi() const {return maxi_;}
                inline UI maxj() const {return maxj_;}
                inline UI size() const {return maxi_*maxj_;}
 
                inline std::valarray<double>& x() const {return (*x_);}
 
                inline virtual double& operator()(const UI i, const UI j)       {return (*x_)[i*maxj_+j];}
                inline virtual double  operator()(const UI i, const UI j) const {return (*x_)[i*maxj_+j];}
 
                a_mat& operator=(const a_mat& m);
 
                bool operator==(const a_mat& a);
                a_mat& operator-();
                a_mat& operator+=(const a_mat& a);
                a_mat operator+(const a_mat& a);
                a_mat& operator-=(const a_mat& a);
                a_mat operator-(const a_mat& a);
                a_mat& operator*=(double v);
                a_mat operator*(double v);
                a_mat& operator/=(double v);
                a_mat operator/(double v);
 
                a_mat transpose();
                void swapc(UI c1, UI c2); //swap the columns c1 and c2
                void swapr(UI r1, UI r2); //swap the rows r1 and r2
                a_mat sub_matrix(const UI i, const UI j) const; //matrix formed from (*this) by suppressing row i and column j
                //: input/output
                virtual void read(std::istream &i);
                virtual void write(std::ostream &o) const;
 
                //      friend   std::istream&    operator>> (std::istream& i, a_mat& m);
                //      friend   std::ostream&    operator<< (std::ostream& o, const a_mat& m);
 
                class range_error {};
                class compatibility_error {};
 
        protected:
                UI                    maxi_;
                UI                    maxj_;
 
                std::valarray<double> * x_;
};
 
#endif