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

---

HistogramAccumulation.h

Go to the documentation of this file.

#ifndef Impala_Core_Vector_HistogramAccumulation_h
#define Impala_Core_Vector_HistogramAccumulation_h

#include "Core/Vector/VectorSet.h"
#include "Core/Vector/DivAssign.h"
#include "Core/Vector/Sum.h"
#include "Core/Vector/ElemMax.h"

namespace Impala
{
namespace Core
{
namespace Vector
{


template <class ArrayT>
inline void
HistogramAccumulation(VectorSet<ArrayT>* vs, bool doNormalizeHeight, 
                      bool doNormalizeWeight, bool* filter)
{
    for (int i=0 ; i<vs->Size() ; i++)
    {
        if (filter && (!filter[i]))
            continue;
        typename VectorSet<ArrayT>::VectorT v = vs->GetVector(i, true);
        typename VectorSet<ArrayT>::ElemT* p = v.GetData();
        for (int i=1 ; i<v.Size() ; i++)
            p[i] += p[i-1];
        if (doNormalizeHeight || doNormalizeWeight)
        {
            typename VectorSet<ArrayT>::ElemT factor;
            factor = (doNormalizeHeight) ? p[v.Size()-1] : Sum(v);
            if (factor == 0)
                v = 0;
            else
                DivAssign(v, factor);
        }
    }
}

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

#endif

---