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

---

MatTranspose.h

Go to the documentation of this file.

#ifndef Impala_Core_Matrix_MatTranspose_h
#define Impala_Core_Matrix_MatTranspose_h

#include "Core/Matrix/MatFunc.h"

namespace Impala
{
namespace Core
{
namespace Matrix
{


template<class ArrayT>
inline ArrayT*
MatTranspose(ArrayT* m)
{
    ArrayT* res = MatCreate<ArrayT>(MatNrCol(m), MatNrRow(m));
    for (int i=0 ; i<MatNrRow(m) ; i++)
        for (int j=0 ; j<MatNrCol(m) ; j++)
            *MatE(res, j, i) = *MatE(m, i, j);
    return res;
}

template<class ArrayT>
inline void
MatTransposeInplace(ArrayT* m)
{
    ILOG_VAR(Core.Matrix.MatTransposeInplace);
    if(MatNrRow(m) != MatNrCol(m))
        ILOG_ERROR("only works for square matrices!");
    for (int i=0 ; i<MatNrRow(m) ; i++)
        for (int j=i+1 ; j<MatNrCol(m) ; j++)
        {
            typename ArrayT::ArithType temp;
            temp = *MatE(m, j, i);
            *MatE(m, j, i) = *MatE(m, i, j);
            *MatE(m, i, j) = temp;
        }
}

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

#endif

---