From c4b5e0a48835b6bdbf4d3d92fd42b29a22662ed4 Mon Sep 17 00:00:00 2001 From: Leonard Kugis Date: Tue, 16 Jan 2018 01:21:27 +0100 Subject: Mandelbrot drawing possible --- src/sets.c | 46 ++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 46 insertions(+) create mode 100644 src/sets.c (limited to 'src/sets.c') diff --git a/src/sets.c b/src/sets.c new file mode 100644 index 0000000..6787a6e --- /dev/null +++ b/src/sets.c @@ -0,0 +1,46 @@ +/* + * sets.c + * + * Created on: 15.01.2018 + * Author: Superleo1810 + */ + +#include "sets.h" + +u32 mandelbrot_s(double x, double y, u32 iterations) +{ + double cx = x, cy = y, x2; + u32 m = 0; + while(m <= iterations && (x*x)+(y*y) <= 4) + { + x2 = x; + x = (x*x) - (y*y) + cx; + y = 2.0*x2*y + cy; + m++; + } + return m; +} + +u32 mandelbrot_r(double x, double y, u32 iterations) +{ + return _mandelbrot_r(x, y, 0.0, 0.0, 0, iterations, 4.0); +} + +u32 _mandelbrot_r(double x, double y, double zx, double zy, u32 n, u32 iterations, double threshold) +{ + if ((n <= iterations) && ((zx * zx + zy * zy) < threshold)) { + double zx_new = (zx * zx - zy * zy + x); + double zy_new = (2 * zx * zy + y); + if ((zx_new == zx) && (zy_new == zy)) { + return iterations; + } + n = _mandelbrot_r(x, y, zx_new, zy_new, n + 1, iterations, threshold); + } + return n; +} + +u32 julia(double x, double y, u32 iterations) +{ + // TODO: Julia-Menge + return 0; +} -- cgit v1.2.1