nortega/indivisible-legacy
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Cancel if `primes.list[i] >= numRoot'.

Browse Source 
It has been solved!!!
This commit is contained in:
Nicolás A. Ortega 2017-02-01 23:12:23 +01:00
parent b9cafadf8e
commit a5b68282b3
No known key found for this signature in database
GPG Key ID: 614272579C2070D1
1 changed files with 8 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
src/main.c
View File

@ -190,23 +190,18 @@ int main(int argc, char *argv[]) {
* number we're analyzing. Also, skip `2' since we're not testing even
* numbers.
*/
/*for(size_t i = 1; mpz_cmp(primes.list[i], numRoot) <= 0; ++i) {
// If `num' is divisible by a prime then go to the next number
if(mpz_divisible_p(num, primes.list[i]) != 0) {
isPrime = false;
break;
}
}*/
#pragma omp parallel
#pragma omp parallel
{
#pragma omp for
#pragma omp for
for(size_t i = 1; i < primes.end; ++i) {
if(mpz_divisible_p(num, primes.list[i])) {
#pragma omp atomic write
if(mpz_cmp(primes.list[i], numRoot) >= 0) {
#pragma omp cancel for
} else if(mpz_divisible_p(num, primes.list[i])) {
#pragma omp atomic write
isPrime = false;
#pragma omp cancel for
#pragma omp cancel for
}
#pragma omp cancellation point for
#pragma omp cancellation point for
}
}