#include "Kernel.hpp"

int main()
{
    Vector v1(10);

    v1[0]=1.0;
    v1[1]=1.0;
    v1 = v1+v1-v1;
    LinearKernel  lk;
    PolynomialKernel pk;
    RadialKernel<Vector >  rk(1.0f);

    std::cout << lk(v1,v1) << std::endl;
    std::cout << pk(v1,v1) << std::endl;
    std::cout << rk(v1,v1) << std::endl;

    SparseVector v2;
    v2[0]=1.0;
    v2[1]=1.0;
    v2 = v2+v2-v2;
    LinearKernel  lk2;
    PolynomialKernel pk2;
    RadialKernel<SparseVector>  rk2(1.0f);

    std::cout << lk2(v2,v2) << std::endl;
    std::cout << pk2(v2,v2) << std::endl;
    std::cout << rk2(v2,v2) << std::endl;
    return 1;
}