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_2_8_1-0-gea32989045 $
 * @(#)$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 -------------------------------*/