Changeset 1011 for trunk/test/math/remez.cppTweet

Timestamp:

Oct 6, 2011, 9:18:17 AM (10 years ago)

Author:

sam

Message:

test: more Remez exchange experimentations.

File:

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

trunk/test/math/remez.cpp

@@ -27 +27 @@
 static real myfun(real const &x)
 {
+    static real const a0 = real::R_1;
+    static real const a1 = real(-11184811) >> 26;
+    static real const b1 = real(-1) / real(6);
+    static real const b2 = real(1) / real(120);
     real y = sqrt(x);
     if (!y)
-        return real::R_PI_2;
-    return sin(real::R_PI_2 * y) / y;
+        return b2;
+    return ((sin(y) / y - a0) / x - a1) / x;
 }
 
@@ -40 +44 @@
 int main(int argc, char **argv)
 {
-    RemezSolver<4> solver;
-    solver.Run(0, 1, myfun, myfun, 15);
+    RemezSolver<3> solver;
+    //solver.Run(0, 1, myfun, myfun, 15);
+    solver.Run(0, real::R_PI * real::R_PI >> 2, myfun, myfun, 15);
     //solver.Run(-1, 1, myfun, myfun, 15);
     //solver.Run(0, real::R_PI * real::R_PI >> 4, myfun, myfun, 15);