ANOVAKernel.h

Go to the documentation of this file.

/*
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 3 of the License, or
 * (at your option) any later version.
 *
 * Written (W) 2011 Andrew Tereskin
 * Copyright (C) 2011 Berlin Institute of Technology and Max-Planck-Society
 */

#include <shogun/lib/config.h>

#ifndef ANOVAKERNEL_H_
#define ANOVAKERNEL_H_

#include <shogun/lib/common.h>
#include <shogun/kernel/DotKernel.h>
#include <shogun/features/Features.h>
#include <shogun/features/SimpleFeatures.h>

namespace shogun
{

class CDistance;

class CANOVAKernel: public CDotKernel
{
public:
    CANOVAKernel();

    CANOVAKernel(int32_t cache, int32_t d);

    CANOVAKernel(
        CSimpleFeatures<float64_t>* l, CSimpleFeatures<float64_t>* r, int32_t d, int32_t cache);

    virtual ~CANOVAKernel();

    virtual bool init(CFeatures* l, CFeatures* r);

    inline virtual EKernelType get_kernel_type() { return K_ANOVA; }

    inline virtual EFeatureType get_feature_type() { return F_DREAL; }

    inline virtual EFeatureClass get_feature_class() { return C_SIMPLE; }

    inline virtual const char* get_name() const { return "ANOVAKernel"; }

    inline int32_t get_cardinality() { return this->cardinality; }
    
    inline void set_cardinality(int32_t value) { this->cardinality = value; }

    float64_t compute_rec1(int32_t idx_a, int32_t idx_b);
    float64_t compute_rec2(int32_t idx_a, int32_t idx_b);
protected:

    virtual float64_t compute(int32_t idx_a, int32_t idx_b);

    void init();

    void allocate_arrays();

    virtual void cleanup();

    void register_params();

    virtual void load_serializable_post() throw (ShogunException);

private:
    float64_t compute_recursive1(float64_t* avec, float64_t* bvec, int32_t len);
    float64_t compute_recursive2(float64_t* avec, float64_t* bvec, int32_t len);

protected:
    int32_t cardinality;
    
    float64_t* DP;
    int32_t DP_len;
    
    float64_t* KD;
    float64_t* KS;
    float64_t* vec_pow;
};
}

#endif /* ANOVAKERNEL_H_ */