Changeset 404 for library/bdm/stat/exp_family.cpp

Timestamp:

07/02/09 22:16:05 (15 years ago)

Author:

smidl

Message:

Change in epdf: evallog returns -inf for points out of support. Merger is aware of it now.

Files:

• library/bdm/stat/exp_family.cpp (modified) (3 diffs)

library/bdm/stat/exp_family.cpp

@@ -24 +24 @@
         if ( dimx==1 ) { //same as the following, just quicker.
                 double r = val ( vend ); //last entry!
-                if (r<0) return -1e+100;
+                if (r<0) return -inf;
                 vec Psi ( nPsi+dimx );
                 Psi ( 0 ) = -1.0;
@@ -36 +36 @@
                 fsqmat R ( reshape ( val ( nPsi*dimx,vend ),dimx,dimx ) );
                 double ldetR=R.logdet();
-                if (ldetR) return -1e+100;
+                if (ldetR) return -inf;
                 mat Tmp=concat_vertical ( -eye ( dimx ),Th );
                 fsqmat iR ( dimx );
@@ -190 +190 @@
         int i;
 
-        if (any(val<=0.)) return -1e100;
-        if (any(beta<=0.)) return -1e100;
+        if (any(val<=0.)) return -inf;
+        if (any(beta<=0.)) return -inf;
         for ( i=0; i<dim; i++ ) {
                 res += ( alpha ( i ) - 1 ) *std::log ( val ( i ) ) - beta ( i ) *val ( i );