Tengo un código en C++ que intento paralelizar con OpenMPI, consiste en buscar los números primos hasta un tope indicado:
#include <vector>
#include <cmath>
#include <iostream>
#include <omp.h>
#include <ctime>
using namespace std;
int main() {
unsigned t0, t1;
t0 = clock();
omp_set_num_threads(4);
auto maximo = 1000000ul; //numero maximo a buscar
// Criba secuencial
std::vector<int> primos{ 2 }; // Añadimos el 2 al array para poder saltarnos los pares
primos.reserve(maximo);
#pragma omp parallel for
for (auto posible_primo = 3ul; posible_primo < maximo; posible_primo += 2) {
/* Comprobamos si es primo buscando si en el array hay algun
numero primo que lo divida, pero solo hasta raiz de n. */
const auto limite = sqrt(posible_primo);
bool es_divisible = false;
for (auto divisor : primos) {
if (limite < divisor)
break;
if (es_divisible = !(posible_primo % divisor))
break;
}
if (!es_divisible) {
primos.push_back(posible_primo);
}
}
for (auto primo : primos)
std::cout << primo << ' ';
t1 = clock();
double time = (double(t1 - t0) / CLOCKS_PER_SEC);
cout << "\n\n\tExecution Time: " << time << endl;
return 0;
}
Tengo que paralelizar el bucle exterior (el que recorre los números desde 3 hasta el límite superior), para ello he utilizado #pragma omp parallel for
Cuando ejecuto el algoritmo, tarda lo mismo que si quito la sentencia paralela.
¿Qué hago mal? tampoco se cómo hacer para evitar las condiciones de carrera en la escritura del número primo en el vector.
std::vector<int> primos
y se bloquean entre ellos.