1 | rank = MPI_COMM_RANK( comm ) |
2 | if (rank .eq. 0) then |
3 | read *, n |
4 | end if |
5 | call MPI_BCAST(n, 1, MPI_INTEGER, 0, comm ) |
6 | lo = rank*n+1 |
7 | hi = lo+n-1 |
8 | sum = 0.0d0 |
9 | do i = lo, hi |
10 | sum = sum + 1.0d0 / i |
11 | end do |
12 | call MPI_REDUCEALL( sum, sumout, 1, MPI_DOUBLE, |
13 | & MPI_ADD_DOUBLE, comm) |