Changeset 524

Timestamp:

08/13/09 12:23:30 (15 years ago)

Author:

vbarta

Message:

moved test of mepdf of emix from testSmp to mpdf_test

Files:

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

library/tests/mpdf_test.cpp

@@ -9 +9 @@
 using namespace bdm;
 
-static void check_mean(mmix &mMix, const vec &mu0, int nsamples, const vec &mean, double tolerance);
+static void check_mean(mpdf &distrib_obj, const vec &mu0, int nsamples, const vec &mean, double tolerance);
 
-static void check_covariance(mmix &mMix, const vec &mu0, int nsamples, const mat &R, double tolerance);
+static void check_covariance(mmix &distrib_obj, const vec &mu0, int nsamples, const mat &R, double tolerance);
 
 TEST ( test_mepdf ) {
@@ -68 +68 @@
 }
 
-static void check_mean(mmix &mMix, const vec &mu0, int nsamples, const vec &mean, double tolerance) {
+// not using mpdf_harness because emix isn't configurable (yet?)
+TEST ( test_mepdf_emix ) {
+        int N = 10000; //number of samples
+        vec mu0 ( "1.5 1.7" );
+        mat V0 ( "1.2 0.3; 0.3 5" );
+        ldmat R = ldmat ( V0 );
+
+        shared_ptr<enorm<ldmat> > eN = new enorm<ldmat>();
+        eN->set_parameters ( mu0, R );
+
+        vec a = "100000,10000";
+        vec b = a / 10.0;
+        shared_ptr<egamma> eG = new egamma();
+        eG->set_parameters ( a, b );
+
+        shared_ptr<emix> eMix = new emix();
+        Array<shared_ptr<epdf> > Coms ( 2 );
+        Coms ( 0 ) = eG;
+        Coms ( 1 ) = eN;
+
+        eMix->set_parameters ( vec_2 ( 0.5, 0.5 ), Coms );
+        mepdf meMix ( eMix );
+        check_mean ( meMix, mu0, N, eMix->mean(), 0.1 );
+}
+
+static void check_mean(mpdf &distrib_obj, const vec &mu0, int nsamples, const vec &mean, double tolerance) {
         int tc = 0;
         Array<vec> actual(CurrentContext::max_trial_count);
         do {
-                mat smp = mMix.samplecond_m ( mu0, nsamples );
+                mat smp = distrib_obj.samplecond_m ( mu0, nsamples );
                 vec emu = smp * ones ( nsamples ) / nsamples ;
                 actual( tc ) = emu;
@@ -88 +113 @@
 }
 
-static void check_covariance(mmix &mMix, const vec &mu0, int nsamples, const mat &R, double tolerance) {
+static void check_covariance(mmix &distrib_obj, const vec &mu0, int nsamples, const mat &R, double tolerance) {
         int tc = 0;
         Array<mat> actual(CurrentContext::max_trial_count);
         do {
-                mat smp = mMix.samplecond_m ( mu0, nsamples );
+                mat smp = distrib_obj.samplecond_m ( mu0, nsamples );
                 vec emu = smp * ones ( nsamples ) / nsamples ;
                 mat er = ( smp * smp.T() ) / nsamples - outer_product ( emu, emu );