source: trunk/lib/classifier/Kernel_SEV.cc @ 542

// $Id: Kernel_SEV.cc 542 2006-03-05 15:46:47Z peter $

#include <c++_tools/classifier/Kernel_SEV.h>

#include <c++_tools/classifier/DataLookup1D.h>
#include <c++_tools/classifier/Kernel.h>
#include <c++_tools/classifier/KernelFunction.h>
#include <c++_tools/classifier/MatrixLookup.h>
#include <c++_tools/gslapi/matrix.h>
#include <c++_tools/gslapi/vector.h>

namespace theplu {
namespace classifier {  


  Kernel_SEV::Kernel_SEV(const MatrixLookup& data, const KernelFunction& kf)
    : Kernel(data,kf)
  {
    kernel_matrix_ = gslapi::matrix(data.columns(),data.columns());
    for (size_t i=0; i<kernel_matrix_.rows(); i++) 
      for (size_t j=i; j<kernel_matrix_.columns(); j++)
        kernel_matrix_(i,j) = kernel_matrix_(j,i) =
          (*kf_)(DataLookup1D(data_,i,false),DataLookup1D(data_,j,false));
  }


  double Kernel_SEV::element(const DataLookup1D& vec, const size_t i) const
  { 
    return kf_->operator()(vec, DataLookup1D(data_,i)); 
  }

  double Kernel_SEV::element(const DataLookup1D& vec, const DataLookup1D& w,
                             const size_t i) const
  { 
    return (*kf_)(vec, DataLookup1D(data_,i), w, DataLookup1D(w.size(),1.0));
  }


}} // of namespace classifier and namespace theplu