OMP Parallel

OMP functional parallel here

Author: pouyalj

.vscode/c_cpp_properties.json

@@ -0,0 +1,21 @@
+{
+    "configurations": [
+        {
+            "name": "Win32",
+            "includePath": [
+                "${workspaceFolder}/**"
+            ],
+            "defines": [
+                "_DEBUG",
+                "UNICODE",
+                "_UNICODE"
+            ],
+            "compilerPath": "/usr/bin/gcc",
+            "cStandard": "gnu11",
+            "cppStandard": "gnu++14",
+            "intelliSenseMode": "linux-gcc-x64",
+            "compileCommands": "${workspaceFolder}/compile_commands.json"
+        }
+    ],
+    "version": 4
+}

.vscode/settings.json

@@ -0,0 +1,5 @@
+{
+    "files.associations": {
+        "complex": "cpp"
+    }
+}

CMakeLists.txt

@@ -4,9 +4,8 @@ cmake_minimum_required(VERSION 3.2)
 set(CMAKE_CXX_STANDARD 17)
 set(CMAKE_CXX_STANDARD_REQUIRED ON)
 set(CMAKE_CXX_EXTENSIONS OFF)
-SET(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -std=c++11 -pthread")
-
-project(fractals)
+SET(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -std=c++17 -pthread")
+SET(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -std=c++17 -fopenmp")
 
 # Add dependencies
 find_package(Qt5 COMPONENTS Widgets PrintSupport REQUIRED)

exported/qcustomplot.cmake

@@ -4,7 +4,7 @@ if("${CMAKE_MAJOR_VERSION}.${CMAKE_MINOR_VERSION}" LESS 2.5)
    message(FATAL_ERROR "CMake >= 2.6.0 required")
 endif()
 cmake_policy(PUSH)
-cmake_policy(VERSION 2.6)
+cmake_policy(VERSION 2.6...3.18)
 #----------------------------------------------------------------
 # Generated CMake target import file.
 #----------------------------------------------------------------

include/fractals.hpp

@@ -3,6 +3,7 @@
 #include <iostream>
 #include <vector>
 #include <thread>
+#include <omp.h>
 
 namespace fractals {
 
@@ -27,37 +28,36 @@ ComplexDouble mandelbrot_process(ComplexDouble c, int boundary = 2,
 }
 
 void mandelbrot(double start_val, double end_val, double step_val) {
+    FractalSet temp_set;
     double step = step_val; // 0.01;
     double start = start_val; //-2.0;
     double end = end_val; //2.0;
-    // FractalSet m_set;
     for (double i = start; i < end; i += step) {
         for (double j = start; j < end; j += step) {
             ComplexDouble c(i, j);
-            m_set.push_back(mandelbrot_process(c));
+            temp_set.push_back(mandelbrot_process(c));
             }
     }
+    #pragma omp critical
+    m_set.insert(m_set.end(), temp_set.begin(), temp_set.end());
 
 
     // return m_set;
 }  // namespace fractals
 
 
 FractalSet thread(){
-    // m_set.push_back(0);
-    // Pre loop thread setup
-    int num_threads = std::thread::hardware_concurrency(); 
-    std::vector<std::thread> threads;
-    // std::mutex critical;
+    // // Pre loop thread setup
+    int num_threads = omp_get_max_threads();
     std::cout << num_threads << std::endl;
+    omp_set_num_threads(omp_get_max_threads());
+    #pragma omp parallel
+    #pragma omp for
     for(int t=0; t<num_threads; t++) {
         double start_here = -2.0+(t*4/num_threads);
         double end_here = -2.0+((t+1)*4/num_threads);
-        threads.push_back(std::thread(mandelbrot, start_here, end_here, 0.01));
-        threads[t].join();
+        mandelbrot(start_here, end_here, 0.001);
     }
-    // std::thread first(mandelbrot, -2.0, 2.0, 0.01, m_set);
-    // first.join();
 
     return m_set;
 }

src/mandelbrot.cpp

@@ -1,7 +1,8 @@
 #include <iostream>
-#include <thread>
 #include <fstream>
 #include <complex>
+#include <chrono>
+
 
 
 #include "../external/qcustomplot/qcustomplot.h"
@@ -25,7 +26,8 @@ using ComplexDouble = std::complex<double>;
 using FractalSet = std::vector<ComplexDouble>;
 
 int main() {
-    FractalSet m_set = fractals::thread();
+    auto start = std::chrono::high_resolution_clock::now();
+    auto m_set = fractals::thread();
     // std::cout << "{ ";
     // for (auto i : m_set) {
     //     std::cout << i << " ";
@@ -44,6 +46,9 @@ int main() {
     // custom_plot.addGraph();
     // custom_plot.graph(0)->setPen(QPen(Qt::blue)); // line color blue for
     // first graph custom_plot.graph(0)->setBrush(QBrush(QColor(0, 0, 255,20))); // first graph will be filled with translucent blue
+    auto stop = std::chrono::high_resolution_clock::now();
+    auto duration = std::chrono::duration_cast<std::chrono::milliseconds>(stop - start);
+    std::cout << "Runtime is " << duration.count() << " milliseconds" <<std::endl;
     return 0;
 }