Changeset 536 for library/tests
- Timestamp:
- 08/16/09 18:13:31 (15 years ago)
- Location:
- library/tests
- Files:
-
- 2 modified
Legend:
- Unmodified
- Added
- Removed
-
library/tests/arx_elem_test.cpp
r477 r536 12 12 mat mu ( 1, 1 ); 13 13 mat R ( 1, 1 ); 14 Ar. _e()->mean_mat ( mu, R );14 Ar.posterior().mean_mat ( mu, R ); 15 15 cout << "Prior moments: mu=" << mu << ", R=" << R << endl; 16 16 … … 24 24 // Ar is now filled with estimates of N(0,1); 25 25 cout << "Empirical moments: mu=" << sum ( smp ) / ndat << ", R=" << sum_sqr ( smp ) / ndat - pow ( sum ( smp ) / ndat, 2 ) << endl; 26 Ar. _e()->mean_mat ( mu, R );26 Ar.posterior().mean_mat ( mu, R ); 27 27 cout << "Posterior moments: mu=" << mu << ", R=" << R << endl; 28 28 -
library/tests/epdf_harness.cpp
r529 r536 210 210 Array<vec> actual(CurrentContext::max_trial_count); 211 211 do { 212 mat smp = mep.samplecond ( vec ( 0 ), nsamples );212 mat smp = mep.samplecond_m ( vec ( 0 ), nsamples ); 213 213 vec emu = sum ( smp, 2 ) / nsamples; 214 214 actual( tc ) = emu; … … 231 231 Array<mat> actual(CurrentContext::max_trial_count); 232 232 do { 233 mat smp = mep.samplecond ( vec ( 0 ), nsamples );233 mat smp = mep.samplecond_m ( vec ( 0 ), nsamples ); 234 234 vec emu = sum ( smp, 2 ) / nsamples; 235 235 mat er = ( smp * smp.T() ) / nsamples - outer_product ( emu, emu );