src/mathlib_src/GMatrix_Util.h

説明を見る。

#ifndef GMatrix_UTIL_H
#define GMatrix_UTIL_H

#include "GObject.h"
#include "GMatrix.h"
#include "io/GError_Output.h"
#include <string>

namespace GMathLib
{

namespace GMatrix_Util{

inline void Identity_Matrix(GMatrix& mx){
    int column = mx.Column();
    int row = mx.Row();
    
    //正方行列でない場合はエラーを出力する
    if( column != row ){
        const std::string context = "A given matrix is not a square matrix.";
        IO::GError_Output::Puts("GMatrix_Util::Identity_Matrix", context);   

    }else{
        //行列を単位行列に書き換える
        for(int i = 0; i < row; i++){
            for(int j=0; j < column; j++){
                if(i == j){
                    mx.SetElement(i, j, 1);
                }else{
                    mx.SetElement(i, j, 0);
                }
            }
        }

    }
}

} //end namespace GMatrix_Util
} //end namespace GMathLib

#endif

---