OpenMP set_num_threads () não está funcionando

Eu estou escrevendo um programa paralelo usando o OpenMP em C ++.

Eu quero controlar o número de threads no programa usandoomp_set_num_threads(), mas não funciona.

#include <iostream>
#include <omp.h>
#include "mpi.h"

using namespace std;

int myrank;
int groupsize;
double sum;
double t1,t2;
int n = 10000000;

int main(int argc, char *argv[])
{
    MPI_Init( &argc, &argv);
    MPI_Comm_rank( MPI_COMM_WORLD, &myrank );
    MPI_Comm_size(MPI_COMM_WORLD,&groupsize);

    omp_set_num_threads(4);

    sum = 0;
    #pragma omp for  reduction(+:sum)
    for (int i = 0; i < n; i++)
        sum+= i/(n/10);

    cout<<"sum="<<sum<<endl;
    cout<<"threads="<<omp_get_num_threads()<<endl;

    MPI_Finalize();
    return 0;
}

O programa produz:

sum = 4.5e+007
threads=1

Como controlar o número de threads?

questionAnswers(5)

yourAnswerToTheQuestion