Параллельное программирование в стандарте MPI. Баканов В.М - 52 стр.

UptoLike

ВУЗ: 

МИРЭА | Москва

Составители: 

Рубрика: 

- 52 -
rc = MPI_Send(&homepi, 1, MPI_DOUBLE, MASTER, i, MPI_COMM_WORLD);
if (rc != MPI_SUCCESS) // if error MPI_Send function
printf(“%d: Send failure on round %d\n”, taskid, i);
}
else // I am MASTER
{
pisum = 0;
for (n=1; n<numtasks; n++)
{
rc = MPI_Recv(&pirecv, 1, MPI_DOUBLE, MPI_ANY_SOURCE, i,
MPI_COMM_WORLD, &status);
if (rc != MPI_SUCCESS) // if error MPI_Recv function
printf(“%d: Receive failure on round %d\n”, taskid, i);
pisum = pisum + pirecv;
}
/* Master calculates the average value of pi for this iteration */
pi = (pisum + homepi)/numtasks;
/* Master calculates the average value of pi over all iterations */
avepi = ((avepi * i) + pi) / (i + 1);
t2=MPI_Wtime(); // fix end time calculated
printf(“%9d throws aver.value of PI= %.12lf (rel.error= %.5lf %%, time= %.3lf sec)\n”,
(DARTS * (i+1)), avepi, 1.0e2*(pi_prec-avepi)/pi_prec, t2-t1);
}
} // end for (i=0; i<ROUNDS; i++)
MPI_Finalize();
return 0;
} // end of MAIN function
// end of PI_02.C program
Выдача программы (вычислительный кластер кафедры ИТ-4, число ‘вы-
стрелов’ от 1000 до 10’000) приведена ниже (напомним, что каждое вычис-
ленное значение следует рассматривать как одну из реализаций случайной
величины, нормально распределенной вокруг истинного значения
π
):
LAM 7.0.6/MPI 2 C++ - Indiana University
10000 throws aver.value of PI= 3.119600000000 (rel.error= 0.70005 %, time= 0.003 sec)
20000 throws aver.value of PI= 3.137800000000 (rel.error= 0.12072 %, time= 0.006 sec)
30000 throws aver.value of PI= 3.140000000000 (rel.error= 0.05070 %, time= 0.037 sec)
40000 throws aver.value of PI= 3.144100000000 (rel.error= -0.07981 %, time= 0.040 sec)
50000 throws aver.value of PI= 3.143440000000 (rel.error= -0.05880 %, time= 0.043 sec)
60000 throws aver.value of PI= 3.139333333333 (rel.error= 0.07192 %, time= 0.046 sec)
70000 throws aver.value of PI= 3.141200000000 (rel.error= 0.01250 %, time= 0.048 sec)
80000 throws aver.value of PI= 3.145250000000 (rel.error= -0.11642 %, time= 0.051 sec)
90000 throws aver.value of PI= 3.144933333333 (rel.error= -0.10634 %, time= 0.053 sec)
…
100000 throws aver.value of PI= 3.143760000000 (rel.error= -0.06899 %, time= 0.056 sec)
…
200000 throws aver.value of PI= 3.143840000000 (rel.error= -0.07154 %, time= 0.082 sec)
…
400000 throws aver.value of PI= 3.142310000000 (rel.error= -0.02283 %, time= 0.133 sec) 
                                           - 52 -

 rc = MPI_Send(&homepi, 1, MPI_DOUBLE, MASTER, i, MPI_COMM_WORLD);
 if (rc != MPI_SUCCESS) // if error MPI_Send function
  printf(“%d: Send failure on round %d\n”, taskid, i);
}
else // I am MASTER
{
 pisum = 0;
 for (n=1; n

Страницы