Changeset 725 for library/tests/epdf_harness.cpp

Timestamp:

11/17/09 00:53:54 (15 years ago)

Author:

smidl

Message:

Sampling from egiw

Files:

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

library/tests/epdf_harness.cpp

@@ -59 +59 @@
         }
 
-        if ( R.rows() > 0 ) {
+        if ( variance.length() > 0 ) {
                 check_sample_mean();
+        }
+        if (R.rows()>0){
                 check_covariance();
         }
@@ -114 +116 @@
         vec yb = support.get_row ( 1 );
 
-        int tc = 0;
-        Array<vec> actual ( CurrentContext::max_trial_count );
-        do {
-                vec emu = num_mean2 ( hepdf.get(), xb, yb, nbins ( 0 ), nbins ( 1 ) );
-                actual ( tc ) = emu;
-                ++tc;
-        } while ( ( tc < CurrentContext::max_trial_count ) &&
-                  !UnitTest::AreClose ( mean, actual ( tc - 1 ), tolerance ) );
-        if ( ( tc == CurrentContext::max_trial_count ) &&
-                ( !UnitTest::AreClose ( mean, actual ( CurrentContext::max_trial_count - 1 ), tolerance ) ) ) {
-                UnitTest::MemoryOutStream stream;
-                stream << CurrentContext::format_context ( __LINE__ ) << "expected " << mean << " +/- " << tolerance << " but was " << actual;
-
-                UnitTest::TestDetails details ( *UnitTest::CurrentTest::Details(), 0, false );
-
-                UnitTest::CurrentTest::Results()->OnTestFailure ( details, stream.GetText() );
-        }
+        vec actual;
+        actual = num_mean2 ( hepdf.get(), xb, yb, nbins ( 0 ), nbins ( 1 ) );
+
+        CHECK_CLOSE(mean, actual, tolerance);
 }
 
@@ -137 +126 @@
         vec yb = support.get_row ( 1 );
 
-        int tc = 0;
-        Array<double> actual ( CurrentContext::max_trial_count );
-        do {
-                double nc = normcoef ( hepdf.get(), xb, yb, nbins ( 0 ), nbins ( 1 ) );
-                actual ( tc ) = nc;
-                ++tc;
-        } while ( ( tc < CurrentContext::max_trial_count ) &&
-                  !UnitTest::AreClose ( 1.0, actual ( tc - 1 ), tolerance ) );
-        if ( ( tc == CurrentContext::max_trial_count ) &&
-                ( !UnitTest::AreClose ( 1.0, actual ( CurrentContext::max_trial_count - 1 ), tolerance ) ) ) {
-                UnitTest::MemoryOutStream stream;
-                stream << CurrentContext::format_context ( __LINE__ ) << "expected " << mean << " +/- " << tolerance << " but was " << actual;
-
-                UnitTest::TestDetails details ( *UnitTest::CurrentTest::Details(), 0, false );
-
-                UnitTest::CurrentTest::Results()->OnTestFailure ( details, stream.GetText() );
-        }
+        double nc = normcoef ( hepdf.get(), xb, yb, nbins ( 0 ), nbins ( 1 ) );
+        CHECK_CLOSE(1.0,nc,0.01);
 }
 
@@ -168 +142 @@
         } while ( ( tc < CurrentContext::max_trial_count ) &&
                   !UnitTest::AreClose ( mean, actual ( tc - 1 ), delta ) );
+                          
         if ( ( tc == CurrentContext::max_trial_count ) &&
                 ( !UnitTest::AreClose ( mean, actual ( CurrentContext::max_trial_count - 1 ), delta ) ) ) {