Użytkownik "slawek" <s...@h...pl> napisał w wiadomości grup
dyskusyjnych:4f50f420$0$1231$6...@n...neostrada
.pl...
> Wyniki, w sekundach, są takie jak niżej. Program taki jak jeszcze niżej.
>
> gfortran -Ofast test-omp.f95 CPU_TIME = 0.20312500
> gfortran test-omp.f95 CPU_TIME = 1.9062500
> gfortran -fopenmp test-omp.f95 CPU_TIME = 7.0000000
> gfortran -Ofast -fopenmp .... CPU_TIME = 1.0625000 - jeżeli
> wątki tworzy się tylko raz (tj. $omp parallel obejmuje CAŁY program)

Po przeniesieniu do C i kompilacji MSVC 2010 jest nieco lepiej - ale i tak w
rebusie

bez OpenMP jest np. 8 sekund
z OpenMP jest 12 sekund

(Program jest trochę inny, bo nie ma complex i pętle z większą ilością
iteracji.)




/***************************************************
****************************************************
************************
!
! Program test-omp
!
! /openmp
!
****************************************************
****************************************************
***********************/

#include <stdio.h>
#include <stdlib.h>
#include <time.h>

// stale matematyczne i fizyczne (CODATA 2006)

const double pi = 3.14159265358979323846264338327950288419716939938;
const double epsilon0 =
8.85418781762038985053656303171075026060837016660E-1
2; // przenikalnosc
elektryczna prozni [F/m]
const double c = 299792458.0;
// predkosc swiatla [m/s]

static int n, m;

void setup(double vec[])
{
int i;
for(i = 0; i < n; i++)
vec[i] = pi*c*epsilon0;
}

void solve(double v1[], double v2[])
{
int i;

#pragma omp parallel
{

#pragma omp for schedule(static,100)
for(i = 0; i < n; i++)
v2[i] = v1[i]/(c*c) + epsilon0 + pi;
}
}

int main(void)
{
clock_t stamp0; // pomiar czasu - stempel 0
clock_t stamp1; // pomiar czasu - stempel 1

// takie duze tablice sa
# define nmax 10000

// problem - nie ma typu complex w MSVC - zamiast
complex
// będzie używany double :(
double vec[2][nmax];
int i1,i2,i,j;
double sum;

stamp0 = clock();

n = 1000;
m = 1000000;

if(n < nmax)
{
i1 = 0;
i2 = 1;

setup(vec[i1]);
for(j = 0; j < m; j++)
{
solve(vec[i1],vec[i2]);
i1 = (i1+1) % 2;
i2 = (i2+1) % 2;
}
}
else
puts("error: dimension(s)");

stamp1 = clock();

printf("CPU time = %lf\n", (double)(stamp1-stamp0)/CLOCKS_PER_SEC);

sum = 0.;
for(i = 0; i < n; i++)
sum += 1./(1.+vec[i1][i]*vec[i1][i]);

printf("just for fun %lf\n",sum);

getchar();
return 0;
}