indivisible-legacy/src/Main.cpp

53 lines
868 B
C++
Raw Normal View History

2016-11-30 13:32:06 +00:00
#include <iostream>
#include <omp.h>
2016-11-30 14:15:52 +00:00
#include <vector>
#include <csignal>
#include <cassert>
2016-11-30 14:15:52 +00:00
static bool run;
2016-11-30 13:32:06 +00:00
int main(void) {
std::cout << "Indivisible v0.1\n";
2016-11-30 14:15:52 +00:00
run = true;
2016-12-01 11:35:15 +00:00
signal(SIGINT, [](int signum) {
std::cout << "Exiting (" << signum << ")\n";
run = false;
});
2016-11-30 14:15:52 +00:00
std::vector<unsigned long long> primes;
primes.push_back(2);
unsigned long long num = 3;
2016-11-30 14:15:52 +00:00
// Use for to accomodate for OpenMP
#pragma omp parallel
{
do {
unsigned long long myNum;
#pragma omp critical
{
myNum = num;
num += 2;
}
bool isPrime = true;
2016-11-30 20:37:10 +00:00
#pragma omp barrier
for(auto i : primes) {
if(i > myNum / 2) break;
if(myNum % i == 0) {
isPrime = false;
break;
}
}
if(isPrime) {
#pragma omp critical
{
primes.push_back(num);
std::cout << num << std::endl;
}
}
} while(run);
2016-11-30 14:15:52 +00:00
}
2016-11-30 13:32:06 +00:00
return 0;
}