Changeset 441 for library/tests/epdf_harness.cpp

Timestamp:

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

Author:

vbarta

Message:

added configurable tests for egamma

Files:

• library/tests/epdf_harness.cpp (modified) (3 diffs)

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);
+    }
 }