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

---

MatReplicateMatrix.h

Go to the documentation of this file.

#ifndef Impala_Core_Matrix_MatReplicateMatrix_h
#define Impala_Core_Matrix_MatReplicateMatrix_h

#include "Core/Matrix/MatFunc.h"

namespace Impala
{
namespace Core
{
namespace Matrix
{

template<class ArrayT>
inline ArrayT*
MatReplicateMatrixOneRow(ArrayT* m, int columnCount)
{
    ArrayT* ones = MatCreate<ArrayT>(1, columnCount);
    for(int a = 0; a < columnCount; a++)
    {
        *MatE(ones, 0, a) = 1.0;
    }
    ArrayT* res = MatMul(m, ones);
    delete ones;
    return res;
}

template<class ArrayT>
inline ArrayT*
MatReplicateMatrix(ArrayT* m, int rowCount, int columnCount)
{
    if((rowCount == 1) && (MatNrCol(m) == 1))
    {
        // this operation has poor cache coherence 
        // (the mirror case does not need this optimization)
        return MatReplicateMatrixOneRow(m, columnCount);
    }
    ArrayT* res = MatCreate<ArrayT>(MatNrRow(m) * rowCount, MatNrCol(m) * columnCount);
    for (int a = 0; a < rowCount; a++)
    {
        for (int b = 0; b < columnCount; b++)
        {
            int c = MatNrRow(m) * a;
            int d = MatNrCol(m) * b;
            for (int i=0; i<MatNrRow(m); i++)
            {
                for (int j=0; j<MatNrCol(m); j++)
                {
                    *MatE(res, c+i, d+j) = *MatE(m, i, j);
                }
            }
        }
    }
    return res;
}

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

#endif

---