src/mathlib_src/GMatrix.h

説明を見る。

#ifndef GMATRIX_H
#define GMATRIX_H

#include "GObject.h"

#include <iostream>
#include <sstream>
#include <string>

namespace GMathLib{

class GMatrix : public GObject
{
public:
        GMatrix();


        GMatrix(int mrow, int mcolumn, double* p_data = 0)
         : GObject()
        {
            allocate_initialize(mrow, mcolumn, p_data);
        }

        GMatrix(std::string obj_name, int mrow, int mcolumn, double* p_data = 0)
         : GObject()
        {
            Object_Name(obj_name);
            allocate_initialize(mrow, mcolumn, p_data);
        }

        virtual ~GMatrix();


        GMatrix(const GMatrix& obj);


        virtual int Copy(const GMatrix& obj);


        virtual int Fill(double a){
            if(p_data == 0){
                return 1;
            }

           for(int i=0; i < row; i++){
                for(int j=0; j < column; j++){
                    Set(i, j, a);
                }
            }   

            return 0;
        }       


        int Add(const GMatrix& a, const GMatrix& b);


        int Sub(const GMatrix& a, const GMatrix& b);


        int Multi(const GMatrix& a, const GMatrix& b);


        inline void Scalar_Product(double sc){
            for(int i=0; i < row; i++){
                for(int j=0; j < column; j++){
                    Set(i, j, sc * Get(i, j) ); 
                }
             }
        }


        inline int Row() const{return row;}


        inline int Column() const{return column;}


        inline double* GetData(){return p_data;}


        inline double GetElement(int row_no, int column_no) const{ 
             return p_data[column_no + row_no * column];
        }
        
        inline double Get(int row_no, int column_no) const{ 
             return GetElement(row_no, column_no);
        }


        inline void SetElement(int row_no, int column_no, double value){ 
             p_data[column_no + row_no * column] = value;
        }


        inline void Set(int row_no, int column_no, double value){ 
             SetElement(row_no, column_no, value);
        }

        virtual std::string ToString()
        {
            std::string text = GObject::ToString();
            std::ostringstream oss;

            text.append("Matrix Size: (");
            oss << Row();
            text.append(oss.str()); text.append(","); 
            
            oss.str(std::string());
            oss << Column();
            text.append(oss.str()); text.append(")\n");

            return text;
        }

        void Print() {
             Print(-1);
        }


        void Print(int turn_down, bool scientific_format=false);


        inline void Vectorize(){
            row = row * column;
            column = 1;
        }

        virtual int Reshape(int new_row, int new_column);

        inline double& operator()(int i, int j){
             return p_data[j + i * column];
        }
        
         GMatrix operator=(const GMatrix &a);
        
         GMatrix operator+(const GMatrix &a);
        
         GMatrix operator-(const GMatrix &a);
        
         GMatrix operator*(const GMatrix &a);

         inline friend GMatrix operator*(double k, const GMatrix &a){
         
             GMatrix r(a.Row(), a.Column());
             r.Copy(a);
             r.Scalar_Product(k);

             return r;
         }
         
         inline friend GMatrix operator*(const GMatrix &a, double k){
             return k * a;
         }

protected:
       
       void allocate_initialize(int mrow, int mcolumn, double* p_data = 0);

       double* p_data;                   
private:
        int row;                         
        int column;                      
        bool internal_new_flag;          
};

}
#endif

---