Changeset 713 for library/tests/pdf_test.cpp
- Timestamp:
- 11/10/09 11:03:09 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
library/tests/pdf_test.cpp
r706 r713 13 13 static void check_covariance(mmix &distrib_obj, const vec &mu0, int nsamples, const mat &R, double tolerance); 14 14 15 TEST ( mepdf_test ) {16 pdf_harness::test_config ( "mepdf.cfg" );17 }18 19 15 TEST ( mgamma_test ) { 20 16 pdf_harness::test_config ( "mgamma.cfg" ); … … 29 25 } 30 26 31 // not using pdf_harness because mmix isn't configurable (yet?)32 27 TEST ( mmix_test ) { 28 pdf_harness::test_config ( "mmix.cfg" ); 29 //vec tmu = 0.5 * eN->mean() + 0.5 * mu0; 30 31 /* 33 32 RV x ( "{mmixx }", "2" ); 34 33 RV y ( "{mmixy }", "2" ); … … 63 62 mat observedR ( "1.27572 0.778247; 0.778247 3.33129" ); 64 63 check_covariance( mMix, mu0, N, observedR, 0.2 ); 64 */ 65 65 } 66 66 67 67 // not using pdf_harness because emix isn't configurable (yet?) 68 TEST ( mepdf_emix_test ) { 68 TEST ( emix_test ) { 69 pdf_harness::test_config ( "emix.cfg" ); 70 /* 69 71 int N = 10000; //number of samples 70 72 vec mu0 ( "1.5 1.7" ); … … 87 89 eMix->set_parameters ( vec_2 ( 0.5, 0.5 ), Coms ); 88 90 check_mean ( *eMix, mu0, N, eMix->mean(), 0.1 ); 91 */ 89 92 } 90 93 … … 93 96 Array<vec> actual(CurrentContext::max_trial_count); 94 97 do { 95 mat smp = distrib_obj.samplecond_m ( mu0, nsamples );98 mat smp = distrib_obj.samplecond_mat ( mu0, nsamples ); 96 99 vec emu = smp * ones ( nsamples ) / nsamples ; 97 100 actual( tc ) = emu; … … 113 116 Array<mat> actual(CurrentContext::max_trial_count); 114 117 do { 115 mat smp = distrib_obj.samplecond_m ( mu0, nsamples );118 mat smp = distrib_obj.samplecond_mat ( mu0, nsamples ); 116 119 vec emu = smp * ones ( nsamples ) / nsamples ; 117 120 mat er = ( smp * smp.T() ) / nsamples - outer_product ( emu, emu );