D:/User/Nacho/Mecatronica/C2-CSSD/ProgramacionCpp/Ej1Matrix/MyMatrix.h

Go to the documentation of this file.

#ifndef _INC_MYMATRIX_H
#define _INC_MYMATRIX_H
 
#include <vector>
#include <string>
#include <stdexcept>
#include <math.h>
 
typedef double dataType;    
 
class MyMatrix {
public:
   MyMatrix();
 
   MyMatrix( int nrows, int ncols,dataType value=NAN );
 
   MyMatrix(const std::string& dataTxt );
 
   MyMatrix( const std::vector<dataType>& data, bool rowVector=true );
 
   MyMatrix( const MyMatrix& other );
 
   ~MyMatrix();
 
   inline bool IsValid() const {
       return nRows>0 && nCols>0 && data!=nullptr;
   }
 
   inline dataType GetAt( int iRow, int iCol ) const {
        return data[ GetOffset(iRow, iCol ) ];
    }
 
    inline void SetAt( int iRow, int iCol, dataType value ) {
        data[ GetOffset(iRow, iCol) ] = value;
    }
 
    inline int GetRows() const {
        return nRows;
    }
 
    inline int GetCols() const {
        return nCols;
    }
 
    inline dataType& operator()(int iRow, int iCol ) {
        return data[ GetOffset(iRow,iCol) ];
    }
 
        inline const dataType& operator()(int iRow, int iCol ) const {
        return data[ GetOffset(iRow,iCol) ];
    }
 
    MyMatrix& operator=( const MyMatrix& other );
 
    MyMatrix operator+(dataType scalar ) const;
 
    MyMatrix operator+(const MyMatrix& other ) const;
 
    MyMatrix operator-(dataType scalar ) const;
 
    MyMatrix operator-(const MyMatrix& other ) const;
 
    MyMatrix operator*(dataType scalar ) const;
 
 
    MyMatrix operator*(const MyMatrix& other ) const;
 
    MyMatrix operator/(dataType scalar ) const;
 
    MyMatrix ProdElByEl(const MyMatrix& other) const;
 
    MyMatrix DivElByEl(const MyMatrix& other) const;
 
    MyMatrix Transpose() const;
 
 
    bool IsEqual(const MyMatrix& other, double tolerance=1.0e-6 ) const;
 
    static MyMatrix Eye(int size);
 
 
    MyMatrix Adjunct(int i_row,int i_col) const;
// TODO    int Rank() const;
    dataType Determinant() const;
 
    MyMatrix GetRow(int i_row) const;
    MyMatrix GetCol(int i_col) const;
    MyMatrix GetSubmatrix(int i_row,int n_rows,int i_col,int n_cols) const;
 
private:
    // Documentation for these functions and data should not be in doxygen format, so it doesn't appear in the output
    int nRows=-1,nCols=-1;
    dataType* data=nullptr;
 
    void AssignMemo(int i_rows,int i_cols);
    void FreeMemo();
    inline int GetOffset(int i_row,int i_col) const
    {
        if (i_row>=0 && i_row<nRows && i_col>=0 && i_col<=nCols)
            return i_row*nCols+i_col;
        throw std::runtime_error("Matrix error access");
    }
 
    void CopyFrom(const MyMatrix& other);
    bool ReadFromString(const char *pt, int *ptRows = nullptr, int *ptCols = nullptr);
};
 
 
std::ostream& SetPrecisionAndFormat(std::ostream& out,int nDigits,int nDecimals);
 
std::ostream& operator<<(std::ostream& out,const MyMatrix& m );
 
std::istream& operator>>(std::istream& in,MyMatrix& m );
 
inline MyMatrix operator*(dataType scalar, const MyMatrix& m ) {
    return m*scalar;
}
 
inline MyMatrix operator+(dataType scalar, const MyMatrix& m ) {
    return m+scalar;
}
 
inline MyMatrix operator-(dataType scalar, const MyMatrix& m ) {
    return m*(-1)+scalar;
}
 
#endif // _INC_MYMATRIX_H