Merged regiongraph.* and daialg.* from SVN head,
[libdai.git] / include / dai / lc.h
index 57037d9..86034f7 100644 (file)
@@ -28,6 +28,7 @@
 #include <dai/enum.h>
 #include <dai/factorgraph.h>
 #include <dai/properties.h>
+#include <dai/exceptions.h>
 
 
 namespace dai {
@@ -48,9 +49,9 @@ class LC : public DAIAlgFG {
             size_t verbose;
             size_t maxiter;
             double tol;
-            ENUM6(CavityType,FULL,PAIR,PAIR2,PAIRINT,PAIRCUM,UNIFORM)
+            DAI_ENUM(CavityType,FULL,PAIR,PAIR2,PAIRINT,PAIRCUM,UNIFORM)
             CavityType cavity;
-            ENUM3(UpdateType,SEQFIX,SEQRND,NONE)
+            DAI_ENUM(UpdateType,SEQFIX,SEQRND,NONE)
             UpdateType updates;
             std::string cavainame;
             PropertySet cavaiopts;
@@ -65,6 +66,8 @@ class LC : public DAIAlgFG {
         LC(const LC & x) : DAIAlgFG(x), _pancakes(x._pancakes), _cavitydists(x._cavitydists), _phis(x._phis), _beliefs(x._beliefs), props(x.props), maxdiff(x.maxdiff) {}
         /// Clone function
         LC* clone() const { return new LC(*this); }
+        /// Create (virtual constructor)
+        virtual LC* create() const { return new LC(); }
         /// Construct LC object from a FactorGraph and parameters
         LC( const FactorGraph & fg, const PropertySet &opts );
         /// Assignment operator
@@ -87,26 +90,44 @@ class LC : public DAIAlgFG {
         long SetCavityDists( std::vector<Factor> &Q );
 
         void init();
+        /// Clear messages and beliefs corresponding to the nodes in ns
+        virtual void init( const VarSet &/*ns*/ ) {
+            DAI_THROW(NOT_IMPLEMENTED);
+        }
         Factor NewPancake (size_t i, size_t _I, bool & hasNaNs);
         double run();
 
         std::string identify() const;
         Factor belief (const Var &n) const { return( _beliefs[findVar(n)] ); }
-        Factor belief (const VarSet &/*ns*/) const { assert( 0 == 1 ); }
+        Factor belief (const VarSet &/*ns*/) const { 
+            DAI_THROW(NOT_IMPLEMENTED);
+            return Factor(); 
+        }
         std::vector<Factor> beliefs() const { return _beliefs; }
-        Complex logZ() const { return NAN; }
+        Real logZ() const { 
+            DAI_THROW(NOT_IMPLEMENTED);
+            return 0.0; 
+        }
         void CalcBelief (size_t i);
         const Factor &belief (size_t i) const { return _beliefs[i]; };
         const Factor &pancake (size_t i) const { return _pancakes[i]; };
         const Factor &cavitydist (size_t i) const { return _cavitydists[i]; };
 
-        void clamp( const Var &/*n*/, size_t /*i*/ ) { assert( 0 == 1 ); }
-        void undoProbs( const VarSet &/*ns*/ ) { assert( 0 == 1 ); }
-        void saveProbs( const VarSet &/*ns*/ ) { assert( 0 == 1 ); }
-        virtual void makeCavity(const Var & /*n*/) { assert( 0 == 1 ); }
-
+        void clamp( const Var &/*n*/, size_t /*i*/ ) { 
+            DAI_THROW(NOT_IMPLEMENTED);
+        }
+        void restoreFactors( const VarSet &/*ns*/ ) { 
+            DAI_THROW(NOT_IMPLEMENTED);
+        }
+        void backupFactors( const VarSet &/*ns*/ ) { 
+            DAI_THROW(NOT_IMPLEMENTED);
+        }
+        virtual void makeCavity(const Var & /*n*/) { 
+            DAI_THROW(NOT_IMPLEMENTED);
+        }
         void setProperties( const PropertySet &opts );
         PropertySet getProperties() const;
+        std::string printProperties() const;
         double maxDiff() const { return maxdiff; }
 };