Gaussian.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 Alesis Novik
 * Written (W) 2014 Parijat Mazumdar
 * Copyright (C) 2011 Berlin Institute of Technology and Max-Planck-Society
 */

#ifndef _GAUSSIAN_H__
#define _GAUSSIAN_H__

#include <shogun/lib/config.h>

#ifdef HAVE_LAPACK

#include <shogun/distributions/Distribution.h>
#include <shogun/features/DotFeatures.h>
#include <shogun/lib/common.h>
#include <shogun/mathematics/Math.h>

#include <shogun/lib/SGVector.h>
#include <shogun/lib/SGMatrix.h>

namespace shogun
{
class CDotFeatures;

enum ECovType
{
    FULL,
    DIAG,
    SPHERICAL
};

class CGaussian : public CDistribution
{
    public:
        CGaussian();
        CGaussian(const SGVector<float64_t> mean, SGMatrix<float64_t> cov, ECovType cov_type=FULL);
        virtual ~CGaussian();

        void init();

        virtual bool train(CFeatures* data=NULL);

        virtual int32_t get_num_model_parameters();

        virtual float64_t get_log_model_parameter(int32_t num_param);

        virtual float64_t get_log_derivative(
            int32_t num_param, int32_t num_example);

        virtual float64_t get_log_likelihood_example(int32_t num_example);

        virtual float64_t update_params_em(float64_t* alpha_k, int32_t len);

        virtual float64_t compute_PDF(SGVector<float64_t> point)
        {
            return CMath::exp(compute_log_PDF(point));
        }

        virtual float64_t compute_log_PDF(SGVector<float64_t> point);

        virtual SGVector<float64_t> get_mean();

        virtual void set_mean(const SGVector<float64_t> mean);

        virtual SGMatrix<float64_t> get_cov();

        virtual void set_cov(SGMatrix<float64_t> cov);

        inline ECovType get_cov_type()
        {
            return m_cov_type;
        }

        inline void set_cov_type(ECovType cov_type)
        {
            m_cov_type = cov_type;
        }

        inline SGVector<float64_t> get_d()
        {
            return m_d;
        }

        void set_d(const SGVector<float64_t> d);

        inline SGMatrix<float64_t> get_u()
        {
            return m_u;
        }

        inline void set_u(SGMatrix<float64_t> u)
        {
            m_u = u;
        }

        SGVector<float64_t> sample();

        static CGaussian* obtain_from_generic(CDistribution* distribution);

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

    private:
        void register_params();

        void decompose_cov(SGMatrix<float64_t> cov);

    protected:
        float64_t m_constant;
        SGVector<float64_t> m_d;
        SGMatrix<float64_t> m_u;
        SGVector<float64_t> m_mean;
        ECovType m_cov_type;
};
}
#endif //HAVE_LAPACK
#endif //_GAUSSIAN_H__