Changeset 211 for bdm/stat/libEF.h

Show
Ignore:
Timestamp:
11/13/08 20:00:53 (16 years ago)
Author:
smidl
Message:

prejmenovani evalpdflog a evalcond

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • bdm/stat/libEF.h

    r205 r211  
    4646        virtual void dupdate ( mat &v ) {it_error ( "Not implemented" );}; 
    4747        //!Evaluate normalized log-probability 
    48         virtual double evalpdflog_nn ( const vec &val ) const{it_error ( "Not implemented" );return 0.0;}; 
     48        virtual double evallog_nn ( const vec &val ) const{it_error ( "Not implemented" );return 0.0;}; 
    4949        //!Evaluate normalized log-probability 
    50         virtual double evalpdflog ( const vec &val ) const {double tmp;tmp= evalpdflog_nn ( val )-lognc();it_assert_debug(std::isfinite(tmp),"why?"); return tmp;} 
     50        virtual double evallog ( const vec &val ) const {double tmp;tmp= evallog_nn ( val )-lognc();it_assert_debug(std::isfinite(tmp),"why?"); return tmp;} 
    5151        //!Evaluate normalized log-probability for many samples 
    52         virtual vec evalpdflog ( const mat &Val ) const { 
     52        virtual vec evallog ( const mat &Val ) const { 
    5353                vec x ( Val.cols() ); 
    54                 for ( int i=0;i<Val.cols();i++ ) {x ( i ) =evalpdflog_nn ( Val.get_col ( i ) ) ;} 
     54                for ( int i=0;i<Val.cols();i++ ) {x ( i ) =evallog_nn ( Val.get_col ( i ) ) ;} 
    5555                return x-lognc(); 
    5656        } 
     
    129129        mat sample ( int N ) const; 
    130130        double eval ( const vec &val ) const ; 
    131         double evalpdflog_nn ( const vec &val ) const; 
     131        double evallog_nn ( const vec &val ) const; 
    132132        double lognc () const; 
    133133        vec mean() const {return mu;} 
     
    193193        void mean_mat ( mat &M, mat&R ) const; 
    194194        //! In this instance, val= [theta, r]. For multivariate instances, it is stored columnwise val = [theta_1 theta_2 ... r_1 r_2 ] 
    195         double evalpdflog_nn ( const vec &val ) const; 
     195        double evallog_nn ( const vec &val ) const; 
    196196        double lognc () const; 
    197197 
     
    227227        vec mean() const {return beta/sum ( beta );}; 
    228228        //! In this instance, val is ... 
    229         double evalpdflog_nn ( const vec &val ) const {return ( beta-1 ) *log ( val );}; 
     229        double evallog_nn ( const vec &val ) const {return ( beta-1 ) *log ( val );}; 
    230230        double lognc () const { 
    231231                double gam=sum ( beta ); 
     
    318318        //! TODO: is it used anywhere? 
    319319//      mat sample ( int N ) const; 
    320         double evalpdflog ( const vec &val ) const; 
     320        double evallog ( const vec &val ) const; 
    321321        double lognc () const; 
    322322        //! Returns poiter to alpha and beta. Potentially dengerous: use with care! 
     
    358358        euni ( const RV rv ) :epdf ( rv ) {} 
    359359        double eval ( const vec &val ) const  {return nk;} 
    360         double evalpdflog ( const vec &val ) const  {return lnk;} 
     360        double evallog ( const vec &val ) const  {return lnk;} 
    361361        vec sample() const { 
    362362                vec smp ( rv.count() ); 
     
    538538        vec sample() const {it_error ( "Not implemented" );return 0;} 
    539539        //! inherited operation : NOT implemneted 
    540         double evalpdflog ( const vec &val ) const {it_error ( "Not implemented" );return 0.0;} 
     540        double evallog ( const vec &val ) const {it_error ( "Not implemented" );return 0.0;} 
    541541        vec mean() const { 
    542542                vec pom=zeros ( rv.count() ); 
     
    599599double enorm<sq_T>::eval ( const vec &val ) const { 
    600600        double pdfl,e; 
    601         pdfl = evalpdflog ( val ); 
     601        pdfl = evallog ( val ); 
    602602        e = exp ( pdfl ); 
    603603        return e; 
     
    605605 
    606606template<class sq_T> 
    607 double enorm<sq_T>::evalpdflog_nn ( const vec &val ) const { 
     607double enorm<sq_T>::evallog_nn ( const vec &val ) const { 
    608608        // 1.83787706640935 = log(2pi) 
    609609        double tmp=-0.5* ( R.invqform ( mu-val ) );// - lognc();