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

---

WeibullSim.h

Go to the documentation of this file.

#ifndef Impala_Core_Vector_WeibullSim_h
#define Impala_Core_Vector_WeibullSim_h

#include "Core/Vector/VectorTem.h"
#include "Basis/NativeTypes.h"

namespace Impala
{
namespace Core
{
namespace Vector
{


inline Real64
WeibullSim(Real64 b1, Real64 g1, Real64 b2, Real64 g2)
{
    if (g1==0 || g2==0 || b1==0 || b2==0)
        return 0;
    Real64 g = (g1>g2 ? g2/g1 : g1/g2);
    Real64 b = (b1>b2 ? b2/b1 : b1/b2);
    return g*b;
}

template <class ElemT>
inline Real64
WeibullSim(const VectorTem<ElemT>& v1, const VectorTem<ElemT>& v2)
{
    Real64 score = 0;
    for (int i=0 ; i<v1.Size() ; i+=2)
    {
        //score += WeibullSim((Real64) v1[i], (Real64) v1[i+1],
        //                    (Real64) v2[i], (Real64) v2[i+1]);
        // This seems to be about 10% faster on windows...
        Real64 b1 = v1[i];
        Real64 g1 = v1[i+1];
        Real64 b2 = v2[i];
        Real64 g2 = v2[i+1];
        if (g1!=0 && g2!=0 && b1!=0 && b2!=0)
        {
            Real64 g = (g1>g2 ? g2/g1 : g1/g2);
            Real64 b = (b1>b2 ? b2/b1 : b1/b2);
            score += g*b;
        }
    }
    return score / (v1.Size() / 2);
}

template <class ElemT>
class WeibullSimFtor
{
public:

    Real64
    DoIt(const VectorTem<ElemT>& v1, const VectorTem<ElemT>& v2)
    {
        Real64 score = 0;
        for (int i=0 ; i<v1.Size() ; i+=2)
        {
            //score += WeibullSim((Real64) v1[i], (Real64) v1[i+1],
            //                    (Real64) v2[i], (Real64) v2[i+1]);
            // This seems to be about 10% faster on windows...
            Real64 b1 = v1[i];
            Real64 g1 = v1[i+1];
            Real64 b2 = v2[i];
            Real64 g2 = v2[i+1];
            if (g1!=0 && g2!=0 && b1!=0 && b2!=0)
            {
                Real64 g = (g1>g2 ? g2/g1 : g1/g2);
                Real64 b = (b1>b2 ? b2/b1 : b1/b2);
                score += g*b;
            }
        }
        return score;
        //return score / (v1.Size() / 2);
    }
};

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

#endif

---