Changes between Version 5 and Version 6 of doc/maths/remez/tutorial-fixing-parameters Tweet

Timestamp:

Jan 9, 2012, 11:38:19 PM (11 years ago)

Author:

sam

Comment:

fix tutorial for upcoming LolRemez 0.2

doc/maths/remez/tutorial-fixing-parameters

@@ -13 +13 @@
 {{{
 #!cpp
-const double a1 = 9.999999766327004083636926544075452569800e-1;
+const double a1 = 9.999999765898820673279342160490060830302e-1;
 }}}
 
@@ -68 +68 @@
 {
     RemezSolver<3, real> solver;
-    solver.Run("1e-1000", real::R_PI_2 * real::R_PI_2, f, g, 30);
+    solver.Run("1e-1000", real::R_PI_2 * real::R_PI_2, f, g, 40);
     return 0;
 }
@@ -87 +87 @@
 
 {{{
-Final error: 4.613764035583927492630791542066631050746e-9
+Step 8 error: 4.618689007546850899022101933442449327546e-9
 Polynomial estimate:
-x**0*-1.666665708792613295297026208726412532085e-1
-+x**1*8.333017134192478312472752663154642556843e-3
-+x**2*-1.980660618990173078599247774582918018345e-4
-+x**3*2.600038299369332084157603784882057661080e-6
+x**0*-1.666665709650470145824129400050267289858e-1
++x**1*8.333017291562218127986291618761571373087e-3
++x**2*-1.980661520135080504411629636078917643846e-4
++x**3*2.600054767890361277123254766503271638682e-6
 }}}
 
@@ -101 +101 @@
 double fastsin2(double x)
 {
-    const double a3 = -1.666665708792613295297026208726412532085e-1;
-    const double a5 = 8.333017134192478312472752663154642556843e-3;
-    const double a7 = -1.980660618990173078599247774582918018345e-4;
-    const double a9 = 2.600038299369332084157603784882057661080e-6;
+    const double a3 = -1.666665709650470145824129400050267289858e-1;
+    const double a5 = 8.333017291562218127986291618761571373087e-3;
+    const double a7 = -1.980661520135080504411629636078917643846e-4;
+    const double a9 = 2.600054767890361277123254766503271638682e-6;
 
     return x + x*x*x * (a3 + x*x * (a5 + x*x * (a7 + x*x * a9))));