fixed the series approximation

Author: jackyarndley

output.png

@@ -61,7 +61,7 @@ fn main() {
         center_im,
         0.01,
         false,
-        16
+        64
     );
 
     let time = Instant::now();

src/bin/main.rs

@@ -67,6 +67,8 @@ impl SeriesApproximation {
 
             self.next_coefficients[0] = to_extended(&self.z);
             self.next_coefficients[1] = self.coefficients[0] * self.coefficients[1] * 2.0 + add_value;
+            self.next_coefficients[0].reduce();
+            self.next_coefficients[1].reduce();
 
             // Calculate the new coefficents
             for k in 2..=self.order {
@@ -109,7 +111,8 @@ impl SeriesApproximation {
 
                 // Check to make sure that the derivative is greater than or equal to 1
                 if derivative.to_float() < 1.0 {
-                    derivative = FloatExtended::new(1.0, 0);
+                    derivative.mantissa = 1.0;
+                    derivative.exponent = 0;
                 }
 
                 // Check that the error over the derivative is less than the pixel spacing
@@ -189,6 +192,21 @@ impl SeriesApproximation {
         approximation
     }
 
+    fn evaluate_next(&self, point_delta: ComplexExtended) -> ComplexExtended {
+        // 1907 ms packing opus 4K
+        // Horner's rule
+        let mut approximation = self.next_coefficients[self.order];
+
+        for k in (1..=(self.order - 1)).rev() {
+            approximation *= point_delta;
+            approximation += self.coefficients[k];
+        }
+
+        approximation *= point_delta;
+        approximation.reduce();
+        approximation
+    }
+
     // pub fn evaluate_derivative(&self, point_delta: ComplexExtended) -> FloatExtended {
     //     let mut original_point_derivative_n = ComplexExtended::new(1.0, 0, 0.0, 0);
     //     let mut approximation_derivative = ComplexExtended::new(0.0, 0, 0.0, 0);

src/math/series_approximation.rs

@@ -11,25 +11,20 @@ pub struct ComplexExtended {
 }
 
 impl ComplexExtended {
-    // maybe set up different versions for these function if the reduce is needed
     #[inline]
     pub fn new(mantissa: Complex<f64>, exponent: i32) -> Self {
-        let mut temp = ComplexExtended {
+        ComplexExtended {
             mantissa,
             exponent
-        };
-        temp.reduce();
-        temp
+        }
     }
 
     #[inline]
     pub fn new2(re: f64, im: f64, exponent: i32) -> Self {
-        let mut temp = ComplexExtended {
+        ComplexExtended {
             mantissa: Complex::<f64>::new(re, im),
             exponent
-        };
-        temp.reduce();
-        temp
+        }
     }
 
     #[inline]