matrix.h

Go to the documentation of this file.

 
#ifndef MATRIX_H
#define MATRIX_H
 
#include "vector.h"
#include <iomanip>
 
enum MATRIX_TRANSPOSE {TRANSPOSE, NO_TRANSPOSE};
 
class Matrix {
    public: 
        /* CONSTRUCTORS*/
        Matrix(); 
        Matrix(uint16_t , uint16_t );
        Matrix(const Matrix & );
 
        /* DESTRUCTOR */
        virtual ~Matrix();
 
        /* ACCESSORS */
        uint16_t get_n_rows() const; 
        uint16_t get_n_cols() const;
        uint32_t get_n_elements() const; 
 
        /* METHODS */
        
        bool is_square() const; 
        Vector row(uint16_t ) const;
        Vector col(uint16_t ) const;
        double norm1() const;
        double normInf() const; 
        double normFrob() const; 
        Matrix transpose() const;
 
        /* OPERATORS */
        double & operator()(uint32_t ) const;
        double & operator()(uint16_t, uint16_t) const;
 
        Matrix & operator=(const Matrix & );
        Matrix & operator+=(const Matrix & ); 
        Matrix & operator-=(const Matrix & ); 
 
        Matrix operator-() const; 
 
        friend bool operator==(const Matrix & , const Matrix & );
 
        friend Matrix operator*(const Matrix & , const Matrix & ); 
        friend Vector operator*(const Vector & , const Matrix & ); 
        friend Vector operator*(const Matrix & , const Vector & ); 
        friend Matrix operator*(const double , const Matrix & ); 
        friend Matrix operator*(const Matrix & , const double );
 
        friend Matrix operator+(const Matrix & , const Matrix & ); 
        friend Matrix operator+(const double , const Matrix & ); 
        friend Matrix operator+(const Matrix & , const double );
 
        friend Matrix operator-(const Matrix & , const Matrix & ); 
        friend Matrix operator-(const Matrix & , const double );
        friend Matrix operator-(const double , const Matrix & ); 
        
        friend Matrix operator/(const Matrix & , const double );
 
        friend std::ostream & operator<<(std::ostream & , const Matrix & );
 
        /* STATIC METHODS */
        static Matrix zeros(uint16_t , uint16_t ); 
        static Matrix ones(uint16_t , uint16_t ); 
        static Matrix rand(uint16_t , uint16_t );
        static Matrix hilbert(uint16_t , uint16_t ); 
        static Matrix vandermonde(const Vector & , uint16_t );
        static Matrix outer(const Vector & , const Vector & );
        static Matrix matmul(const Matrix & , MATRIX_TRANSPOSE , const Matrix & , MATRIX_TRANSPOSE );
 
    protected:
        /* ATTRIBUTES */
        uint16_t n_rows;     
        uint16_t n_cols;     
        uint32_t n_elements; 
        double * data;       
        /* METHODS */
        double & at(uint32_t ) const;
        double & at(uint16_t , uint16_t ) const;
        void show() const;
        void default_data(); 
};
 
#endif