/*  
*  Для прогрева всех ядрывшек, теперь c OpenMP
*   gcc-4.4.2 primesum.c -std=gnu99 -O3 -ffast-math -lgmp -fopenmp -o pimesum
*/

#include <gmp.h>

#define POW2_128 "340282366920938463463374607431768211456"

int main (int argc, char **argv)
{
  int k;
  mpz_t i, sum;
  mpz_t LIMIT;

  mpz_init (i);
  mpz_init (sum);
  mpz_init_set_str(LIMIT, POW2_128, 10);

#pragma omp parallel
{
#pragma omp do private (i, k)
                do {
                    mpz_add_ui(i, i, 1ULL);
                    k = mpz_cmp(LIMIT, i);
                    if (mpz_probab_prime_p (i, 10) == 2 ) {
#pragma omp shared (sum)
#pragma omp critical (sum)
                      mpz_add(sum, i, sum);
                      gmp_printf("Prime: %Zd, Sum  %Zd\n", i, sum);
#pragma omp flush(sum)
                   }
               } while ( k != 0 );
#pragma omp end parallel
}
  gmp_printf("%Zd\n", sum);
  mpz_clear (i);
  mpz_clear (sum);

return 0;
}