fixed the glitch correction

Author: jackyarndley

output.png

@@ -59,9 +59,9 @@ fn main() {
         iterations.parse::<usize>().unwrap(),
         center_re,
         center_im,
-        0.1,
-        true,
-        0
+        0.01,
+        false,
+        16
     );
 
     let time = Instant::now();

src/bin/main.rs

@@ -130,7 +130,7 @@ impl FractalRenderer {
         println!("{:<14}{:>6} ms", "Iteration", time.elapsed().as_millis());
 
         let time = Instant::now();
-        Colouring::Iteration.run(&pixel_data, &mut self.image, self.maximum_iteration, delta_pixel, &reference);
+        Colouring::IterationSmooth.run(&pixel_data, &mut self.image, self.maximum_iteration, delta_pixel, &reference);
         println!("{:<14}{:>6} ms", "Coloring", time.elapsed().as_millis());
 
         let time = Instant::now();
@@ -166,7 +166,7 @@ impl FractalRenderer {
 
             Perturbation::iterate(&mut pixel_data, &r, r.current_iteration);
 
-            Colouring::Iteration.run(&pixel_data, &mut self.image, self.maximum_iteration, delta_pixel, &r);
+            Colouring::IterationSmooth.run(&pixel_data, &mut self.image, self.maximum_iteration, delta_pixel, &r);
 
             // Remove all non-glitched points from the remaining points
             pixel_data.retain(|packet| {

src/renderer.rs

@@ -132,7 +132,7 @@ impl Sub<ComplexExtended> for ComplexExtended {
         if self.exponent > other.exponent {
             ComplexExtended::new(self.mantissa - other.mantissa * 1.0f64.ldexp(other.exponent - self.exponent), self.exponent)
         } else {
-            ComplexExtended::new(other.mantissa - self.mantissa * 1.0f64.ldexp(self.exponent - other.exponent), other.exponent)
+            ComplexExtended::new(self.mantissa * 1.0f64.ldexp(self.exponent - other.exponent) - other.mantissa, other.exponent)
         }
     }
 }