ms_matrix.hpp

/*
 * Copyright (C) 2002-2013, Matrix Science Limited, London, UK.
 * All rights reserved.
 *
 * @(#)$Source: parser/inc/ms_matrix.hpp $
 * @(#)$Revision: 436aa2996c6a3d2d2f0cdb88fab902c2780ae29a | MSPARSER_REL_3_0_0-2024-09-24-0-g93ebaeb4f4 $
 * @(#)$Date: 2019-05-16 11:53:33 +0100 $
 */
 
#ifndef MS_MATRIX_HPP
#define MS_MATRIX_HPP
 
 
#include <vector>
 
namespace matrix_science
{
 
    class MS_MASCOTRESFILE_API ms_matrix
    {
        public:
 
            typedef std::vector<double> column_vector_t;
            typedef std::vector<column_vector_t> matrix_t;
 
            static ms_matrix diagonalMatrix(size_t s, double v);
 
            static ms_matrix identityMatrix(size_t s);
 
            static ms_matrix rowMatrix(const std::vector<double> & v);
 
            static ms_matrix columnMatrix(const std::vector<double> & v);
 
            ms_matrix();
 
            ms_matrix(const matrix_t & m);
 
            ms_matrix(size_t si, size_t sj, double v = 0.0);
 
            size_t sizei() const;
 
            size_t sizej() const;
 
#ifndef SWIG
            column_vector_t & operator[](size_t i);
 
            const column_vector_t & operator[](size_t i) const;
#endif
            double getCR(size_t i, size_t j) const;
            void setCR(size_t i, size_t j, double value);
 
            double getRC(size_t i, size_t j) const { return getCR(j, i); }
            void setRC(size_t i, size_t j, double value) {  setCR(j, i, value); }
 
            const matrix_t & data() const;
#ifndef SWIG
            bool operator==(const ms_matrix & m) const;
 
            bool operator!=(const ms_matrix & m) const;
 
            ms_matrix operator*(const ms_matrix & m) const;
#endif
            bool isZero() const;
 
            void resize(size_t si, size_t sj, double v = 0.0);
 
            double multiplyByRow(const std::vector<double> & column, size_t row) const;
 
            void setColumnDown(size_t row, size_t col, const std::vector<double> & vec);
 
            void setColumnUp(size_t row, size_t col, const std::vector<double> & vec);
 
            void transpose();
 
            void logContents(ms_errs & errs, const char * description) const;
 
        private:
            matrix_t matrix_;
    };
 // end of tools_group
} // namespace matrix_science
 
#endif // MS_MATRIX_HPP
 
/*------------------------------- End of File -------------------------------*/