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

---

HistogramIntersection.h

Go to the documentation of this file.

#ifndef Impala_Core_Vector_HistogramIntersection_h
#define Impala_Core_Vector_HistogramIntersection_h

#include "Core/Vector/VectorTem.h"

namespace Impala
{
namespace Core
{
namespace Vector
{


template <class ElemT>
inline ElemT
HistogramIntersection(const VectorTem<ElemT>& v1, const VectorTem<ElemT>& v2)
{
    ElemT score = 0;
    for (int i=0 ; i<v1.Size() ; i++)
        score += Impala::Min(v1[i], v2[i]);
    return score;
}

inline double
HistogramIntersectionDouble(const Vector::VectorTem<double>& v1,
                            const Vector::VectorTem<double>& v2)
{
    double accumulator = 0;
    for(int i = 0; i < v1.Size(); i++)
    {
        accumulator += Impala::Min(v1[i], v2[i]);
    }
    return accumulator;
}

template <class ElemT>
class HistogramIntersectionFtor
{
public:

    ElemT
    DoIt(const VectorTem<ElemT>& v1, const VectorTem<ElemT>& v2)
    {
        ElemT score = 0;
        for (int i=0 ; i<v1.Size() ; i++)
            score += Impala::Min(v1[i], v2[i]);
        return score;
    }
};

template <class FType>
inline void
DoSlabHistogramIntersection(FType* C, const FType* A, const FType* B,
                            unsigned int slabSizeA, unsigned int slabSizeB,
                            unsigned int vectorSize)
{
    #pragma omp parallel for
    for (unsigned int a = 0; a < slabSizeA; a++)
    {
        for (unsigned int b = 0; b < slabSizeB; b++)
        {
            const FType* baseA = A + (a * vectorSize);
            const FType* baseB = B + (b * vectorSize);
            FType sum = 0;
            for (unsigned int k = 0; k < vectorSize; k++)
            {
                if (baseA[k] < baseB[k])
                {
                    sum += baseA[k];
                }
                else
                {
                    sum += baseB[k];
                }
            }
            C[b * slabSizeA + a] = (FType)sum;
        }
    }
}

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

#endif

---