Como converter um loop OpenMP paralelo para loop para usar o tempo de execução de simultaneidade
Este exemplo demonstra 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 à esquerda 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 permite o tipo de índice ser um tipo integral um tipo integral assinado ou não assinado. O algoritmo de parallel_for também garante que o intervalo especificado não excederá um tipo assinado. As versões 2,0 e 2,5 OpenMP permitem tipos integrais assinados de índice. 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 da política de atomic para incrementar o valor do contador sem exigir a sincronização.
Para obter mais informações sobre como parallel_for e outros algoritmos paralelos, consulte Algoritmos paralelos. Para obter mais informações sobre a classe combinable, consulte Contêineres e objetos em paralelo.
Este exemplo altera anterior para atuar em um objeto de std::array em vez de uma matriz nativo. Como as versões 2,0 e 2,5 OpenMP permitem o índice integral assinado tipo somente em uma compilação de parallelfor , você não pode usar iteradores em paralelo para acessar os elementos de um contêiner de (STL) da biblioteca do modelo padrão de. A biblioteca (PPL) de padrões de paralelo fornece o algoritmo de concurrency::parallel_for_each , que executa tarefas, em paralelo, em um contêiner iterativo como os que são fornecidos por STL. Usa a mesma lógica de particionamento que o algoritmo de parallel_for usa. O algoritmo de parallel_for_each se assemelha ao algoritmo 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-o em um arquivo chamado concrt-omp-count-primes.cpp e execute o comando a seguir em uma janela de prompt de comando do Visual Studio.
cl.exe /EHsc /openmp concrt-omp-count-primes.cpp