#include "mpi.h"
#include <stdio.h>
#include <math.h>

#define N 8192
// 行列のサイズ

int main(int argc, char **argv){
  int n, myid, numprocs, i, j;
  double A[N][N], x[N], y[N];
  int f;

  double Stime, Etime;

  MPI_Init(&argc, &argv);
  MPI_Comm_size(MPI_COMM_WORLD, &numprocs);
  MPI_Comm_rank(MPI_COMM_WORLD, &myid);

  for(i=0; i<N; i++){
    for(j=0; j<N; j++){
      A[i][j] = 1;
    }
    x[i] = 1;
  }
    // とりあえず行列とベクトルの中身は全部1にしてみる

  Stime = MPI_Wtime();
  
  for(i=0; i<N; i++){
    y[i] = 0;
    for(j=0; j<N; j++){
      y[i] += A[i][j]*x[j];
    }
  }

  Etime = MPI_Wtime();

  if(myid==0){
      // myidが0, つまり0番目のプロセッサのみこの命令を実行
    f=1;
    for(i=0; i<N; i++){
      if(fabs(y[i]-N)>1e-7) f=0;
    }
      // 答えの検証。合ってたらfは1

    printf("%lf %lf\n",(Etime-Stime), MPI_Wtick() );
    printf("f=%d\n", f);
  }
  MPI_Finalize();
  
  return 0;
}