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

// $Id: KernelWeighted_SEV.cc 545 2006-03-06 13:35:45Z peter $

#include <c++_tools/classifier/KernelWeighted_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 {  


  KernelWeighted_SEV::KernelWeighted_SEV(const MatrixLookup& data, 
                                         const KernelFunction& kf, 
                                         const MatrixLookup& weights)
    : Kernel(data,kf,weights)
  {
    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), 
                 DataLookup1D(*weights_,i,false),
                 DataLookup1D(*weights_,j,false));
  }
  

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

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



}} // of namespace classifier and namespace theplu