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

---

KullbackDivergence.h

Go to the documentation of this file.

#ifndef Impala_Core_Vector_KullbackDivergence_h
#define Impala_Core_Vector_KullbackDivergence_h

#include "Core/Vector/VectorTem.h"

namespace Impala
{
namespace Core
{
namespace Vector
{

template <class ElemT>
inline ElemT
KullbackDivergence(const VectorTem<ElemT>& v1, const VectorTem<ElemT>& v2)
{
    ElemT score = 0;
    for (int i=0 ; i<v1.Size() ; i++)
        if((v2[i]!=0)&&(v1[i]!=0))
            score += v1[i]*log(v1[i]/v2[i]);
    return score;
}

template <class ElemT>
class KullbackDivergenceFtor
{
public:

    ElemT
    DoIt(const VectorTem<ElemT>& v1, const VectorTem<ElemT>& v2)
    {
        ElemT score = 0;
        for (int i=0 ; i<v1.Size() ; i++)
            if((v2[i]!=0)&&(v1[i]!=0))
                score += v1[i]*log(v1[i]/v2[i]);
        return score;
    }
};


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

#endif

---