/*
 * 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;
}