From 7f0f089f4dc98d8c6cc145f7355acf1a5ef9e53a Mon Sep 17 00:00:00 2001 From: David Phillips Date: Thu, 13 Aug 2015 10:37:31 +1200 Subject: Decomposed horrible loaded for loops, removed tricorn --- algorithms/mandelbrot.c | 13 ++++++++++--- 1 file changed, 10 insertions(+), 3 deletions(-) (limited to 'algorithms/mandelbrot.c') diff --git a/algorithms/mandelbrot.c b/algorithms/mandelbrot.c index 28bb188..8225d8c 100644 --- a/algorithms/mandelbrot.c +++ b/algorithms/mandelbrot.c @@ -37,12 +37,17 @@ void *generate_mandelbrot_section(void *section) double top = -1.75f; double left = -2.5f; - for (y = d->core, b = (d->core*(size_units/size)+top); y < size; b+=((cores*size_units)/size), y+=cores) + + /* FIXME document this */ + b = (d->core*(size_units/size)+top); + + for (y = d->core; y < size; y += cores) { - for (x = clust_id, a = (clust_id*(size_units/size)+left); x < size; a+=((clust_total*size_units)/size), x+=clust_total) + a = clust_id*(size_units/size)+left; /* FIXME document this */ + for (x = clust_id; x < size; x+=clust_total) { z = 0; - c = a+I*b; + c = a + I*b; for (i = 0; i < iterat; i++) { if (cabsf(z) >= 2) @@ -51,7 +56,9 @@ void *generate_mandelbrot_section(void *section) z = cpow(z , power) + c; } d->data[d->idx++] = (255*i)/iterat; + a += (clust_total*size_units)/size; } + b += (cores*size_units)/size; } return NULL; } -- cgit v1.1