GEMPLP.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) 2014 Jiaolong Xu
 * Copyright (C) 2014 Jiaolong Xu
 */

#ifndef __GEMPLP_H__
#define __GEMPLP_H__

#include <shogun/lib/config.h>
#include <shogun/lib/SGVector.h>
#include <shogun/structure/FactorGraph.h>
#include <shogun/structure/Factor.h>
#include <shogun/structure/MAPInference.h>
#include <shogun/lib/SGNDArray.h>

#include <vector>
#include <set>

namespace shogun
{
class CGEMPLP: public CMAPInferImpl
{
public:
    struct Parameter
    {
#ifndef DOXYGEN_SHOULD_SKIP_THIS
        Parameter(const int32_t max_iter = 1000,
                  const float64_t obj_del_thr = 0.0002,
                  const float64_t int_gap_thr = 0.0002)
            : m_max_iter(max_iter),
              m_obj_del_thr(obj_del_thr),
              m_int_gap_thr(int_gap_thr)
        {}
#endif

        int32_t m_max_iter;
        float64_t m_obj_del_thr;
        float64_t m_int_gap_thr;
    };

public:
    CGEMPLP();

    CGEMPLP(CFactorGraph* fg, Parameter param = Parameter());

    virtual ~CGEMPLP();

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

    virtual float64_t inference(SGVector<int32_t> assignment);

private:
    void init();

    void update_messages(int32_t id_region);

public:
    void max_in_subdimension(SGNDArray<float64_t> tar_arr, SGVector<int32_t> &subset_inds, SGNDArray<float64_t> &max_res) const;

    int32_t find_intersection_index(SGVector<int32_t> region_A, SGVector<int32_t> region_B);

    SGNDArray<float64_t> convert_energy_to_potential(CFactor* factor);

public:
    Parameter m_param;
    CDynamicObjectArray* m_factors;
    std::vector<SGVector<int32_t> > m_all_intersections;
    std::vector<std::set<int32_t> > m_region_intersections;
    std::vector<std::vector<SGVector<int32_t> > > m_region_inds_intersections;
    std::vector<SGNDArray<float64_t> > m_msgs_into_intersections;
    std::vector<std::vector<SGNDArray<float64_t> > > m_msgs_from_region;
    std::vector<SGNDArray<float64_t> > m_theta_region;
};
}
#endif