Optimization and norming a bitHEADmaster

1 files changed, 14 insertions, 18 deletions

diff --git a/src/fractals/tricorn.c b/src/fractals/tricorn.c
index 052cc2f..2abe01b 100644
--- a/src/fractals/tricorn.c
+++ b/src/fractals/tricorn.c
@@ -6,37 +6,33 @@
 /*   By: cacharle <marvin@42.fr>                    +#+  +:+       +#+        */
 /*                                                +#+#+#+#+#+   +#+           */
 /*   Created: 2020/02/24 16:00:46 by cacharle          #+#    #+#             */
-/*   Updated: 2020/02/25 07:34:01 by cacharle         ###   ########.fr       */
+/*   Updated: 2020/02/27 12:24:35 by cacharle         ###   ########.fr       */
 /*                                                                            */
 /* ************************************************************************** */
 
 #include "fractol.h"
 
-#define TRICORN_MAX_ITERATION 20
-#define TRICORN_ESCAPE_RADIUS_SQUARED 100
+#define TRICORN_ESCAPE_RADIUS_SQUARED 4.0
 
-int	tricorn(t_state *state, t_complex z)
+int	tricorn(t_state *state, t_complex c)
 {
-	int		n;
-	double	zr;
-	double	zi;
-	double	zr_square;
-	double	zi_square;
-	double	tmp;
+	int			n;
+	t_complex	z;
+	t_complex	z_square;
+	double		tmp;
 	
 	(void)state;
-	zr = z.r;
-	zi = z.i;
+	z = c;
 	n = -1;
 	while (++n < state->iterations)
 	{
-		zi_square = zi * zi;
-		zr_square = zr * zr;
-		if (zr_square + zi_square > TRICORN_ESCAPE_RADIUS_SQUARED)
+		z_square.i = z.i * z.i;
+		z_square.r = z.r * z.r;
+		if (z_square.r + z_square.i > TRICORN_ESCAPE_RADIUS_SQUARED)
 			break;
-		tmp = zr_square - zi_square + z.r;
-		zi = -2.0 * zr * zi + z.i;
-		zr = tmp;
+		tmp = z_square.r - z_square.i + z.r;
+		z.i = -2.0 * z.r * z.i + c.i;
+		z.r = tmp;
 	}
 	return (n);
 }