Changeset 713 for library/tests

Show
Ignore:
Timestamp:
11/10/09 11:03:09 (15 years ago)
Author:
mido
Message:

_m changed to _mat

emix.cfg prepared, but it is not yet debugged!

Location:
library/tests
Files:
1 added
10 modified

Legend:

Unmodified
Added
Removed
  • library/tests/CMakeLists.txt

    r700 r713  
    4141# using UnitTest++ 
    4242 
    43 SET(unit_test_configurations  egiw.cfg mlnorm.cfg edirich.cfg mprod.cfg generator.cfg epdfds.cfg pdfds.cfg test_user_info_input.cfg merger.cfg erroneous.cfg egamma.cfg test_user_info_external.cfg merger_error.cfg mgamma.cfg enorm.cfg) 
     43SET(unit_test_configurations  egiw.cfg mlnorm.cfg edirich.cfg mprod.cfg generator.cfg epdfds.cfg pdfds.cfg test_user_info_input.cfg merger.cfg erroneous.cfg egamma.cfg test_user_info_external.cfg merger_error.cfg mgamma.cfg enorm.cfg mmix.cfg emix.cfg) 
    4444 
    4545add_executable(testsuite arx_straux_test.cpp datalink_test.cpp  datasource_test.cpp egiw_test.cpp emix_test.cpp epdf_test.cpp logger_test.cpp LQG_test.cpp merger_test.cpp  
  • library/tests/arx_elem_test.cpp

    r682 r713  
    3838 
    3939        mlstudent* Ap = Ar.predictor_student(); 
    40         vec Ap_x = Ap->evallogcond_m ( X, empty_vec ); 
    41         vec ll_x = Ar.logpred_m ( X2 ); 
     40        vec Ap_x = Ap->evallogcond_mat ( X, empty_vec ); 
     41        vec ll_x = Ar.logpred_mat ( X2 ); 
    4242 
    4343        cout << "normalize : " << xstep*sum ( exp ( Ap_x ) ) << endl; 
  • library/tests/emix_test.cpp

    r693 r713  
    6868 
    6969        int N = 3; 
    70         mat Smp = M2.sample_m ( N ); 
     70        mat Smp = M2.sample_mat ( N ); 
    7171 
    7272        vec exp_ll ( "-5.0 -2.53563 -2.62171" ); 
    73         vec ll = M2.evallog_m ( Smp ); 
     73        vec ll = M2.evallog_mat ( Smp ); 
    7474        CHECK_CLOSE ( exp_ll, ll, 5.0 ); 
    7575 
     
    115115        Array<vec> actual(CurrentContext::max_trial_count); 
    116116        do { 
    117                 mat smp = distrib_obj.sample_m ( nsamples ); 
     117                mat smp = distrib_obj.sample_mat ( nsamples ); 
    118118                vec emu = sum ( smp, 2 ) / nsamples; 
    119119                actual( tc ) = emu; 
     
    135135        Array<mat> actual(CurrentContext::max_trial_count); 
    136136        do { 
    137                 mat smp = distrib_obj.sample_m ( nsamples ); 
     137                mat smp = distrib_obj.sample_mat ( nsamples ); 
    138138                vec emu = sum ( smp, 2 ) / nsamples; 
    139139                mat er = ( smp * smp.T() ) / nsamples - outer_product ( emu, emu ); 
  • library/tests/epdf_harness.cpp

    r706 r713  
    9595                } 
    9696 
    97                 vec lpzv_act = hepdf->evallog_m ( zero1n ); 
     97                vec lpzv_act = hepdf->evallog_mat ( zero1n ); 
    9898                CHECK_CLOSE_EX ( lpzv, lpzv_act, tolerance ); 
    9999 
     
    105105                } 
    106106 
    107                 lpzv_act = hepdf->evallog_m ( zeroa ); 
     107                lpzv_act = hepdf->evallog_mat ( zeroa ); 
    108108                CHECK_CLOSE_EX ( lpzv, lpzv_act, tolerance ); 
    109109        } 
     
    162162        Array<vec> actual(CurrentContext::max_trial_count); 
    163163        do { 
    164                 mat smp = hepdf->sample_m ( nsamples ); 
     164                mat smp = hepdf->sample_mat ( nsamples ); 
    165165                vec emu = smp * ones ( nsamples ) / nsamples; 
    166166                actual( tc ) = emu; 
     
    183183        Array<mat> actual(CurrentContext::max_trial_count); 
    184184        do { 
    185                 mat smp = hepdf->sample_m ( nsamples ); 
     185                mat smp = hepdf->sample_mat ( nsamples ); 
    186186                vec emu = smp * ones ( nsamples ) / nsamples; 
    187187                mat er = ( smp * smp.T() ) / nsamples - outer_product ( emu, emu ); 
     
    208208        Array<vec> actual(CurrentContext::max_trial_count); 
    209209        do { 
    210                 mat smp = mep.samplecond_m ( vec ( 0 ), nsamples ); 
     210                mat smp = mep.samplecond_mat ( vec ( 0 ), nsamples ); 
    211211                vec emu = sum ( smp, 2 ) / nsamples; 
    212212                actual( tc ) = emu; 
     
    229229        Array<mat> actual(CurrentContext::max_trial_count); 
    230230        do { 
    231                 mat smp = mep.samplecond_m ( vec ( 0 ), nsamples ); 
     231                mat smp = mep.samplecond_mat ( vec ( 0 ), nsamples ); 
    232232                vec emu = sum ( smp, 2 ) / nsamples; 
    233233                mat er = ( smp * smp.T() ) / nsamples - outer_product ( emu, emu ); 
  • library/tests/merger_2d_test.cpp

    r693 r713  
    5858        emix* MP = MM.epredictor(); //xy 
    5959 
    60         vec Res1 = M.evallog_m ( Grid ); 
    61         mat Res2 = ( MP )->evallog_M ( Grid ); 
     60        vec Res1 = M.evallog_mat ( Grid ); 
     61        mat Res2 = ( MP )->evallog_mat ( Grid ); 
    6262 
    6363        it_file it ( "merger_2d_test.it" ); 
     
    6666        it << Name ( "Res1" ) << Res1; 
    6767        it << Name ( "Res2" ) << Res2; 
    68         it << Name ( "S1" ) << f1->evallog_m ( Grid ); 
    69         it << Name ( "S2" ) << f2->evallog_m ( Grid ); 
     68        it << Name ( "S1" ) << f1->evallog_mat ( Grid ); 
     69        it << Name ( "S2" ) << f2->evallog_mat ( Grid ); 
    7070        cout << ( ( enorm<ldmat>* ) ( MP->_Coms ( 0 ).get() ) )->_R().to_mat() << endl; 
    7171} 
  • library/tests/merger_iter_test.cpp

    r693 r713  
    6060        epdf* MP = MM.epredictor();//xy 
    6161 
    62         vec Res1 = M.evallog_m ( Grid ); 
    63         mat Res2 = ( ( emix* ) MP )->evallog_M ( Grid ); 
     62        vec Res1 = M.evallog_mat ( Grid ); 
     63        mat Res2 = ( ( emix* ) MP )->evallog_mat ( Grid ); 
    6464 
    6565        it_file it ( "merger_iter_test.it" ); 
  • library/tests/merger_test.cpp

    r693 r713  
    2929        x_grid.set_row ( 0, linspace ( -10.0, 10.0 ) ); 
    3030 
    31         vec l_f1 = f1->evallog_m ( x_grid ); 
    32         vec l_f2 = f2->evallog_m ( x_grid ); 
     31        vec l_f1 = f1->evallog_mat ( x_grid ); 
     32        vec l_f2 = f2->evallog_mat ( x_grid ); 
    3333        mat lW ( 2, Npoints ); 
    3434        lW.set_row ( 0, l_f1 ); 
  • library/tests/mixtures_test.cpp

    r679 r713  
    6969        // Sample parameters 
    7070        int ndat = 100; 
    71         mat Smp = Simul.sample_m ( ndat ); 
     71        mat Smp = Simul.sample_mat ( ndat ); 
    7272 
    7373        cout << "Simulated means: " << m1 << " and " << m2 << endl; 
  • library/tests/pdf_harness.cpp

    r706 r713  
    5050        Array<vec> actual(CurrentContext::max_trial_count); 
    5151        do { 
    52                 mat smp = hpdf->samplecond_m ( cond, nsamples ); 
     52                mat smp = hpdf->samplecond_mat ( cond, nsamples ); 
    5353                vec emu = smp * ones ( nsamples ) / nsamples; 
    5454                actual( tc ) = emu; 
     
    7171        Array<mat> actual(CurrentContext::max_trial_count); 
    7272        do { 
    73                 mat smp = hpdf->samplecond_m ( cond, nsamples ); 
     73                mat smp = hpdf->samplecond_mat ( cond, nsamples ); 
    7474                vec emu = smp * ones ( nsamples ) / nsamples; 
    7575                mat er = ( smp * smp.T() ) / nsamples - outer_product ( emu, emu ); 
  • library/tests/pdf_test.cpp

    r706 r713  
    1313static void check_covariance(mmix &distrib_obj, const vec &mu0, int nsamples, const mat &R, double tolerance); 
    1414 
    15 TEST ( mepdf_test ) { 
    16         pdf_harness::test_config ( "mepdf.cfg" ); 
    17 } 
    18  
    1915TEST ( mgamma_test ) { 
    2016        pdf_harness::test_config ( "mgamma.cfg" ); 
     
    2925} 
    3026 
    31 // not using pdf_harness because mmix isn't configurable (yet?) 
    3227TEST ( mmix_test ) { 
     28        pdf_harness::test_config ( "mmix.cfg" ); 
     29        //vec tmu = 0.5 * eN->mean() + 0.5 * mu0; 
     30 
     31        /* 
    3332        RV x ( "{mmixx }", "2" ); 
    3433        RV y ( "{mmixy }", "2" ); 
     
    6362        mat observedR ( "1.27572 0.778247; 0.778247 3.33129" ); 
    6463        check_covariance( mMix, mu0, N, observedR, 0.2 ); 
     64*/ 
    6565} 
    6666 
    6767// not using pdf_harness because emix isn't configurable (yet?) 
    68 TEST ( mepdf_emix_test ) { 
     68TEST ( emix_test ) { 
     69        pdf_harness::test_config ( "emix.cfg" ); 
     70        /* 
    6971        int N = 10000; //number of samples 
    7072        vec mu0 ( "1.5 1.7" ); 
     
    8789        eMix->set_parameters ( vec_2 ( 0.5, 0.5 ), Coms ); 
    8890        check_mean ( *eMix, mu0, N, eMix->mean(), 0.1 ); 
     91        */ 
    8992} 
    9093 
     
    9396        Array<vec> actual(CurrentContext::max_trial_count); 
    9497        do { 
    95                 mat smp = distrib_obj.samplecond_m ( mu0, nsamples ); 
     98                mat smp = distrib_obj.samplecond_mat ( mu0, nsamples ); 
    9699                vec emu = smp * ones ( nsamples ) / nsamples ; 
    97100                actual( tc ) = emu; 
     
    113116        Array<mat> actual(CurrentContext::max_trial_count); 
    114117        do { 
    115                 mat smp = distrib_obj.samplecond_m ( mu0, nsamples ); 
     118                mat smp = distrib_obj.samplecond_mat ( mu0, nsamples ); 
    116119                vec emu = smp * ones ( nsamples ) / nsamples ; 
    117120                mat er = ( smp * smp.T() ) / nsamples - outer_product ( emu, emu );