Changeset 737 for library/tests

Show
Ignore:
Timestamp:
11/25/09 12:14:38 (14 years ago)
Author:
mido
Message:

ASTYLER RUN OVER THE WHOLE LIBRARY, JUPEE

Location:
library/tests
Files:
6 modified

Legend:

Unmodified
Added
Removed
  • library/tests/epdf_harness.cpp

    r730 r737  
    4545        if ( support ) { // support is given 
    4646                grid_fnc ep_disc; 
    47                 ep_disc.set_support(*support); 
    48                 ep_disc.set_values(*hepdf); 
     47                ep_disc.set_support ( *support ); 
     48                ep_disc.set_values ( *hepdf ); 
    4949                // ep_disc is discretized at support points 
    5050 
    51                 double point_volume =prod(support->_steps()); 
    52                 CHECK_CLOSE(1.0, sum(ep_disc._values())*point_volume, 0.01); 
    53  
    54                 vec pdf=ep_disc._values(); 
    55                 pdf /=sum(pdf); // normalize 
    56                  
    57                 vec mea=pdf(0) * support->first_vec(); 
    58                 mat Remp=pdf(0) * outer_product(support->act_vec(), support->act_vec()); 
     51                double point_volume = prod ( support->_steps() ); 
     52                CHECK_CLOSE ( 1.0, sum ( ep_disc._values() ) *point_volume, 0.01 ); 
     53 
     54                vec pdf = ep_disc._values(); 
     55                pdf /= sum ( pdf ); // normalize 
     56 
     57                vec mea = pdf ( 0 ) * support->first_vec(); 
     58                mat Remp = pdf ( 0 ) * outer_product ( support->act_vec(), support->act_vec() ); 
    5959 
    6060                // run through all points 
    61                 for (int i=1; i<support->points(); i++){ 
    62                         mea += pdf(i)*support->next_vec(); 
    63                         Remp += pdf(i) * outer_product(support->act_vec(), support->act_vec()); 
    64                 } 
    65                 CHECK_CLOSE(mean, mea, tolerance); 
    66                 CHECK_CLOSE(R, Remp-outer_product(mea,mea), tolerance); 
     61                for ( int i = 1; i < support->points(); i++ ) { 
     62                        mea += pdf ( i ) * support->next_vec(); 
     63                        Remp += pdf ( i ) * outer_product ( support->act_vec(), support->act_vec() ); 
     64                } 
     65                CHECK_CLOSE ( mean, mea, tolerance ); 
     66                CHECK_CLOSE ( R, Remp - outer_product ( mea, mea ), tolerance ); 
    6767        } 
    6868 
     
    7070                check_sample_mean(); 
    7171        } 
    72         if (R.rows()>0){ 
     72        if ( R.rows() > 0 ) { 
    7373                check_covariance(); 
    7474        } 
     
    133133        } while ( ( tc < CurrentContext::max_trial_count ) && 
    134134                  !UnitTest::AreClose ( mean, actual ( tc - 1 ), delta ) ); 
    135                            
     135 
    136136        if ( ( tc == CurrentContext::max_trial_count ) && 
    137137                ( !UnitTest::AreClose ( mean, actual ( CurrentContext::max_trial_count - 1 ), delta ) ) ) { 
  • library/tests/testsuite/LQG_test.cpp

    r723 r737  
    55using namespace bdm; 
    66 
    7 TEST(LQG_test) { 
     7TEST ( LQG_test ) { 
    88        LQG reg; 
    9         shared_ptr<StateSpace<chmat> > stsp=new StateSpace<chmat>; 
     9        shared_ptr<StateSpace<chmat> > stsp = new StateSpace<chmat>; 
    1010        // 2 x 1 x 1 
    11         stsp-> set_parameters(eye(2), ones(2,1), ones(1,2), ones(1,1), /* Q,R */ eye(2), eye(1)); 
    12         reg.set_system(stsp); // A, B, C 
    13         reg.set_control_parameters(eye(1), eye(1),  vec_1(1.0), 6); //Qy, Qu, horizon 
     11        stsp-> set_parameters ( eye ( 2 ), ones ( 2, 1 ), ones ( 1, 2 ), ones ( 1, 1 ), /* Q,R */ eye ( 2 ), eye ( 1 ) ); 
     12        reg.set_system ( stsp ); // A, B, C 
     13        reg.set_control_parameters ( eye ( 1 ), eye ( 1 ),  vec_1 ( 1.0 ), 6 ); //Qy, Qu, horizon 
    1414        reg.validate(); 
    1515 
    1616        reg.redesign(); 
    17         double reg_apply=reg.ctrlaction("0.5, 1.1","0.0")(0); /*convert vec to double*/ 
    18         CHECK_CLOSE(reg_apply, -0.248528137234392, 0.0001); 
     17        double reg_apply = reg.ctrlaction ( "0.5, 1.1", "0.0" ) ( 0 ); /*convert vec to double*/ 
     18        CHECK_CLOSE ( reg_apply, -0.248528137234392, 0.0001 ); 
    1919} 
    2020 
    21 TEST(to_state_test) { 
     21TEST ( to_state_test ) { 
    2222        mlnorm<fsqmat> ml; 
    23         mat A="1.1, 2.3"; 
    24         ml.set_parameters(A, vec_1(1.3), eye(1)); 
    25         RV yr=RV("y",1); 
    26         RV ur=RV("u",1); 
    27         ml.set_rv(yr); 
    28         yr.t_plus(-1); 
    29         ml.set_rvc(concat(yr, ur)); 
    30          
    31         shared_ptr<StateCanonical > Stsp=new StateCanonical; 
    32         Stsp->connect_mlnorm(ml); 
    33          
    34         /* results from  
     23        mat A = "1.1, 2.3"; 
     24        ml.set_parameters ( A, vec_1 ( 1.3 ), eye ( 1 ) ); 
     25        RV yr = RV ( "y", 1 ); 
     26        RV ur = RV ( "u", 1 ); 
     27        ml.set_rv ( yr ); 
     28        yr.t_plus ( -1 ); 
     29        ml.set_rvc ( concat ( yr, ur ) ); 
     30 
     31        shared_ptr<StateCanonical > Stsp = new StateCanonical; 
     32        Stsp->connect_mlnorm ( ml ); 
     33 
     34        /* results from 
    3535        [A,B,C,D]=tf2ss([2.3 0],[1 -1.1]) 
    3636        */ 
    37         CHECK_CLOSE_EX(Stsp->_A().get_row(0), vec("1.1"), 0.0001); 
    38         CHECK_CLOSE_EX(Stsp->_C().get_row(0), vec("2.53"), 0.0001); 
    39         CHECK_CLOSE_EX(Stsp->_D().get_row(0), vec("2.30"), 0.0001); 
     37        CHECK_CLOSE_EX ( Stsp->_A().get_row ( 0 ), vec ( "1.1" ), 0.0001 ); 
     38        CHECK_CLOSE_EX ( Stsp->_C().get_row ( 0 ), vec ( "2.53" ), 0.0001 ); 
     39        CHECK_CLOSE_EX ( Stsp->_D().get_row ( 0 ), vec ( "2.30" ), 0.0001 ); 
    4040} 
    4141 
    42 TEST(to_state_arx_test) { 
     42TEST ( to_state_arx_test ) { 
    4343        mlnorm<chmat> ml; 
    44         mat A="1.1, 2.3, 3.4"; 
    45         ml.set_parameters(A, vec_1(1.3), eye(1)); 
    46         RV yr=RV("y",1); 
    47         RV ur=RV("u",1); 
    48         ml.set_rv(yr); 
    49         ml.set_rvc(concat(yr.copy_t(-1), concat(ur, ur.copy_t(-1)))); 
    50          
    51         shared_ptr<StateFromARX> Stsp=new StateFromARX; 
    52         RV xrv; RV urv; 
    53         Stsp->connect_mlnorm(ml,xrv,urv); 
    54          
    55         /* results from  
     44        mat A = "1.1, 2.3, 3.4"; 
     45        ml.set_parameters ( A, vec_1 ( 1.3 ), eye ( 1 ) ); 
     46        RV yr = RV ( "y", 1 ); 
     47        RV ur = RV ( "u", 1 ); 
     48        ml.set_rv ( yr ); 
     49        ml.set_rvc ( concat ( yr.copy_t ( -1 ), concat ( ur, ur.copy_t ( -1 ) ) ) ); 
     50 
     51        shared_ptr<StateFromARX> Stsp = new StateFromARX; 
     52        RV xrv; 
     53        RV urv; 
     54        Stsp->connect_mlnorm ( ml, xrv, urv ); 
     55 
     56        /* results from 
    5657        [A,B,C,D]=tf2ss([2.3 0],[1 -1.1]) 
    5758        */ 
    58         CHECK_CLOSE_EX(Stsp->_A().get_row(0), vec("1.1, 3.4, 1.3"), 0.0001); 
    59         CHECK_CLOSE_EX(Stsp->_B().get_row(0), vec("2.3"), 0.0001); 
    60         CHECK_CLOSE_EX(Stsp->_C().get_row(0), vec("1, 0, 0"), 0.0001); 
     59        CHECK_CLOSE_EX ( Stsp->_A().get_row ( 0 ), vec ( "1.1, 3.4, 1.3" ), 0.0001 ); 
     60        CHECK_CLOSE_EX ( Stsp->_B().get_row ( 0 ), vec ( "2.3" ), 0.0001 ); 
     61        CHECK_CLOSE_EX ( Stsp->_C().get_row ( 0 ), vec ( "1, 0, 0" ), 0.0001 ); 
    6162} 
    6263 
    63 TEST(arx_LQG_test){ 
     64TEST ( arx_LQG_test ) { 
    6465        mlnorm<chmat> ml; 
    65         mat A="1.81, -.81, .00468, .00438"; 
    66         ml.set_parameters(A, vec_1(0.0), 0.00001*eye(1)); 
    67         RV yr=RV("y",1); 
    68         RV ur=RV("u",1); 
    69         RV rgr = yr.copy_t(-1); 
    70         rgr.add(yr.copy_t(-2)); 
    71         rgr.add(yr.copy_t(-2)); 
    72         rgr.add(ur.copy_t(-1)); 
    73         rgr.add(ur); 
    74          
    75         ml.set_rv(yr); 
    76         ml.set_rvc(rgr); 
     66        mat A = "1.81, -.81, .00468, .00438"; 
     67        ml.set_parameters ( A, vec_1 ( 0.0 ), 0.00001*eye ( 1 ) ); 
     68        RV yr = RV ( "y", 1 ); 
     69        RV ur = RV ( "u", 1 ); 
     70        RV rgr = yr.copy_t ( -1 ); 
     71        rgr.add ( yr.copy_t ( -2 ) ); 
     72        rgr.add ( yr.copy_t ( -2 ) ); 
     73        rgr.add ( ur.copy_t ( -1 ) ); 
     74        rgr.add ( ur ); 
     75 
     76        ml.set_rv ( yr ); 
     77        ml.set_rvc ( rgr ); 
    7778        ml.validate(); 
    78          
    79         shared_ptr<StateFromARX> Stsp=new StateFromARX; 
    80         RV xrv; RV urv; 
    81         Stsp->connect_mlnorm(ml,xrv,urv); 
    82          
     79 
     80        shared_ptr<StateFromARX> Stsp = new StateFromARX; 
     81        RV xrv; 
     82        RV urv; 
     83        Stsp->connect_mlnorm ( ml, xrv, urv ); 
     84 
    8385        LQG L; 
    84         L.set_system(Stsp); 
    85         L.set_control_parameters(eye(1), sqrt(1.0/1000)*eye(1), vec_1(0.0), 100); 
     86        L.set_system ( Stsp ); 
     87        L.set_control_parameters ( eye ( 1 ), sqrt ( 1.0 / 1000 ) *eye ( 1 ), vec_1 ( 0.0 ), 100 ); 
    8688        L.validate(); 
    87          
     89 
    8890        L.redesign(); 
    89         cout << L.to_string()<<endl; 
     91        cout << L.to_string() << endl; 
    9092} 
  • library/tests/testsuite/datalink_test.cpp

    r723 r737  
    101101 
    102102        vec val_up ( "1 2" ); 
    103         dl.store_data(val_up); 
     103        dl.store_data ( val_up ); 
    104104        val_up = "3 4"; 
    105105 
  • library/tests/testsuite/logger_test.cpp

    r728 r737  
    4242 
    4343        int pol = 11; 
    44          
     44 
    4545        int rid = logger.add_vector ( r, "" ); 
    4646        int thid = logger.add_vector ( th, "th" ); 
    47         int polid = logger.add_setting("pols"); 
     47        int polid = logger.add_setting ( "pols" ); 
    4848 
    4949        logger.init(); 
     
    5252                logger.log_vector ( rid, vec_2 ( ( double ) i, ( double ) ( i + 1 ) ) ); 
    5353                logger.log_vector ( thid, vec_2 ( ( double ) ( 100 - i ), ( double ) ( i - 50 ) ) ); 
    54                 UI::save(pol, logger.log_to_setting(polid), "polozka"); 
     54                UI::save ( pol, logger.log_to_setting ( polid ), "polozka" ); 
    5555                logger.step(); 
    5656        } 
     
    5858        logger.finalize(); 
    5959 
    60         CHECK_EQUAL ( get_file_size ( "logger_test.matrix" ), get_file_size ( (ls+".it").c_str() ) ); 
    61         CHECK_EQUAL ( get_file_size ( (ls+".cfg.check").c_str() ), get_file_size ( (ls+".cfg" ).c_str() )); 
     60        CHECK_EQUAL ( get_file_size ( "logger_test.matrix" ), get_file_size ( ( ls + ".it" ).c_str() ) ); 
     61        CHECK_EQUAL ( get_file_size ( ( ls + ".cfg.check" ).c_str() ), get_file_size ( ( ls + ".cfg" ).c_str() ) ); 
    6262} 
    6363 
  • library/tests/unittest-cpp/TestRunner.cpp

    r736 r737  
    7171 
    7272void TestRunner::ClearEnvironment() const { 
    73         bdm::RV::clear_all();    
     73        bdm::RV::clear_all(); 
    7474} 
    7575 
  • library/tests/unittest-cpp/TestRunner.h

    r736 r737  
    3131 
    3232                while ( curTest != 0 ) { 
    33                         if ( IsTestInSuite ( curTest, suiteName ) && predicate ( curTest ) ) 
    34                         { 
     33                        if ( IsTestInSuite ( curTest, suiteName ) && predicate ( curTest ) ) { 
    3534                                ClearEnvironment(); 
    3635                                RunTest ( m_result, curTest, maxTestTimeInMs ); 
     
    5150        int Finish() const; 
    5251        bool IsTestInSuite ( const Test* const curTest, char const* suiteName ) const; 
    53         void RunTest ( TestResults* const result, Test* const curTest, int const maxTestTimeInMs ) const;        
     52        void RunTest ( TestResults* const result, Test* const curTest, int const maxTestTimeInMs ) const; 
    5453        void ClearEnvironment() const; 
    5554};