O código abaixo é fruto de um trabalho que estou desenvolvendo, basicamente é a multiplicação de uma matriz quadrada, porém, os resultados que eu tive paralelizando a aplicação com a API OpenMP foram superiores aos resultados que obtive usando SIMD da mesma API.
O que estou fazendo de errado? é a sintaxe?
Algumas informações que podem ser pertinente em identificar o problema: Estou usando o compilador da intel atravéz da IDE do visual studio, o OpenMP do visual studio é versão 2.0 (que não suporta SIMD) mas acho que é o 4.0 que vem com o compilador que está sendo usado. Enfim, pra mim é uma atividade nova (processamento paralelo) então se puderem esclarecer as coisas agradeceria de coração. Segue o código:
#include "stdafx.h"
#include <iostream>
#include <time.h>
#include <omp.h>
using namespace std;
int lin = 800, col = 800; // Valores de linha e coluna
int main()
{
// --------------------------------------
// Cria a matriz 1
int** m1 = new int*[lin];
for (int i = 0; i < lin; ++i)
m1[i] = new int[col];
// --------------------------------------
// --------------------------------------
// Cria a matriz 2
int** m2 = new int*[lin];
for (int i = 0; i < lin; ++i)
m2[i] = new int[col];
// --------------------------------------
// --------------------------------------
// Cria a matriz resposta
int** res = new int*[lin];
for (int i = 0; i < lin; ++i)
res[i] = new int[col];
// --------------------------------------
cout << "criou matrizes" << endl;
//PREENCHE m1 e m2
// ----------------------------------------------------------------------------
// BLOCO PARALELO
#pragma omp simd collapse (2)
for (int i = 0; i < lin; ++i) {
for (int j = 0; j < lin; ++j) {
m1[i][j] = (i + 1);
}
}
// FIM DO BLOCO PARALELO
// BLOCO PARALELO
#pragma omp simd collapse (2)
for (int i = 0; i < lin; ++i) {
for (int j = 0; j < lin; ++j) {
m2[i][j] = (i + 1);
}
}
// FIM DO BLOCO PARALELO
cout << "preencheu" << endl;
// ----------------------------------------------------------------------------
//faz a magica rolar
clock_t timer = clock(); //valores de marcação de tempo
// ----------------------------------------------------------------------------
cout << "iniciou" << endl;
#pragma omp simd collapse (2)
for (int i = 0; i < lin; i++)
{
for (int j = 0; j < lin; j++)
{
res[i][j] = 0;
for (int k = 0; k < lin; k++)
res[i][j] += m1[i][k] * m2[k][j];
}
}
cout << "finalizou" << endl;
// ----------------------------------------------------------------------------
//marca tempo final e exibe
timer = clock() - timer;
cout << "Programa Finalizado em " << ((float)timer) / CLOCKS_PER_SEC << " Segundos" << endl;
system("Pause");
}
// This code is contributed
// by Soumik Mondal