Changeset 441

Timestamp:

07/29/09 08:29:35 (16 years ago)

Author:

vbarta

Message:

added configurable tests for egamma

Location:

library/tests

Files:

• CMakeLists.txt (modified) (1 diff)
• egamma.cfg (added)
• egamma_test.cpp (added)
• egiw.cfg (modified) (1 diff)
• epdf_harness.cpp (modified) (3 diffs)
• epdf_harness.h (modified) (1 diff)

library/tests/CMakeLists.txt

@@ -38 +38 @@
 
 # using UnitTest++
-add_executable(testsuite datalink_test.cpp egiw_test.cpp enorm_test.cpp loggers_test.cpp rv_test.cpp square_mat_test.cpp testsuite.cpp test_user_info.cpp test_shared_ptr.cpp)
+add_executable(testsuite datalink_test.cpp egamma_test.cpp egiw_test.cpp enorm_test.cpp loggers_test.cpp rv_test.cpp square_mat_test.cpp testsuite.cpp test_user_info.cpp test_shared_ptr.cpp)
 target_link_libraries(testsuite bdm itpp testutil unittest)
 

library/tests/egiw.cfg

@@ -18 +18 @@
   nsamples = 200;
   integral = 1.0;
+},
+{
+  class = "egiw_harness";
+  epdf = {
+    class = "egiw";
+    V = ( "matrix", 2, 2, [ 20.0, 8.0, 8.0, 4.0 ] );
+    nu = 4.0;
+    dimx = 1.0;
+    rv : 
+    {
+      class = "RV";
+      names = ( "a", "b" );
+    };
+  };
+  mean = [ 2, -4 ];
+  lognc = 0.451583;
+  variance = [ -1, -16 ];
 } );
 

library/tests/epdf_harness.cpp

@@ -44 +44 @@
     UI::get(mean, set, "mean");
     UI::get(variance, set, "variance");
-    UI::get(support, set, "support");
+
+    if (set.exists("support")) {
+        UI::get(support, set, "support");
+    }
 
     if (set.exists("nsamples")) {
@@ -50 +53 @@
     }
 
-    UI::get(integral, set, "integral");
+    if (set.exists("integral")) {
+        UI::get(integral, set, "integral");
+    }
+
+    if (set.exists("R")) {
+        UI::get(R, set, "R");
+    }
 }
 
@@ -58 +67 @@
     CHECK_CLOSE(variance, hepdf->variance(), epsilon);
 
-    vec xb = support.get_row(0);
-    vec yb = support.get_row(1);
-    CHECK_CLOSE(mean, num_mean2(hepdf.get(), xb, yb, nsamples, nsamples), 0.1);
-    CHECK_CLOSE(integral, normcoef(hepdf.get(), xb, yb, nsamples, nsamples), 0.1);
+    if (support.rows() == 2) {
+        vec xb = support.get_row(0);
+        vec yb = support.get_row(1);
+        CHECK_CLOSE(mean, num_mean2(hepdf.get(), xb, yb, nsamples, nsamples), 0.1);
+        CHECK_CLOSE(integral, normcoef(hepdf.get(), xb, yb, nsamples, nsamples), 0.1);
+    }
+
+    if (R.rows() > 0) {
+        mat Smp = hepdf->sample_m(10000);
+        int N = Smp.cols();
+        vec Emu = Smp*ones(N) /N ;
+        mat Er = (Smp*Smp.transpose())/N - outer_product(Emu,Emu);
+        CHECK_CLOSE(mean, Emu, 0.1);
+        CHECK_CLOSE(R, Er, 0.1);
+    }
 }
 

library/tests/epdf_harness.h

@@ -34 +34 @@
     int nsamples;
     double integral;
+    mat R;
 
 public: