fixed regression with perturbation rescaling

Author: jackyarndley

default.toml

@@ -1,3 +1,3 @@
 image_width = 1000
 image_height = 1000
-frames = 1
+frames = 99999

src/math/perturbation.rs

@@ -60,6 +60,8 @@ impl Perturbation {
 
                                 // check the counter, if it is > 250, do a normalisation
                                 if scaled_iterations > 250 {
+                                    pixel.delta_current.reduce();
+
                                     scaled_scale_factor_1 = 1.0f64.ldexp(pixel.delta_current.exponent);
                                     scaled_delta_reference = 1.0f64.ldexp(pixel.delta_reference.exponent - pixel.delta_current.exponent) * pixel.delta_reference.mantissa;
 

src/math/reference.rs

@@ -100,6 +100,7 @@ impl Reference {
             )
         }
 
+        // only needed on central reference
         self.approximation_data.push(z_extended);
 
         if self.start_iteration == 1 {

src/math/series_approximation.rs

@@ -14,7 +14,7 @@ pub struct SeriesApproximation {
     approximation_probes: Vec<Vec<ComplexExtended>>,
     approximation_probes_derivative: Vec<Vec<ComplexExtended>>,
     pub delta_top_left: ComplexExtended,
-    valid_coefficients: Vec<ComplexExtended>,
+    pub valid_coefficients: Vec<ComplexExtended>,
     pub valid_iteration: usize,
 }
 
@@ -57,10 +57,8 @@ impl SeriesApproximation {
             // Calculate the new coefficents
             for k in 2..=self.order {
                 let mut sum = coefficients[0] * coefficients[k];
-                sum.reduce();
 
                 for j in 1..=((k - 1) / 2) {
-                    sum.reduce();
                     sum += coefficients[j] * coefficients[k - j];
                 }
                 sum *= 2.0;

src/renderer.rs

@@ -70,7 +70,7 @@ impl FractalRenderer {
             zoom,
             maximum_iteration,
             glitch_tolerance,
-            data_export: DataExport::new(image_width, image_height, display_glitches, DataType::RAW),
+            data_export: DataExport::new(image_width, image_height, display_glitches, DataType::COLOUR),
             start_render_time: Instant::now(),
             remaining_frames,
             zoom_scale_factor,
@@ -112,6 +112,7 @@ impl FractalRenderer {
                 let j = index / self.image_width;
                 let element = ComplexFixed::new(i as f64 * delta_pixel + delta_top_left.re, j as f64 * delta_pixel + delta_top_left.im);
                 let point_delta = ComplexExtended::new(element, -self.zoom.exponent);
+
                 let new_delta = self.series_approximation.evaluate(point_delta, None);
 
                 PixelData {