PolynomialKernel.h

// PolynomialKernel.h:
// -------------------
//

#ifndef __IPSDKMATH_POLYNOMIALKERNEL_H__
#define __IPSDKMATH_POLYNOMIALKERNEL_H__

#include <IPSDKMath/IPSDKMathExports.h>
#include <IPSDKUtil/Tools/ForceInline.h>
#include <IPSDKMath/Clustering/ClusteringTypes.h>
#include <IPSDKMath/LinearAlgebra/LinearAlgebraUtils.h>

namespace ipsdk {
namespace math {
namespace clustering {


class IPSDKMATH_API PolynomialKernel
{
protected:
    PolynomialKernel(const ipReal64 alpha,
                     const ipReal64 beta,
                     const ipReal64 degree);

public:
    ~PolynomialKernel();

// methods
public:
    static PolynomialKernel     createKernel(const ipReal64 alpha,
                                             const ipReal64 beta,
                                             const ipReal64 degree);

    inline ipReal64             getAlpha() const
    {
        return _alpha;
    }

    inline ipReal64             getBeta() const
    {
        return _beta;
    }

    inline ipReal64             getDegree() const
    {
        return _degree;
    }

    IPSDK_FORCEINLINE ipReal64  value(const Vector& v) const
    {
        return value(v, v);
    }

    IPSDK_FORCEINLINE ipReal64  value(const Vector& v0,
                                      const Vector& v1) const
    {
        // compute dot product between elements
        const ipReal64 dotProd = boost::numeric::ublas::inner_prod(v0, v1);

        return std::pow(_alpha*dotProd + _beta, _degree);
    }

// attributes
protected:
    ipReal64        _alpha;

    ipReal64        _beta;

    ipReal64        _degree;
};


} // end of namespace clustering
} // end of namespace math
} // end of namespace ipsdk

#endif // __IPSDKMATH_POLYNOMIALKERNEL_H__