From 8db6918f320b792698f385bb2ea58ddd94cba6fe Mon Sep 17 00:00:00 2001 From: Charles Date: Mon, 24 Feb 2020 13:48:17 +0100 Subject: window_complex position buffer, basic palette, mandelbrot working --- src/fractals/mandelbrot.c | 41 +++++++++++++++++++++++++++++++++++++++++ 1 file changed, 41 insertions(+) (limited to 'src/fractals') diff --git a/src/fractals/mandelbrot.c b/src/fractals/mandelbrot.c index e69de29..4211cfa 100644 --- a/src/fractals/mandelbrot.c +++ b/src/fractals/mandelbrot.c @@ -0,0 +1,41 @@ +/* ************************************************************************** */ +/* */ +/* ::: :::::::: */ +/* mandelbrot.c :+: :+: :+: */ +/* +:+ +:+ +:+ */ +/* By: cacharle +#+ +:+ +#+ */ +/* +#+#+#+#+#+ +#+ */ +/* Created: 2020/02/24 11:07:41 by cacharle #+# #+# */ +/* Updated: 2020/02/24 13:46:22 by cacharle ### ########.fr */ +/* */ +/* ************************************************************************** */ + +#include "fractol.h" + +#define MANDEL_MAX_ITERATION 20 +#define MANDEL_ESCAPE_RADIUS_SQUARED 100 + +int mandelbrot(t_complex z) +{ + int n; + double zr; + double zi; + double zr_square; + double zi_square; + + zr = z.r; + zi = z.i; + n = -1; + while (++n < MANDEL_MAX_ITERATION) + { + zi_square = zi * zi; + zr_square = zr * zr; + if (zr_square + zi_square > MANDEL_ESCAPE_RADIUS_SQUARED) + return (n); + zi = 2.0 * zr * zi; + zr = zr_square - zi_square; + zi += z.i; + zr += z.r; + } + return (n); +} -- cgit