KMeans.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) 1999-2008 Gunnar Raetsch
 * Written (W) 2007-2009 Soeren Sonnenburg
 * Copyright (C) 1999-2009 Fraunhofer Institute FIRST and Max-Planck-Society
 */

#ifndef _KMEANS_H__
#define _KMEANS_H__

#include <stdio.h>
#include "lib/common.h"
#include "lib/io.h"
#include "features/SimpleFeatures.h"
#include "distance/Distance.h"
#include "classifier/DistanceMachine.h"

namespace shogun
{
class CDistanceMachine;

class CKMeans : public CDistanceMachine
{
    public:
        CKMeans();

        CKMeans(int32_t k, CDistance* d);
        virtual ~CKMeans();

        virtual inline EClassifierType get_classifier_type() { return CT_KMEANS; }

        virtual bool train(CFeatures* data=NULL);

        virtual bool load(FILE* srcfile);

        virtual bool save(FILE* dstfile);

        inline void set_k(int32_t p_k)
        {
            ASSERT(p_k>0);
            this->k=p_k;
        }

        inline int32_t get_k()
        {
            return k;
        }

        inline void set_max_iter(int32_t iter)
        {
            ASSERT(iter>0);
            max_iter=iter;
        }

        inline float64_t get_max_iter()
        {
            return max_iter;
        }

        inline void get_radi(float64_t*& radi, int32_t& num)
        {
            radi=R;
            num=k;
        }

        inline void get_centers(float64_t*& centers, int32_t& dim, int32_t& num)
        {
            centers=mus;
            dim=dimensions;
            num=k;
        }

        inline void get_radiuses(float64_t** radii, int32_t* num)
        {
            size_t sz=sizeof(*R)*k;
            *radii=(float64_t*) malloc(sz);
            ASSERT(*radii);

            memcpy(*radii, R, sz);
            *num=k;
        }

        inline void get_cluster_centers(
            float64_t** centers, int32_t* dim, int32_t* num)
        {
            size_t sz=sizeof(*mus)*dimensions*k;
            *centers=(float64_t*) malloc(sz);
            ASSERT(*centers);

            memcpy(*centers, mus, sz);
            *dim=dimensions;
            *num=k;
        }

        inline int32_t get_dimensions()
        {
            return dimensions;
        }

    protected:
        void clustknb(bool use_old_mus, float64_t *mus_start);

        virtual CLabels* classify()
        {
            SG_NOTIMPLEMENTED;
            return NULL;
        }

        virtual CLabels* classify(CFeatures* data)
        {
            SG_NOTIMPLEMENTED;
            return NULL;
        }



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

    protected:
        int32_t max_iter;

        int32_t k;

        int32_t dimensions;

        float64_t* R;
        
        float64_t* mus;

    private:
        float64_t* Weights;
};
}
#endif