source: trunk/lib/classifier/Kernel.h @ 542

// $Id: Kernel.h 542 2006-03-05 15:46:47Z peter $

#ifndef _theplu_classifier_kernel_
#define _theplu_classifier_kernel_

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

#include <cctype>

namespace theplu {
namespace classifier {

  ///
  ///   @brief Base Class for Kernels. 
  ///
  ///   Class taking care of the \f$NxN\f$ kernel matrix, where
  ///   \f$N\f$ is number of samples. Each element in the Kernel
  ///   matrix is the scalar product of the corresponding pair of
  ///   samples. Type of Kernel is defined by a KernelFunction.
  ///   
  ///   @note If the KernelFunction is destroyed, the Kernel is no
  ///   longer defined.
  ///
  class Kernel
  {
    
  public:

    ///
    ///   Constructor taking the data matrix and KernelFunction as
    ///   input.Each column in the data matrix corresponds to one
    ///   sample. 
    ///
    ///   @note Can not handle NaNs. 
    ///
    Kernel(const MatrixLookup& data, const KernelFunction& kf) 
      : data_(data), kf_(&kf) {};
    
    ///
    ///   Destructor
    ///
    virtual ~Kernel(void) {};

    ///
    /// @return element at position (\a row, \a column) of the Kernel
    /// matrix
    ///
    virtual double operator()(const size_t row, const size_t column) const=0;

    ///
    /// @return number columns in Kernel
    ///
    inline size_t columns(void) const { return size(); } 

    ///
    /// @return number of rows in Kernel
    ///
    inline size_t rows(void) const { return size(); } 

    ///
    /// @brief number of samples
    ///
    inline size_t size(void) const { return data_.columns(); } 

    
    virtual double element(const DataLookup1D& vec, const size_t i) const=0;
    virtual double element(const DataLookup1D& vec, const DataLookup1D& w,
                           const size_t i) const=0;

  protected:
    // Peter should be a copy
    const MatrixLookup& data_;
    const KernelFunction* kf_;

  private:
    ///
    /// Copy constructor (not implemented)
    ///
    Kernel(const Kernel&);


  }; // class Kernel

}} // of namespace classifier and namespace theplu

#endif