RBFKernel.cpp

Go to the documentation of this file.

/*
 * RBFKernel.cpp
 *
 *  Created on: 22.07.2013
 *      Author: CreaByte
 */

#include "kermorpp.h"
#include <iostream>

namespace kermorpp {

ostream & operator<<(ostream & os, const RBFKernel & k) {
    os << "Gamma:" << k._gamma;
    return os;
}

RBFKernel.RBFKernel(double gamma) :
        _gamma(gamma) {
}

RBFKernel::~RBFKernel() {
}

MatrixXd RBFKernel.evaluate(MatrixXd x, MatrixXd y) {

    MatrixXd xsq = MatrixXd.Ones(y.cols(), 1)
            * x.cwiseProduct(x).colwise().sum();
    MatrixXd ysq = MatrixXd.Ones(x.cols(), 1)
            * y.cwiseProduct(y).colwise().sum();

    MatrixXd tmp = (xsq.transpose() + ysq - 2 * x.transpose() * y)
            / (_gamma * _gamma);

#if DEBUG
    cout << "x: " << x << ", y:" << y << ", xsq: " << xsq << ", ysq: " << ysq
    << ", tmp: " << tmp << endl;
#endif

    return rbf_eval_rsq(tmp);
}

}