Como: Converter um OpenMP paralelo para que o loop use o tempo de execução de concorrência
Este exemplo mostra como converter um loop básico que usa o OpenMP paralela e políticas de para para usar o algoritmo de concurrency::parallel_for de tempo de execução de simultaneidade.
Exemplo
Este exemplo usa OpenMP e o tempo de execução de simultaneidade para calcular a contagem de números primos em uma matriz de valores aleatórios.
// concrt-omp-count-primes.cpp
// compile with: /EHsc /openmp
#include <ppl.h>
#include <random>
#include <array>
#include <iostream>
using namespace concurrency;
using namespace std;
// Determines whether the input value is prime.
bool is_prime(int n)
{
if (n < 2)
return false;
for (int i = 2; i < n; ++i)
{
if ((n % i) == 0)
return false;
}
return true;
}
// Uses OpenMP to compute the count of prime numbers in an array.
void omp_count_primes(int* a, size_t size)
{
if (size == 0)
return;
size_t count = 0;
#pragma omp parallel for
for (int i = 0; i < static_cast<int>(size); ++i)
{
if (is_prime(a[i])) {
#pragma omp atomic
++count;
}
}
wcout << L"found " << count
<< L" prime numbers." << endl;
}
// Uses the Concurrency Runtime to compute the count of prime numbers in an array.
void concrt_count_primes(int* a, size_t size)
{
if (size == 0)
return;
combinable<size_t> counts;
parallel_for<size_t>(0, size, [&](size_t i)
{
if (is_prime(a[i])) {
counts.local()++;
}
});
wcout << L"found " << counts.combine(plus<size_t>())
<< L" prime numbers." << endl;
}
int wmain()
{
// The length of the array.
const size_t size = 1000000;
// Create an array and initialize it with random values.
int* a = new int[size];
mt19937 gen(42);
for (size_t i = 0; i < size; ++i) {
a[i] = gen();
}
// Count prime numbers by using OpenMP and the Concurrency Runtime.
wcout << L"Using OpenMP..." << endl;
omp_count_primes(a, size);
wcout << L"Using the Concurrency Runtime..." << endl;
concrt_count_primes(a, size);
delete[] a;
}
O exemplo produz a seguinte saída.
O algoritmo de parallel_for e o OpenMP 3,0 permitem o tipo do índice ser um tipo integral sinal ou um tipo integral sem sinal.O algoritmo de parallel_for também garante que o intervalo especificado não transborda um tipo com sinal.As versões 2,0 e 2,5 de OpenMP permitem tipos integral com sinal o índice somente.OpenMP também não valida o intervalo de índice.
A versão deste exemplo que usa o tempo de execução de simultaneidade também usa um objeto de concurrency::combinable no lugar de política de atômico para incrementar o valor do contador sem exigir a sincronização.
Para obter mais informações sobre parallel_for e outros algoritmos paralelos, consulte Algoritmos paralelos.Para obter mais informações sobre a classe combinable, consulte Contêiner e objetos paralelos.
Este exemplo altera anterior para atuar no em um objeto de std::array em vez de uma matriz nativo.Porque as versões 2,0 e 2,5 de OpenMP permitem o índice inteiro com sinal necessários somente em uma compilação de parallelfor , você não pode usar iteradores para acessar paralelamente os elementos de um contêiner de (STL) da biblioteca do modelo de padrão.A biblioteca (PPL) dos padrões de paralela fornece o algoritmo de concurrency::parallel_for_each , que executa tarefas, paralelamente, em um contêiner iterativo como aqueles fornecidos pelo STL.Usa a mesma lógica do particionamento que o algoritmo de parallel_for usa.O algoritmo de parallel_for_each lembra o algoritmo de STL std::for_each , exceto que o algoritmo de parallel_for_each executa as tarefas simultaneamente.
// Uses OpenMP to compute the count of prime numbers in an
// array object.
template<size_t Size>
void omp_count_primes(const array<int, Size>& a)
{
if (a.size() == 0)
return;
size_t count = 0;
int size = static_cast<int>(a.size());
#pragma omp parallel for
for (int i = 0; i < size; ++i)
{
if (is_prime(a[i])) {
#pragma omp atomic
++count;
}
}
wcout << L"found " << count
<< L" prime numbers." << endl;
}
// Uses the Concurrency Runtime to compute the count of prime numbers in an
// array object.
template<size_t Size>
void concrt_count_primes(const array<int, Size>& a)
{
if (a.size() == 0)
return;
combinable<size_t> counts;
parallel_for_each(begin(a), end(a), [&counts](int n)
{
if (is_prime(n)) {
counts.local()++;
}
});
wcout << L"found " << counts.combine(plus<size_t>())
<< L" prime numbers." << endl;
}
Compilando o código
Copie o código de exemplo e cole-o em um projeto do Visual Studio, ou cole em um arquivo denominado concrt-omp-count-primes.cpp e execute o seguinte comando em uma janela de prompt de comando do Visual Studio.
cl.exe /EHsc /openmp concrt-omp-count-primes.cpp