[MPI/C] Problema con Recv e Send

[PieceSpencer]

Buonasera a tutti ragazzi, volevo chiedervi una cosa su questo programma che dovrei fare per l'universtita', ma purtroppo va in deadlock alla prima Recv.

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


int main (int argc, char *argv[]){
	float *array_num,*array_sott;
	int i,j,start,n_fuori,z;
	char *ptr;
	int n=argc; //numero di reali immessi in input
	float resto,sum,sum_parz=0;
	int n_somme,numero_piu,num,menum,nproc,tag;
	MPI_Status info;
	
	MPI_Init(&argc,&argv);
	MPI_Comm_rank(MPI_COMM_WORLD,&menum);
	MPI_Comm_size(MPI_COMM_WORLD,&nproc);
	
	
   n_somme=n/nproc; //mi calcola il numero di somme da svolgere 
   
   
   resto=n%nproc;
   MPI_Bcast(&n_somme,1,MPI_INT,0,MPI_COMM_WORLD);
   
   if(resto > 0){ //significa che le somme sono di numero dispari (considerando ovviamente i processori di numero pari)
     numero_piu= (n-1)-(nproc * n_somme); //numero di somme in piu'
     MPI_Bcast(&numero_piu,1,MPI_INT,0,MPI_COMM_WORLD);
   } 
   if(menum==0) {
   	tag=10;
   	sum=(atof(argv[1]));
   	for(i=2;i<=n_somme;i++){
   	sum=sum+(atof(argv[i]));
   	
   	
     	}
     	
    array_num=malloc(n_somme+sizeof(float));
 
    for(j=0;j<n_somme;j++){
      array_num[j]=atof(argv[i]);
      i++;
     	
    } 	
    
     for(z=0;z<n_somme;z++){
      printf("\na[%d] vale %f\n",z,array_num[z]);
     	fflush(stdout);
    }
    
     	/*ptr=argv[i];
     	printf("%s",ptr);
     	printf("\n%p",ptr);
     	printf("\n%p\n",argv[i]);
     	printf("%c\n",(const char)(*ptr+sizeof(char)));*/
     	
     	
     	sum_parz=0;
     	//printf("\n%d\n",n_somme);
     	printf("\n il tag e'  %d\n",tag);
     	
     	MPI_Send(&array_num,n_somme,MPI_FLOAT,0,tag,MPI_COMM_WORLD);
     	
      //MPI_Send(&sum_parz,1,MPI_FLOAT,0,tag,MPI_COMM_WORLD);
      
      /*tag=tag+1;
      MPI_Recv(&sum_parz,1,MPI_FLOAT,1,tag,MPI_COMM_WORLD,&info);
      */
      
     	//sum=sum+sum_parz;
      //printf("\n 1 la somma e' %f\n",sum);
  
   }	
   else {
   	tag=10;
      printf("\n il tag e'  %d\n",tag);
      //printf("\n la somma e' %f\n",sum_parz);
      //printf("\n%d\n",n_somme);
      
     	MPI_Recv(&array_num,n_somme,MPI_FLOAT,0,tag,MPI_COMM_WORLD,&info);
     	
   	
   	printf("\n%d\n",n_somme);
      //MPI_Recv(&sum_parz,1,MPI_FLOAT,0,tag,MPI_COMM_WORLD,&info);
      //printf("\n la somma e' %f\n",sum_parz);
      
      sum_parz=0;
      //printf("\n la somma e' %f\n",sum_parz);
      sum_parz=array_num[0];
      //printf("\n la somma e' %f\n",sum_parz);
      for(j=1;j<=n_somme;j++) {
      	sum_parz=sum_parz+array_num[j];
      	printf("\n%f\n",sum_parz);
      }
      printf("\n 2 la somma e' %f\n",sum_parz);
      /*tag=tag+1;
      MPI_Send(&sum_parz,1,MPI_FLOAT,1,tag,MPI_COMM_WORLD);
     */
   }
	
	MPI_Finalize();
	
	return 0;
	
 }

In pratica alla prima Recv nell'else, mi va in deadlock e non capisco il perche'. Sono arrivato alla conclusione che sia qualcosa legato ad array_num e il suop numero di elementi (n_somme)
Ah dimenticavo....devo calcolare la somma dei primi n numeri reali .

[Kaya]

Originariamente inviato da PieceSpencer

Buonasera a tutti ragazzi, volevo chiedervi una cosa su questo programma che dovrei fare per l'universtita', ma purtroppo va in deadlock alla prima Recv.
[CODE
int n_somme,numero_piu,num,menum,nproc,tag;
MPI_Status info;

MPI_Init(&argc,&argv);
MPI_Comm_rank(MPI_COMM_WORLD,&menum);
MPI_Comm_size(MPI_COMM_WORLD,&nproc);
n_somme=n/nproc; //mi calcola il numero di somme da svolgere
[/code]

.

Purtroppo non conosco MPI quindi non saprei aiutarti, ma non mi torna questo punto
n_somme= n/nproc;
definisci n_somme come int ma la divisione potrebbe ritornarti anche un float.. è solo perchè ti serve la parte intera?

[melko]

non sono esperto di MPI ma:

Originariamente inviato da PieceSpencer

MPI_Send(&array_num,n_somme,MPI_FLOAT,0,tag,MPI_COMM_WORLD);

qui fai una send da rank 0 verso rank 0, mentre immagino quello che vorresti fare sia una send da rank 0 a tutti gli altri (che nel gergo di MPI sarebbe una Scatter).
Oltretutto occhio che ci potrebbero essere altre sviste, ad una rapida occhiata:
- sia nella send che nella recv passi l'indirizzo di array_num, che a sua volta è già un puntatore
- nella parte di ricezione array_num non è inizializzato
- argc conta sempre anche il comando usato per lanciare il programma (argv[0]), quindi il numero totale di reali immessi è argc-1