Changeset 741

Show
Ignore:
Timestamp:
11/26/09 22:46:01 (14 years ago)
Author:
smidl
Message:

Stress tests are passing now. Missing validate calls are filled...

Location:
library
Files:
7 modified

Legend:

Unmodified
Added
Removed

  • library/bdm/base/bdmbase.h

    r738 r741  
    437437        //! Array<vec> version of evallogcond 
    438438        virtual vec evallogcond_mat ( const Array<vec> &Yt, const vec &cond ) { 
    439                 bdm_error ( "Not implemented" ); 
    440                 return vec(); 
     439                vec v ( Yt.length() ); 
     440                for ( int i = 0; i < Yt.length(); i++ ) { 
     441                        v ( i ) = evallogcond ( Yt( i ), cond ); 
     442                } 
     443                return v; 
    441444        } 
    442445 

  • library/bdm/estim/arx.cpp

    r738 r741  
    1212 
    1313        dyad.set_subvector ( 0, yt ); 
    14         dyad.set_subvector ( dimy, cond ); 
     14        if (cond.length()>0) 
     15                dyad.set_subvector ( dimy, cond ); 
    1516        // possible "1" is there from the beginning 
    1617 

  • library/bdm/estim/mixtures.h

    r737 r741  
    7777                        Coms ( i ) = ( BMEF* ) Coms0 ( i )->_copy_(); 
    7878                } 
     79                weights.set_parameters(alpha0); 
     80                weights.validate(); 
    7981                build_est(); 
    8082        } 

  • library/bdm/stat/exp_family.h

    r739 r741  
    518518        void set_parameters ( const vec &beta0 ) { 
    519519                est.set_parameters ( beta0 ); 
     520                est.validate(); 
    520521                if ( evalll ) { 
    521522                        last_lognc = est.lognc(); 

  • library/bdm/stat/merger.cpp

    r739 r741  
    204204        ARX A0; 
    205205        A0.set_statistics ( dim, V0 ); //initial guess of Mix: 
     206        A0.validate(); 
    206207 
    207208        Mix.init ( &A0, Smp_ex, Ncoms ); 

  • library/tests/stresssuite/merger_iter_stress.cpp

    r721 r741  
    2727 
    2828        f1->set_parameters ( "4 3", mat ( "0.4 0.3; 0.3 0.4" ) ); 
     29        f1->validate(); 
    2930        f2->set_parameters ( "1 3", mat ( "0.3 -0.2; -0.2 0.3" ) ); 
     31        f2->validate(); 
    3032        f3->set_parameters ( "2", mat ( "0.4" ) ); 
    31  
     33        f3->validate(); 
     34         
    3235        pdf_array A ( 3 ); 
    3336        A ( 0 ) = f1; 
    … …  
    4952        g0.set_rv ( xy ); 
    5053        g0.set_parameters ( vec ( "4 4" ), mat ( "1 0; 0 1" ) ); 
    51  
     54        g0.validate(); 
     55         
    5256        M.set_parameters ( 5 ); 
    5357        M.set_method ( LOGNORMAL, 1.2 ); 
    5458        M.set_support ( g0, 400 ); 
     59        M.validate(); 
    5560        M.merge(); 
     61         
    5662 
    5763        MixEF &MM = M._Mix(); 

  • library/tests/stresssuite/mixtures_stress.cpp

    r722 r741  
    8888        ARX M1; 
    8989        M1.set_statistics ( 2, V0 - Vp, 8 ); 
     90        M1.validate(); 
    9091        ARX M2; 
    9192        M2.set_statistics ( 2, V0 + Vp, 8 ); 
     93        M2.validate(); 
    9294 
    9395        // Build mixture model 
    … …  
    110112        ARX Aflat; 
    111113        Aflat.set_statistics ( 2, V0, 7 ); 
     114        Aflat.validate(); 
    112115        MixEF RND; 
    113116        RND.init ( &Aflat, Data, 10 ); // already initialized!