//Sorting Using Threads (Assignment).

#include<stdio.h>
#include"struct1.c"
#include"mpi.h"

//main Function

/*void  print(double  a[][m1],int  r,int c)        //function for array printing
{
  int  i,j;
  for(i=0;i<r;++i)
    {
      for(j=0;j<c;++j)
	{
	  fprintf(fp,"%15.12lf\t",a[i][j]);
	 
	}
      fprintf(fp,"%s","\n\n");
    }
}*/

void step0(char f1[],char f2[],char f3[],int row_count,int col_count)
{
  int r,c;
  FILE *fp1,*fp2,*fp3;
  fp1 = fopen(f1,"r");
  fp2 = fopen(f2,"r");
  fp3 = fopen(f3,"r");

  if(fp1 == NULL || fp2 == NULL || fp3 == NULL)
    {
      printf("\n Error in opening file:");
      exit(1);
    }

  for(c=0;c<col_count;c+=3)
    {
      for(r= 0;r<row_count;r++) 
	{
	  fscanf(fp1, "%lf",&a[r][c]);
  	  fscanf(fp2, "%lf",&a[r][c+1]);
	  fscanf(fp3, "%lf",&a[r][c+2]);
	}
    }

  fclose(fp1);
  fclose(fp2);
  fclose(fp3);
}
  
int main(int argc , char *argv[])                                        
{
  int row_count,col_count,r,c,i=0,p;
  char f1[25];
  char f2[25];
  char f3[25];

  int rank,size;
  MPI_Status status;   

  MPI_Init(&argc,&argv);
  MPI_Comm_rank(MPI_COMM_WORLD,&rank);
  MPI_Comm_size(MPI_COMM_WORLD,&size);


  /* if(argc != 6)
    {
      printf("\nError in Input\n\n");
      return 0;
      }*/

  //fp=fopen("output.txt","w");
  //for(i=0;i<argc;i++)
  // printf("%s\n",argv[i]);


  //printf("%s",argv[4]);

  sscanf(argv[1],"%d",&col_count);   // no. of columns
  sscanf(argv[2],"%d",&row_count);   // no. of rows
  strcpy(f1,argv[3]);                // input file 1
  strcpy(f2,argv[4]);                // input file 2
  strcpy(f3,argv[5]);                // input file 3
  
  fprintf(fp,"%s","\nSTEP:=0\n");                //For scanning Element.
  step0(f1,f2,f3,row_count,col_count);
  //print(a,row_count,col_count);  

  void step2()
  {
    int i,j,revcnt,procCNT=col_count,k=0,p,c;
    int send_cnt = (row_count * col_count);
    int recv_cnt = (row_count * col_count);
    
    if(rank == 0)   
      {
	MPI_Bcast(&a,send_cnt,MPI_DOUBLE,0,MPI_COMM_WORLD);
      }

    for(p=0;p<procCNT;p++)
      {
	if(rank = p)
	  {
	    MPI_Recv(&temp1,recv_cnt,MPI_DOUBLE,0,1,MPI_COMM_WORLD,&status);
	    printf("Process:%d\n",rank);

 	    for(i=0;i<row_count;i=i+col_count)
	      {
		for(j=0;j<col_count;j++)
		  {
		    temp1[k][j] = a[i+j][p];
		    //MPI_Gather(&a,send_cnt,MPI_DOUBLE,&temp1[k][j],recv_cnt,MPI_DOUBLE,MPI_COMM_WORLD);
		  }
		k = k+1;
	      }
	  }
	MPI_Gather(&a,send_cnt,MPI_DOUBLE,&temp1,recv_cnt,MPI_DOUBLE,0/*(root process)*/,MPI_COMM_WORLD);
      }
    MPI_Finalize();
    return;
  }
  //  print(temp1,row_count,col_count);
  //fclose(fp);
  MPI_Finalize();
  return 0;
}