Mercurial > hg > what-is-octave

comparison PoliMI2012/examples/Pi.m @ 4:0a1567794b40

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[mq]: folder_reorganization
author Carlo de Falco <cdf@users.sourceforge.net>
date Thu, 08 Nov 2012 09:55:34 +0100
parents
children
comparison
equal deleted inserted replaced
3:acc25d26e9e5 4:0a1567794b40
1 T=clock;
2 MPI_ANY_SOURCE = -1;
3 MPI_Init ();
4 MPI_COMM_WORLD = MPI_Comm_Load ("NEWORLD");
5 rnk = MPI_Comm_rank (MPI_COMM_WORLD);
6 siz = MPI_Comm_size (MPI_COMM_WORLD);
7 SLV = logical(rnk);
8 MST = ~ SLV;
9 width=1/N; lsum=0;
10 i=rnk:siz:N-1;
11 x=(i+0.5)*width;
12 lsum=sum(4./(1+x.^2));
13 TAG=7;
14 if SLV
15 MPI_Send (lsum, 0, TAG, MPI_COMM_WORLD);
16 else
17 Sum =lsum;
18 for slv=1:siz-1
19 lsum = MPI_Recv (MPI_ANY_SOURCE, TAG, MPI_COMM_WORLD);
20 Sum += lsum;
21 endfor
22 endif
23 MPI_Finalize ();
24