Changeset 214 for bdm/stat/emix.h

Timestamp:

11/27/08 23:31:48 (15 years ago)

Author:

smidl

Message:

debug asserts for infinite likelihoods

Files:

• bdm/stat/emix.h (modified) (4 diffs)

bdm/stat/emix.h

@@ -61 +61 @@
         };
         double evallogcond ( const vec &val, const vec &cond ) {
+                double tmp;
                 vec nom_val ( rv.count() +rvc.count() );
                 dl.fill_val_cond ( nom_val,val,cond );
-                return exp ( nom->evallog ( nom_val ) - den->evallog ( cond ) );
+                tmp = exp ( nom->evallog ( nom_val ) - den->evallog ( cond ) );
+                it_assert_debug(std::isfinite(tmp),"Infinite value");
+                return tmp;
         }
         //! Object takes ownership of nom and will destroy it
@@ -106 +109 @@
                 double sum = 0.0;
                 for ( i = 0;i < w.length();i++ ) {sum += w ( i ) * exp ( Coms ( i )->evallog ( val ) );}
-                return log ( sum );
+                if (sum==0.0){sum=std::numeric_limits<double>::epsilon();}
+                double tmp=log ( sum );
+                it_assert_debug(std::isfinite(tmp),"Infinite");
+                return tmp;
         };
         vec evallog_m ( const mat &Val ) const {
@@ -115 +121 @@
                 return log ( x );
         };
+        //! Auxiliary function that returns pdflog for each component
         mat evallog_M ( const mat &Val ) const {
                 mat X ( w.length(), Val.cols() );
@@ -249 +256 @@
                         tmp+=epdfs ( i )->evallog ( dls ( i )->get_val ( val ) );
                 }
+                it_assert_debug(std::isfinite(tmp),"Infinite");
                 return tmp;
         }