Changeset 162 for tests

Timestamp:

09/04/08 20:27:01 (16 years ago)

Author:

smidl

Message:

opravy a dokumentace

Location:

tests

Files:

• CMakeLists.txt (modified) (1 diff)
• arx_test.cpp (modified) (1 diff)
• loggers_test.cpp (modified) (2 diffs)
• rv_test.cpp (modified) (1 diff)
• test0.cpp (modified) (1 diff)
• testEpdf.cpp (modified) (1 diff)
• testKF.cpp (modified) (1 diff)
• testKF_QR.cpp (modified) (1 diff)
• testKF_QRexh.cpp (modified) (1 diff)
• testSmp.cpp (modified) (2 diffs)

tests/CMakeLists.txt

@@ -38 +38 @@
 target_link_libraries (rv_test ${BdmLibs})
 
+add_executable (merger_test merger_test.cpp)
+target_link_libraries (merger_test ${BdmLibs})
+
 add_executable (test0 test0.cpp)
 add_executable (testKF testKF.cpp)

tests/arx_test.cpp

@@ -17 +17 @@
         double nu0 = ord+1;
         
-        RV thr("1","{theta_and_r }",vec_1(th.length()+1),"0");
+        RV thr("{theta_and_r }");
         ARX Ar(thr, V0, nu0);
         epdf& Ar_ep = Ar._epdf();

tests/loggers_test.cpp

@@ -12 +12 @@
         int i;
         
-        RV th = RV ( "1 2","{a b }","1 1","0 0");
-        RV r = RV ( "2","{r }","2","0" );
+        RV th = RV ( "{a b }");
+        RV r = RV ( "{r }","2");
 
         dirfilelog L("exp/dirfile",10);
@@ -30 +30 @@
                 L.step();
         }
-        L.step(true);
+        L.finalize();
 }

tests/rv_test.cpp

@@ -12 +12 @@
 int main()
 {
-        RV a = RV ( "1","{a }","3","0");
-        RV b = RV ( "2","{b }","2","0");
-        RV c = RV ( "3","{c }","1","0");
-        RV ef = RV ( "4 5","{e f }", " 1 1", "0 0" );
+        RV a = RV ( "{a }","3");
+        RV b = RV ( "{b }","2");
+        RV c = RV ( "{c }");
+        RV ef = RV ("{e f }");
         
         cout << "Add a and b"<<endl;

tests/test0.cpp

@@ -12 +12 @@
 int main()
 {
-        RV th = RV ( "1 2","{a b }","1 1","0 0");
-        RV r = RV ( "3 4" );
+        RV th = RV ( "{a b }");
+        RV r = RV ( "{r1 r2 }" );
 
         cout << th << r << endl;
 
-        ldmat ld = ldmat("1 0;0.5 1","1.1 1.3");
-        vec v = "1 -0.1";
+        ldmat ld = ldmat("1 0 0 0; 0.2 1 0 0; 0.1 -0.1 1 0; -0.9 0.3 -0.1 1",
+                                         "3 4 5 6");
+        vec v = "1 -1 3 -1";
         
         cout << "Test to_mat"<<endl;
 //      cout << ld << endl;
-        cout << ld.to_mat()<<endl;
+        mat L = ld.to_mat();
+        cout << L <<endl;
 
         cout << "Test opupdt" <<endl;
+        cout << "w=1" <<endl;
         ldmat ldup = ld;
         ldup.opupdt(v,1);
         cout << ldup.to_mat()<<endl;
+        cout << L + outer_product(v,v) <<endl;
+        ldup = ld;      
+        cout << "w=0.1" <<endl;
+        ldup.opupdt(v,0.1);
+        cout << ldup.to_mat()<<endl;
+        cout << L + 0.1*outer_product(v,v) <<endl;
 
         cout << "Test +="<<endl;
         ldmat ld2x = ld;
         ld2x+=ld;
-        cout << ld.to_mat() << ld2x.to_mat() <<endl;
+        cout << ld.to_mat() <<endl << ld2x.to_mat() <<endl;
 
+        cout << "Test *="<<endl;
+        ld2x = ld;
+        ld2x*=2;
+        cout << ld.to_mat() <<endl << ld2x.to_mat() <<endl;
+        
         cout << "Test ldinv()"<<endl;
         ldmat Il = ld; 

tests/testEpdf.cpp

@@ -14 +14 @@
         vec y = "-10:0.1:10";
 
-        RV rv ( "1","{x }","2","0" );
+        RV rv ( "{x }","2" );
         vec mu0 = "0.0 0.0";
         mat V0 ( "5 -0.05; -0.05 5.20" );

tests/testKF.cpp

@@ -41 +41 @@
         
         // Prepare for Kalman filters in BDM:
-        ivec tmpsize(1);
-        tmpsize(0) = A.cols();
-        RV rx("1","{x }",tmpsize,"0");
-        tmpsize(0) = B.cols();
-        RV ru("2","{u }",tmpsize,"0");
-        tmpsize(0) = C.rows();
-        RV ry("3","{y }",tmpsize,"0");
+        RV rx("{x }",vec_1(A.cols()));
+        RV ru("{u }",vec_1(B.cols()));
+        RV ry("{y }",vec_1(C.rows()));
         
 //      // LDMAT

tests/testKF_QR.cpp

@@ -43 +43 @@
 
 //      cout << KF;
-        RV rx("1","{x}","2","0");
-        RV ru("2","{u}","1","0");
-        RV ry("3","{y}","1","0");
-        RV rQR("4","{Q,R}","3","0");
+        RV rx("{x }","2");
+        RV ru("{u }","1");
+        RV ry("{y }","1");
+        RV rQR("{Q,R }","3");
         //
         KFcondQR KF(rx,ry,ru,rQR);

tests/testKF_QRexh.cpp

@@ -48 +48 @@
         Ndat = Dt.cols();
 
-        RV rx ( "1","{x}","2","0" );
-        RV ru ( "2","{u}","1","0" );
-        RV ry ( "3","{y}","2","0" );
+        RV rx ( "{x}","2" );
+        RV ru ( "{u}","1" );
+        RV ry ( "{y}","2" );
         //
         //

tests/testSmp.cpp

@@ -24 +24 @@
         RNG_randomize();
 
-        RV x("1","{x }","2","0");
-        RV y("2","{y }","2","0");
+        RV x("{x }","2");
+        RV y("{y }","2");
         int N = 10000; //number of samples
         vec mu0 = "1.5 1.7";
@@ -101 +101 @@
         eG._rv()= y; 
         //create array 
-        Array<epdf*> A(2);
-        A(0) = &eN;
-        A(1) = &eG;
+        Array<mpdf*> A(2);
+        mepdf meN(eN);
+        mepdf meG(eG);
+        A(0) = &meN;
+        A(1) = &meG;
         
-        eprod eP(A);
+        mprod eP(A);
         mat epV=zeros(4,4);
         epV.set_submatrix(0,0,V0);
         epV.set_submatrix(2,2,diag(g_var));
-        Smp = eP.sampleN(N);
-        disp(eP.mean(), epV,Smp);
+        
+        vec v0=vec(0);
+        Smp = eP.samplecond(v0,lik,N);
+        disp(concat(eN.mean(),eG.mean()), epV,Smp);
         
         //Exit program: