First upload

Author: EStog

LICENCE

@@ -0,0 +1,8 @@
+# Mandelbrot generation experiments with MPI and OpenMP
+
+This project include:
+- C++ programs for the generation of the Mandelbrot set using MPI and OpenMP, as well as a sequential version.
+- A test running Python program to make performance testing to the C++ programs.
+
+Author: Ernesto Soto Gómez.
+<esto.yinyang@gmail.com>, <esoto@uci.cu>

README.md

@@ -0,0 +1,29 @@
+#include "compute_mandelbrot_subset.h"
+#include <complex>
+
+using namespace std;
+
+const double x_begin = -2.0, x_end = 0.6;
+const double y_begin = -1.2, y_end = 1.2;
+const double x_step = (x_end-x_begin) / (x_resolution-1);
+const double y_step = (y_end-y_begin) / (y_resolution-1);
+
+void compute_mandelbrot_subset(int* result, int iter_limit, int start, int end) {
+    int i, j;
+    complex<double> c, z;
+
+    #pragma omp parallel shared(result, iter_limit, start, end) private(i, j, c, z)
+    #pragma omp for schedule(runtime)
+    for (i = start; i < end; i++) {
+        c = complex<double>(
+                x_begin + (i % x_resolution) * x_step,
+                y_begin + (i / x_resolution) * y_step
+        );
+        z = 0; j = 0;
+        while (norm(z) <= 4 && j < iter_limit) {
+            z = z*z + c;
+            j++;
+        }
+        result[i-start] = j;
+    }
+}

code/common/compute_mandelbrot_subset.cpp

@@ -0,0 +1,8 @@
+#pragma once
+
+const int x_resolution = 1024;
+const int y_resolution = 1024;
+const int result_size = x_resolution * y_resolution;
+
+void compute_mandelbrot_subset(int* result, int iter_limit,
+                               int start = 0, int end = result_size);

code/common/compute_mandelbrot_subset.h

@@ -0,0 +1,36 @@
+#include "misc.h"
+#include <iostream>
+
+using namespace std;
+
+//process command arguments
+void get_iter_limit(int argc, char **argv, int &iter_limit, bool &print) {
+    try {
+        switch (argc) {
+            case 2:
+                iter_limit = stoi(argv[--argc]);
+                print = false;
+                return;
+            case 3:
+                iter_limit = stoi(argv[--argc]);
+                print = stoi(argv[--argc]);
+                return;
+            case 1:
+                cout << endl;
+                cout << "Pass at least one argument" << endl;
+        }
+    }
+    catch (exception) {
+        cout << endl;
+        cout << "Wrong arguments" << endl;
+    }
+    cout << "Parameters syntax after program name: [<print>] <iter_limit>" << endl;
+    cout << "<print>: Whether to print the set to a file or not. This must be somethig that can be converted to boolean" << endl;
+    cout << "<iter_limit>: The iteration limit. This must be a positive integer" << endl;
+    cout << endl;
+    exit(1);
+}
+
+void output_execution_time(double time) {
+    cout << time << endl;
+}

code/common/misc.cpp

@@ -0,0 +1,6 @@
+#pragma once
+
+
+void get_iter_limit(int argc, char **argv, int &iter_limit, bool &print);
+
+void output_execution_time(double time);

code/common/misc.h

@@ -0,0 +1,10 @@
+#include "now.h"
+#include <ctime>
+#include <cstdlib>
+
+double now() {
+    struct timespec tp;
+    if (clock_gettime(CLOCK_MONOTONIC_RAW, &tp) != 0)
+        exit(1);
+    return tp.tv_sec + tp.tv_nsec / (double)1000000000;
+}

code/common/now.cpp

@@ -0,0 +1,4 @@
+#pragma once
+
+
+double now();

code/common/now.h

@@ -0,0 +1,42 @@
+#include "print_result.h"
+#include "compute_mandelbrot_subset.h"
+#include <fstream>
+#include <iostream>
+
+using namespace std;
+
+const unsigned short max_color = 65535;
+const string mandelbrot_file_name ="../../../images/mandelbrot";
+
+// Used to store mandelbrot set into an image file.
+void print_result(int* result, int iter_limit, const char* name) {
+    ofstream file;
+    file.open(mandelbrot_file_name);
+    file << "P3" << endl;
+    file << x_resolution << " " << y_resolution << endl;
+    file << max_color << endl;
+    for (int i = 1; i <= y_resolution; i++) {
+        for (int j = 0; j < x_resolution; j++) {
+            int k = result[(i - 1) * x_resolution + j];
+            file << " ";
+            if (k == iter_limit)
+                file << 0 << " " << 0 << " " << 0;
+            else {
+                unsigned short x = k * max_color / iter_limit;
+                file << x / 4 << " ";
+                file << x / 2 + max_color / 2 << " ";
+                file << x;
+            }
+            file << " ";
+        }
+        file << endl;
+    }
+    file.close();
+
+    string command = "pnmtojpeg -quality=100 -smooth=100 -optimize ./"+mandelbrot_file_name+" > "+mandelbrot_file_name+"_limit="+to_string(iter_limit)+"_"+name+".jpeg";
+
+    if (system(command.c_str()) != 0)
+        exit(1);
+
+    system(("rm "+mandelbrot_file_name).c_str());
+}

code/common/print_result.cpp

@@ -0,0 +1,3 @@
+#pragma once
+
+void print_result(int* result, int iter_limit, const char* name);