Параллельное программирование

[BleStaR]

касательно задания матрица не квадратная и количество процессов эквивалентно количеству строк.
что прога дает при матрице 5 на 5 и при 5-и процессах??
(сам не могу проверить, т.к. сижу на компе без сишки)

[mages]

матрица квадратная, значение N=M в данном случае.
а кол-во процессов должно быть кратно N. т.е. если процессов 6 то матрица должна быть соответственно 2 на 2 ,3 на 3 и 6 на 6.
если процессов 4 то 4на4 и 2на2.
если 5 то 5на5 только.
а вот если взять количество процессов 6 и матрицу 3 на 3 , то программа не может найти максимальные элементы в строках матрицы. аналогично если процессов 8 то не может найти в матрице 4на4 и 2на2.

[mages]

находит только тогда когда N равно кол-ву процессов.

[BleStaR]

Первое: прога скорее всего и не должна работать при N меньшем количества процессов. Т.к. MPI_Scatter раздает части твоего массива процессам, а ведь строк меньше чем процессов! Вопрос что получат остальные? На это я ответить пока не могу(возьможно NULL -что вызывает ошибки, которые обрабатываются самой функцией, поэтому тебе про них не говорят). Но ошибка явно в этом!
Второе:
По заданию N кратно количеству процессов - если не изменяет память, и здравый смысл, то вроде бы N должно делится на цело на количество процессов. Т.е. N больше количества процессов. Например Строк 8 процессов 4,2,1.

[BleStaR]

Ну и касательно того случая когда количество процессов будет меньше чем количество строк.
Тебе надо будет организовать систему которая следила бы за тем какой из процессов закончил свою работу, и снова его загрузить.

[mages]

Иправил прогу.
при равном находит нормально N равно процессам.
а при 3х процессах и матрице 6на6 находит неправильно.
кажетс ошибка в МПИ_Скаттере либо в чунке

[BleStaR]

Выложи исправленый вариант. Возможно ты прав по поводу скатерра. Пробовал отследить прием-оправку?

[mages]

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include <mpi.h>


void Init(FILE *f, int m,int n, double *B)
{int i,j; 
 for(i=0;i<m;i++) {
 for(j=0;j<n;j++)fscanf(f,"%lf",&B[i*m+j]);}
 fscanf(f,"\n"); 
 for(i=0;i<m;i++)
 {
  for(j=0;j<n;j++)
  printf("%4.lf",B[i*m+j]);
  printf("\n");
 }
}
  

  double Max( double *B,int m,int n)
{
 int i,j;
 double  max1;
 for(i = 0; i < n*chunk;) { 
 max1=B[0];
 for(j=i;j<=n;j++) if(B[j]>max1) max1=B[j];
 i=i+n;
 printf("max %i Stroki=%lf\n",rk*chunk+i,max1);
 }
 }  
  
  void print(double *B,int m,int n)
  {
   int i,j;
   for(i=0;i<m;i++){  
   for(j=0;j<n;j++) printf("%4.0lf",B[i*m+j]); printf("\n");
  }
 }
  
 
  int main (int argc,char *argv[])
  { 
   FILE *f=fopen("input.txt","r");
   int i,np,rk,m,n,chunk=1; 
   double *B,*Bloc,t;
  
  MPI_Init(&argc,&argv);
  MPI_Comm_size(MPI_COMM_WORLD,&np);
  MPI_Comm_rank(MPI_COMM_WORLD,&rk);
  if(rk==0)
  { 
   fscanf(f,"%d ",&m); fscanf(f,"%d\n",&n);
   chunk=m/np;
   B=(double *)malloc(m*n*sizeof(double));
  }
  MPI_Bcast(&m,1, MPI_INT,0, MPI_COMM_WORLD);
  MPI_Bcast(&chunk,1, MPI_INT,0, MPI_COMM_WORLD);
  Bloc=(double *)malloc(chunk*n*sizeof(double));
  
  if (rk==0)
  {
   Init(f,m,n,B);
   t=MPI_Wtime();
  }
  MPI_Scatter(B,m*chunk,MPI_DOUBLE,Bloc,m*chunk,MPI_DOUBLE,0,MPI_COMM_WORLD);
  Max(Bloc,chunk,n);
  MPI_Gather(Bloc,m*chunk,MPI_DOUBLE,B,m*chunk,MPI_DOUBLE,0,MPI_COMM_WORLD);
  
  if(rk==0)
  {
   t=MPI_Wtime()-t;
   printf("%lf sec. \n",t);
   }
   return 0;
}

[mages]

теперь программа требует маленькой доработки.
поиск,когда N кратно процессам и не равно. при 2х процессах матрице 6на6 или 4на4.