|| || || || || ||

---

MatEye.h

Go to the documentation of this file.

#ifndef Impala_Core_Matrix_MatEye_h
#define Impala_Core_Matrix_MatEye_h

#include "Core/Matrix/MatFunc.h"

namespace Impala
{
namespace Core
{
namespace Matrix
{

template<class ArrayT>
inline ArrayT*
MatEye(int n)
{
    ArrayT* res = MatCreate<ArrayT>(n, n);
    for(int i = 0; i < n; i++)
    {
        for(int j = 0; j < n; j++)
        {
            if(i == j)
                *MatE(res, i, j) = 1.0;
            else
                *MatE(res, i, j) = 0.0;
        }
    }
    return res;
}

} // namespace Matrix
} // namespace Core
} // namespace Impala

#endif

---