Changeset 1027Tweet

Timestamp:

Oct 17, 2011, 1:45:40 PM (11 years ago)

Author:

sam

Message:

test: minor Remez algorithm tests.

Location:

trunk/test/math

Files:

• remez-solver.h (modified) (3 diffs)
• remez.cpp (modified) (1 diff)

trunk/test/math/remez-solver.h

@@ -155 +155 @@
                 b = zeroes[i];
 
-            printf("Error for [%g..%g]: ", (double)a, (double)b);
             for (;;)
             {
@@ -184 +183 @@
                         final = e;
                     control[i] = (a + b) >> 1;
-                    printf("%g (at %g)\n", (double)e, (double)control[i]);
                     break;
                 }
@@ -291 +289 @@
         }
 
-        printf("Final polynomial:\n");
+        printf("Polynomial estimate:\n");
         for (int j = 0; j < ORDER + 1; j++)
         {

trunk/test/math/remez.cpp

@@ -25 +25 @@
 
 /* The function we want to approximate */
-static real myfun(real const &x)
+real myfun(real const &y)
 {
-    real y = sqrt(x);
-    return (sin(y) - y) / (x * y);
+    real x = sqrt(y);
+    return (sin(x) - x) / (x * y);
 }
 
-static real myerr(real const &x)
+real myerr(real const &y)
 {
-    real y = sqrt(x);
-    return re(x * y);
+    real x = sqrt(y);
+    return sin(x) / (x * y);
 }
 
 int main(void)
 {
-    RemezSolver<4> solver;
-    solver.Run(real::R_1 >> 400, real::R_PI_2 * real::R_PI_2, myfun, myerr, 15);
+    RemezSolver<6> solver;
+    solver.Run(real::R_1 >> 400, real::R_PI_2 * real::R_PI_2, myfun, myerr, 40);
 
     return EXIT_SUCCESS;