diff options
| -rw-r--r-- | TODO.md | 9 | ||||
| -rw-r--r-- | src/serial-fractals.c | 2 | ||||
| -rw-r--r-- | src/shared-fractals.c | 61 |
3 files changed, 67 insertions, 5 deletions
@@ -3,13 +3,14 @@ ## Program * [ ] find way to render data into image files + * [x] read/write grids to/from file * [ ] write serial - * [ ] mandelbrot - * [ ] multibrot + * [x] mandelbrot + * [x] multibrot * [ ] julia * [ ] write shared - * [ ] mandelbrot - * [ ] multibrot + * [x] mandelbrot + * [x] multibrot * [ ] julia * [ ] write gpu * [ ] mandelbrot diff --git a/src/serial-fractals.c b/src/serial-fractals.c index bd4bc02..718fccb 100644 --- a/src/serial-fractals.c +++ b/src/serial-fractals.c @@ -1,6 +1,6 @@ -#include <complex.h> #include "fractals.h" #include "grids.h" + /* * Computes the number of iterations it takes for a point z0 to diverge * if the return value is equal to max_iterations, the point lies within the mandelbrot set diff --git a/src/shared-fractals.c b/src/shared-fractals.c index f0d837e..9f65e80 100644 --- a/src/shared-fractals.c +++ b/src/shared-fractals.c @@ -1 +1,62 @@ #include "fractals.h" + +#include <omp.h> +#include "fractals.h" + +/* + * Computes the number of iterations it takes for a point z0 to diverge + * if the return value is equal to max_iterations, the point lies within the mandelbrot set + * This should be identical to the version found in serial-fractals.c + */ +size_t mandelbrot(const long double complex z0, const size_t max_iterations){ + long double complex z = z0; + size_t iteration = 0; + while(cabsl(z) <= 2 && iteration < max_iterations){ + z = z*z + z0; + iteration++; + } + + return iteration; +} + +/* + * Fills a grid with mandelbrot values + */ +void mandelbrot_grid(grid_t* grid, const size_t max_iterations){ + const size_t size = grid->size; + size_t* data = grid->data; + + #pragma omp parallel for default(none) shared(data, size, grid, max_iterations) schedule(dynamic) + for(size_t i = 0; i < size; i++){ + data[i] = mandelbrot(grid_to_complex(grid, i), max_iterations); + } +} + +/* + * Computes the number of iterations it takes for a point z0 to diverge + * if the return value is equal to max_iterations, the point lies within the multibrot set + * This should be identical to the version found in serial-fractals.c + */ +size_t multibrot(const long double complex z0, const size_t max_iterations, const double d){ + long double complex z = z0; + size_t iteration = 0; + while(cabsl(z) <= 2 && iteration < max_iterations){ + z = cpowl(z, d) + z0; + iteration++; + } + return iteration; +} + + +/* + * Fills a grid with multibrot values + */ +void multibrot_grid(grid_t* grid, const size_t max_iterations, const double d){ + const size_t size = grid->size; + size_t* data = grid->data; + + #pragma omp parallel for default(none) shared(data, size, grid, max_iterations, d) schedule(dynamic) + for(size_t i = 0; i < size; i ++){ + data[i] = multibrot(grid_to_complex(grid, i), max_iterations, d); + } +} |