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

---

Sum.h

Go to the documentation of this file.

#ifndef Impala_Core_Vector_Sum_h
#define Impala_Core_Vector_Sum_h

#include "Core/Vector/VectorTem.h"
#include "Core/Vector/VectorSet.h"
#include "Core/Vector/AddAssign.h"

namespace Impala
{
namespace Core
{
namespace Vector
{


template <class ElemT>
inline ElemT
Sum(VectorTem<ElemT> v)
{
    ElemT* p = v.GetData();
    ElemT res = 0;
    for (int i=0 ; i<v.Size() ; i++)
        res += p[i];
    return res;
}


template <class ArrayT>
inline void
Sum(typename VectorSet<ArrayT>::VectorT* dst, VectorSet<ArrayT>* src,
    int startIdx, int stepSize, int nrSteps)
{
    typedef typename VectorSet<ArrayT>::VectorT VectorT;

    if (!src->HasConstVectorSize())
    {
        std::cout << "Vector::Sum: works for constant vector size only" << std::endl;
        return;
    }
    VectorT res(src->GetVectorLength(0));
    res = 0;
    int endIdx = startIdx + nrSteps * stepSize;
    for (int i=startIdx ; i<endIdx ; i+=stepSize)
        AddAssign(res, src->GetVector(i, true));
    if (!dst)
        dst = new VectorT(res);
    else
        *dst = res;
}

template <class ArrayT>
inline void
Sum(typename VectorSet<ArrayT>::VectorT* dst, VectorSet<ArrayT>* src, int size,
    bool* filter)
{
    typedef typename VectorSet<ArrayT>::VectorT VectorT;

    if (!src->HasConstVectorSize())
    {
        std::cout << "Vector::Sum: works for constant vector size only" << std::endl;
        return;
    }
    VectorT res(src->GetVectorLength(0));
    res = 0;
    for (int i=0 ; i<size ; i++)
    {
        if (!filter || filter[i])
        {
            AddAssign(res, src->GetVector(i, true));
        }
    }
    if (!dst)
        dst = new VectorT(res);
    else
        *dst = res;
}




} // namespace Vector
} // namespace Core
} // namespace Impala

#endif

---