Changeset 1005Tweet

Timestamp:

Oct 4, 2011, 3:35:45 PM (11 years ago)

Author:

sam

Message:

test: add support for relative error in the Remez exchange test program.

File:

• trunk/test/math/remez.cpp (modified) (5 diffs)

trunk/test/math/remez.cpp

@@ -14 +14 @@
 
 #include <cstring>
+#include <cstdio>
 
 #include "core.h"
@@ -20 +21 @@
 
 /* The order of the approximation we're looking for */
-static int const ORDER = 8;
+static int const ORDER = 4;
 
 /* The function we want to approximate */
 static real myfun(real const &x)
 {
-    static real const one = 1.0;
-    if (!x)
-        return real::R_PI_2;
-    return sin(x * real::R_PI_2) / x;
-    //return cos(x) - one;
-    //return exp(x);
+    return exp(x);
+    //if (!x)
+    //    return real::R_PI_2;
+    //return sin(x * real::R_PI_2) / x;
+}
+
+static real myerror(real const &x)
+{
+    return myfun(x);
+    //return real::R_1;
 }
 
@@ -206 +211 @@
 static void remez_findzeroes(real *coeff, real *zeroes, real *control)
 {
-    /* FIXME: this is fake for now */
     for (int i = 0; i < ORDER + 1; i++)
     {
@@ -278 +282 @@
                     final = maxerror;
                 control[i] = (a + b) * (real)0.5;
-                printf("%g (in %g)\n", (double)maxerror, (double)control[i]);
+                printf("%g (at %g)\n", (double)maxerror, (double)control[i]);
                 break;
             }
@@ -314 +318 @@
             mat.m[i][n] = sum;
         }
-        mat.m[i][ORDER + 1] = (real)(-1 + (i & 1) * 2);
+        if (i & 1)
+            mat.m[i][ORDER + 1] = fabs(myerror(control[i]));
+        else
+            mat.m[i][ORDER + 1] = -fabs(myerror(control[i]));
     }