Changeset 294 for tests

Timestamp:

03/12/09 22:18:43 (16 years ago)

Author:

smidl

Message:

tr2245

Location:

tests

Files:

• chmat_test.cpp (modified) (2 diffs)
• testSmp.cpp (modified) (1 diff)

tests/chmat_test.cpp

@@ -17 +17 @@
    cout << "Testing constructors:" << endl
        << "A = " << A << endl
+       << "Ch = " << Ch._Ch() << endl
        << "Ch.to_mat() = " << Ch.to_mat() << endl << endl;
        
@@ -49 +50 @@
    cout << "A+vv' = " << A+outer_product(v,v) << endl
                         <<      "opupdt(Ch,v) = " << Ch2.to_mat() << endl <<endl;
-                                
+
+   vec vCh=vec(Ch._Ch()._data(),9);
+   cout << "vectorized Ch: " << vCh <<endl;
+   chmat nCh(3);
+   nCh.setCh(vCh);
+   cout << "Ch from vec: " << nCh._Ch() <<endl;
+
 }

tests/testSmp.cpp

@@ -116 +116 @@
         disp(concat(eN.mean(),eG.mean()), epV,Smp);
         
+        cout << "======= eWishart ======== " << endl;
+        mat wM="1.0 0.9; 0.9 1.0";
+        eWishartCh eW; eW.set_parameters(wM/100,100);
+        mat mea=zeros(2,2);
+        mat Ch(2,2);
+        for (int i=0;i<100;i++){Ch=eW.sample_mat(); mea+=Ch.T()*Ch;}
+        cout << mea /100 <<endl;
+        
+        cout << "======= rwiWishart ======== " << endl;
+        rwiWishartCh rwW; rwW.set_parameters(2,0.1,"1 1",0.9);
+        mea=zeros(2,2);
+        mat wMch=chol(wM);
+        for (int i=0;i<100;i++){ 
+                vec tmp=rwW.samplecond(vec(wMch._data(),4));
+                 copy_vector(4,tmp._data(), Ch._data()); 
+                 mea+=Ch.T()*Ch;
+        }
+        cout << mea /100 <<endl;
         //Exit program:
         return 0;