From dab9efb7b745fe884fb72017591dce95978c19c4 Mon Sep 17 00:00:00 2001 From: Charles Date: Thu, 27 Feb 2020 14:47:23 +0100 Subject: Optimization and norming a bit --- src/fractals/tricorn.c | 32 ++++++++++++++------------------ 1 file changed, 14 insertions(+), 18 deletions(-) (limited to 'src/fractals/tricorn.c') 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 +#+ +:+ +#+ */ /* +#+#+#+#+#+ +#+ */ /* 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); } -- cgit