修論ツールの履歴(No.3)

履歴一覧
差分を表示
現在との差分を表示
ソースを表示
修論ツールへ行く。
- 1 (2018-10-29 (月) 11:28:19)
- 2 (2018-11-26 (月) 06:05:28)
- 3 (2018-11-26 (月) 06:38:11)
- 4 (2018-12-04 (火) 08:13:18)

ーーーーー以下、ソースコードーーーーー

#include <stdio.h>
#include <math.h>
#include <mpi.h>
void main(int argc, char* argv[]){
  double PI25DT = 3.141592653589793238462643;
  double mypi, pi, h, sum, x, f, a;
  double t1, t2, t0, t_w;
  int n, myid, numprocs, i, rc;
  int ierr;
  ierr = MPI_Init(&argc, &argv);
  ierr = MPI_Comm_rank(MPI_COMM_WORLD, &myid);
  ierr = MPI_Comm_size(MPI_COMM_WORLD, &numprocs);
  if(myid==0){
     printf(“Enter the number of intervals \n”);
     scanf(“%d”, &n);
     printf(“n=%d \n”, n)
  }
  ierr = MPI_Barrier(MPI_COMM_WORLD);
  t1 = MPI_Wtime();
  ierr = MPI_Bcast(&n, 1, MPI_INT, 0, MPI_COMM_WORLD);
  h = 1.0 / n;
  sum = 0.0;
  for(i=myid+1; i<=n; i+=numprocs){
     x = h * (i – 0.5);
     sum = sum + 4.0 / (1.0 + x*x);
  }
  mypi = h * sum;
  MPI_Reduce(&mypi, &pi, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);
  if(myid == 0){
     printf(“pi is approximately: %18.16lf Error is: %18.16lf \n”, pi, 
fabs(pi-PI25DT));
     ierr = MPI_Barrier(MPI_COMM_WORLD);
     t2 = MPI_Wtime();
     t0 = t2 – t1
     ierr = MPI_Reduce(&t0, &t_w, 1, MPI_DOUBLE, MPI_MAX, 0, MPI_COMM_WORLD);
     printf(“execution time = : %8.4lf [sec.] \n”, t_w);
  }
  rc = MPI_Finalize();
}

コンパイル方法

mpicc pisample.c -o pisample

実行方法

mpirun --hostfile host -np 3 ./pisample

↑のnpの後の3は並列数を表す

途中でEnter the number of intervalsと言われて、整数を入力する。しかし、どういう意味があるかまだ理解していない。

参考サイト

http://aero.me.tut.ac.jp/Lectures/CFD/%89%89%8FK%82l%82o%82h.pdf
http://www.r-ccs.riken.jp/r-ccssite/wp-content/uploads/2013/07/ss13_enshu2.pdf
https://www.cc.u-tokyo.ac.jp/public/VOL9/special/10.pdf

修論ツール の履歴(No.3)

参考サイト