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

---

SetRow.h

Go to the documentation of this file.

#ifndef Impala_Core_Matrix_SetRow_h
#define Impala_Core_Matrix_SetRow_h

#include "Basis/ILog.h"
#include "Core/Matrix/MatFunc.h"
#include "Core/Vector/VectorTem.h"

namespace Impala
{
namespace Core
{
namespace Matrix
{

template<class MatT>
inline void
SetRow(MatT* m, int index, const Vector::VectorTem<Real64>& row)
{
    ILOG_VAR(Impala.Core.Matrix.SetRow);
    int length = MatNrCol(m);
    if (length != row.Size())
        ILOG_ERROR_NODE("nonconformant operands");
    if (index < 0 || index >= MatNrRow(m))
        ILOG_ERROR_NODE("faulty row index: " << index);
    typename MatT::StorType* rowPtr = MatE(m, index);
    for (int i=0 ; i<length ; i++)
        rowPtr[i] = row[i];
}

template<class MatT>
inline void
SetRow(MatT* dst, int dstIndex, MatT* src, int srcIndex)
{
    ILOG_VAR(Impala.Core.Matrix.SetRow);
    int length = MatNrCol(dst);
    if (length != MatNrCol(src))
        ILOG_ERROR_NODE("nonconformant operands");
    if (dstIndex < 0 || dstIndex >= MatNrRow(dst))
        ILOG_ERROR_NODE("faulty row dstIndex: " << dstIndex);
    if (srcIndex < 0 || srcIndex >= MatNrRow(src))
        ILOG_ERROR_NODE("faulty row srcIndex: " << srcIndex);
    typename MatT::StorType* dstRowPtr = MatE(dst, dstIndex);
    typename MatT::StorType* srcRowPtr = MatE(src, srcIndex);
    for (int i=0 ; i<length ; i++)
        dstRowPtr[i] = srcRowPtr[i];
}

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

#endif

---