source: trunk/c++_tools/classifier/Kernel_SEV.cc @ 658

// $Id$

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

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

namespace theplu {
namespace classifier {  


  Kernel_SEV::Kernel_SEV(const MatrixLookup& data, const KernelFunction& kf,
                         const bool own)
    : Kernel(data,kf, own)
  {
    build_kernel();
  }


  Kernel_SEV::Kernel_SEV(const MatrixLookupWeighted& data, 
                         const KernelFunction& kf, const bool own)
    : Kernel(data,kf, own)
  {
    kernel_matrix_ = utility::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_)(DataLookupWeighted1D(data,i,false),
                 DataLookupWeighted1D(data,j,false));
  }


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


  void Kernel_SEV::build_kernel(void)
  {
    kernel_matrix_ = utility::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::operator()(const size_t row,const size_t column) const
  { 
    return kernel_matrix_(row,column); 
  }
  

  const Kernel_SEV* Kernel_SEV::make_kernel(const MatrixLookup& data, 
                                            const bool own) const
  {
    return new Kernel_SEV(data, *kf_, own);
  }


  const Kernel_SEV* Kernel_SEV::make_kernel(const MatrixLookupWeighted& data, 
                            const bool own) const
  {
    return new Kernel_SEV(data, *kf_, own);
  }


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


}} // of namespace classifier and namespace theplu