[Teem-commit] SF.net SVN: teem:[5672] teem/trunk/Testing/air/amath.c

[<kin...@us...>]

Revision: 5672
          http://teem.svn.sourceforge.net/teem/?rev=5672&view=rev
Author:   kindlmann
Date:     2012-09-28 19:17:41 +0000 (Fri, 28 Sep 2012)
Log Message:
-----------
looking now at average instead of worst-case error with airCbrt

Modified Paths:
--------------
    teem/trunk/Testing/air/amath.c

Modified: teem/trunk/Testing/air/amath.c
===================================================================
--- teem/trunk/Testing/air/amath.c	2012-09-27 19:35:51 UTC (rev 5671)
+++ teem/trunk/Testing/air/amath.c	2012-09-28 19:17:41 UTC (rev 5672)
@@ -65,22 +65,25 @@
     double aa[2], uu, eps = 8e-27, error;
     rng = airRandMTStateNew(4242);
     airMopAdd(mop, rng, (airMopper)airRandMTStateNix, airMopAlways);
+    error = 0;
     for (ti=0; ti<testnum; ti++) {
+      double dif;
       airNormalRand_r(aa+0, aa+1, rng);
       for (pi=0; pi<2; pi++) {
         aa[pi] *= 10000;
         uu = airCbrt(aa[pi]);
-        error = fabs(aa[pi] - uu*uu*uu);
-        error = (error/aa[pi])*error;
-        if (error > eps) {
-          fprintf(stderr, "airCbrt(%.17g)^3 = %.17g^3 = %.17g != %.17g by "
-                  "(relatively) %g > eps %g \n",
-                  aa[pi], uu, uu*uu*uu, aa[pi],
-                  error, eps);
-          airMopError(mop); return 1;
-        }      
+        dif = fabs(aa[pi] - uu*uu*uu);
+        error += (dif/aa[pi])*dif;
       }
     }
+    error /= testnum;
+    if (error > eps) {
+      fprintf(stderr, "average cbrt error was %.17g > eps %0.17g\n",
+              error, eps);
+      airMopError(mop); return 1;
+    } else {
+      fprintf(stderr, "average cbrt error = %.17g\n", error);
+    }
   }
 
   exit(0);

This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.