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

// $Id: Kernel_SEV.cc 545 2006-03-06 13:35:45Z 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));
  }

  Kernel_SEV::Kernel_SEV(const Kernel_SEV& other, 
                         const std::vector<size_t>& index)
    : Kernel(other, index)
  {
  }

  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));
  }

  const Kernel* Kernel_SEV::selected(const std::vector<size_t>& index) const
  {
    return new Kernel_SEV(*this, index);
  }


}} // of namespace classifier and namespace theplu