Changeset 466

Timestamp:

08/03/09 10:56:48 (15 years ago)

Author:

vbarta

Message:

using computed tolerance also when sampling mprod

Files:

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

library/tests/epdf_harness.cpp

@@ -85 +85 @@
     if (R.rows() > 0) {
         mat smp = hepdf->sample_m(nsamples);
-        int n = smp.cols();
-        vec emu = smp * ones(n) / n;
-        mat er = (smp * smp.T()) / n - outer_product(emu, emu);
+        vec emu = smp * ones(nsamples) / nsamples;
+        mat er = (smp * smp.T()) / nsamples - outer_product(emu, emu);
 
         vec delta = sqrt(variance) / sqrt(nsamples);
@@ -106 +105 @@
         mprod mEp(aa);
 
-        int n = nsamples;
-        mat smp = mEp.samplecond(vec(0), n);
-        vec emu = sum(smp, 2) / n;
-        CHECK_CLOSE_EX(mean, emu, tolerance);
+        mat smp = mEp.samplecond(vec(0), nsamples);
+        vec emu = sum(smp, 2) / nsamples;
+
+        vec delta = sqrt(variance) / sqrt(nsamples);
+        CHECK_CLOSE_EX(mean, emu, delta);
 
         if (R.rows() > 0) {
-            mat er = (smp * smp.T()) / n - outer_product(emu, emu);
+            mat er = (smp * smp.T()) / nsamples - outer_product(emu, emu);
             CHECK_CLOSE_EX(R, er, tolerance);
         }