Changeset 737 for library/tests/epdf_harness.cpp

Timestamp:

11/25/09 12:14:38 (15 years ago)

Author:

mido

Message:

ASTYLER RUN OVER THE WHOLE LIBRARY, JUPEE

Files:

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

library/tests/epdf_harness.cpp

@@ -45 +45 @@
         if ( support ) { // support is given
                 grid_fnc ep_disc;
-                ep_disc.set_support(*support);
-                ep_disc.set_values(*hepdf);
+                ep_disc.set_support ( *support );
+                ep_disc.set_values ( *hepdf );
                 // ep_disc is discretized at support points
 
-                double point_volume =prod(support->_steps());
-                CHECK_CLOSE(1.0, sum(ep_disc._values())*point_volume, 0.01);
-
-                vec pdf=ep_disc._values();
-                pdf /=sum(pdf); // normalize
-                
-                vec mea=pdf(0) * support->first_vec();
-                mat Remp=pdf(0) * outer_product(support->act_vec(), support->act_vec());
+                double point_volume = prod ( support->_steps() );
+                CHECK_CLOSE ( 1.0, sum ( ep_disc._values() ) *point_volume, 0.01 );
+
+                vec pdf = ep_disc._values();
+                pdf /= sum ( pdf ); // normalize
+
+                vec mea = pdf ( 0 ) * support->first_vec();
+                mat Remp = pdf ( 0 ) * outer_product ( support->act_vec(), support->act_vec() );
 
                 // run through all points
-                for (int i=1; i<support->points(); i++){
-                        mea += pdf(i)*support->next_vec();
-                        Remp += pdf(i) * outer_product(support->act_vec(), support->act_vec());
-                }
-                CHECK_CLOSE(mean, mea, tolerance);
-                CHECK_CLOSE(R, Remp-outer_product(mea,mea), tolerance);
+                for ( int i = 1; i < support->points(); i++ ) {
+                        mea += pdf ( i ) * support->next_vec();
+                        Remp += pdf ( i ) * outer_product ( support->act_vec(), support->act_vec() );
+                }
+                CHECK_CLOSE ( mean, mea, tolerance );
+                CHECK_CLOSE ( R, Remp - outer_product ( mea, mea ), tolerance );
         }
 
@@ -70 +70 @@
                 check_sample_mean();
         }
-        if (R.rows()>0){
+        if ( R.rows() > 0 ) {
                 check_covariance();
         }
@@ -133 +133 @@
         } 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 ) ) ) {