MCBayesNetGenerator.h

Go to the documentation of this file.

#ifndef GUM_MC_BAYES_NET_GENERATOR
#define GUM_MC_BAYES_NET_GENERATOR

#include <agrum/agrum.h>

#include <fstream>
#include <iostream>
#include <set>
#include <vector>

#include <sstream>

#ifdef HAVE_DIRENT_H
#  include <dirent.h>
#else
#  include <agrum/core/mvsc/dirent.h>
#endif

#include <agrum/BN/BayesNet.h>
#include <agrum/BN/generator/IBayesNetGenerator.h>
#include <agrum/BN/generator/simpleCPTDisturber.h>
#include <agrum/BN/generator/simpleCPTGenerator.h>
#include <agrum/BN/inference/lazyPropagation.h>
#include <agrum/core/hashTable.h>
#include <agrum/multidim/potential.h>
#include <agrum/variables/labelizedVariable.h>

namespace gum {
  template < typename GUM_SCALAR,
             template < typename > class ICPTGenerator = SimpleCPTGenerator,
             template < typename > class ICPTDisturber = SimpleCPTDisturber >
  class MCBayesNetGenerator
      : public IBayesNetGenerator< GUM_SCALAR, ICPTGenerator >
      , public ICPTDisturber< GUM_SCALAR > {
    public:
    // ############################################################################
    // ############################################################################

    MCBayesNetGenerator(Size nbrNodes,
                        Size maxArcs,
                        Idx  maxModality = 2,
                        Size iteration = 5000,
                        Idx  p = 30,
                        Idx  q = 40);

    explicit MCBayesNetGenerator(BayesNet< GUM_SCALAR > bayesNet,
                                 Size                   iteration = 5000,
                                 Idx                    p = 30,
                                 Idx                    q = 40);

    ~MCBayesNetGenerator() override;


    // ############################################################################
    // ############################################################################

    void generateBN(BayesNet< GUM_SCALAR >& bayesNet) override;

    void disturbBN(BayesNet< GUM_SCALAR >& bayesNetinit, Size iteration = 0);


    // ############################################################################
    // ############################################################################

    Size iteration() const;
    Idx p() const;
    Idx q() const;


    void setIteration(Size iteration);
    void setP(Idx p);
    void setQ(Idx q);

    protected:
    Size                                          _iteration;
    Idx                                           _p, _q;
    bool                                          _disturbing;
    BayesNet< GUM_SCALAR >                        _bayesNettemp;
    HashTable< NodeId, Potential< GUM_SCALAR >* > _hashMarginal;

    bool __isPolytree();
    bool __connect(NodeId i, NodeId j);
    bool __directedPath(NodeId tail, NodeId head);
    void __insertArc(NodeId i, NodeId j);
    void __eraseArc(NodeId i, NodeId j, bool mustbeconnex = true);

    void __PMMx_poly();
    void __PMMx_multi();
    void __jump_poly();

    void __jump_multi();

    void __AorR();
    void __AR();
    virtual bool __checkConditions();

    // NOT USED ? void __createDAG( Size BNSize, Size iniRoot );

    // NOT USED ? std::vector<Idx>* __createPartDAG( Size BNSize, Size iniRoot
    // );

    bool __connect(NodeId i, NodeId j, NodeSet& excluded);

    bool __directedPath(NodeId tail, NodeId head, NodeSet& excluded);

    void __chooseNodes(NodeId& i, NodeId& j);

    void __chooseCloseNodes(NodeId& i, NodeId& j);

    void __transformPoly(Idx nbiter);

    void __createTree(Size BNSize);

    NodeId __createPartTree(Size BNSize, Idx& n);
  };


#ifndef GUM_NO_EXTERN_TEMPLATE_CLASS
  extern template class MCBayesNetGenerator< double >;
#endif

} /*namespace gum*/

#include <agrum/BN/generator/MCBayesNetGenerator_tpl.h>
#endif   // MCBAYESNETGENERATOR