Changeset 1011Tweet

Timestamp:

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

Author:

sam

Message:

test: more Remez exchange experimentations.

Location:

trunk

Files:

• .gitignore (modified) (1 diff)
• test/math/Makefile.am (modified) (3 diffs)
• test/math/poly.cpp (added)
• test/math/remez.cpp (modified) (2 diffs)

trunk/.gitignore

@@ -61 +61 @@
 test/sandbox
 test/math/pi
+test/math/poly
 test/math/remez
 tools/make-font

trunk/test/math/Makefile.am

@@ -4 +4 @@
 all-local: $(noinst_PROGRAMS)
         test x$(MAKE_FSELF) = xno || make_fself pi$(EXEEXT) pi.self
+        test x$(MAKE_FSELF) = xno || make_fself poly$(EXEEXT) pi.self
         test x$(MAKE_FSELF) = xno || make_fself remez$(EXEEXT) remez.self
 
@@ -10 +11 @@
              $(noinst_PROGRAMS:%$(EXEEXT)=%.exe)
 
-noinst_PROGRAMS = pi remez
+noinst_PROGRAMS = pi poly remez
 
 pi_SOURCES = pi.cpp
@@ -17 +18 @@
 pi_DEPENDENCIES = $(top_builddir)/src/liblol.a
 
+poly_SOURCES = poly.cpp
+poly_CPPFLAGS = @LOL_CFLAGS@ @PIPI_CFLAGS@
+poly_LDFLAGS = $(top_builddir)/src/liblol.a @LOL_LIBS@ @PIPI_LIBS@
+poly_DEPENDENCIES = $(top_builddir)/src/liblol.a
+
 remez_SOURCES = remez.cpp remez-matrix.h remez-solver.h
 remez_CPPFLAGS = @LOL_CFLAGS@ @PIPI_CFLAGS@

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);