MPI

Živijo,

Sem gledal po forumu, pa se mi zdi da podobne teme še ni bilo odprte. Skratka, zanima me če ima morda kdo kakšne izkušnje z implementiranjem MPI (wiki)?

Vzpostavil bi rad "emulirano" okolje, tako da bi zadevo osvojil na navadnem enoprocesorskem računalu, potem pa bi iz tega gradil dalje. Predvidevam da ljudje, ki razvijajo programe za izvajanje na clustrih tudi ne počnejo tega v "produkcijskem" okolju.

Skratka na Debianu sem namestil MPICH knjižnice ter prevedel demo na wikipedi-i. Ko pa zadevo poženem, se zadeva ustavi ob klicu MPI_Init(). Google mi ni znal postreči s pravimi odgovori kako naprej, možno pa je, da ga nisem znal pravilno vprašati. Se morda kdo ukvarja s tem in bi mi znal svetovati kako vzpostaviti testno okolje na enoprocesorskem računalu?

Spodaj še prilagam kodo z Wikipedie, na kateri sem zadevo testiral.

 /*
  "Hello World" Type MPI Test Program
 */
 #include <mpi.h>
 #include <stdio.h>
 #include <string.h>
 
 #define BUFSIZE 128
 #define TAG 0
 
 int main(int argc, char *argv[])
 {
   char idstr[32];
   char buff[BUFSIZE];
   int numprocs;
   int myid;
   int i;
   MPI_Status stat; 
 
   MPI_Init(&argc,&argv); /* all MPI programs start with MPI_Init; all 'N' processes exist thereafter */
   MPI_Comm_size(MPI_COMM_WORLD,&numprocs); /* find out how big the SPMD world is */
   MPI_Comm_rank(MPI_COMM_WORLD,&myid); /* and this processes' rank is */
 
   /* At this point, all the programs are running equivalently, the rank is used to
      distinguish the roles of the programs in the SPMD model, with rank 0 often used
      specially... */
   if(myid == 0)
   {
     printf("%d: We have %d processors\n", myid, numprocs);
     for(i=1;i<numprocs;i++)
     {
       sprintf(buff, "Hello %d! ", i);
       MPI_Send(buff, BUFSIZE, MPI_CHAR, i, TAG, MPI_COMM_WORLD);
     }
     for(i=1;i<numprocs;i++)
     {
       MPI_Recv(buff, BUFSIZE, MPI_CHAR, i, TAG, MPI_COMM_WORLD, &stat);
       printf("%d: %s\n", myid, buff);
     }
   }
   else
   {
     /* receive from rank 0: */
     MPI_Recv(buff, BUFSIZE, MPI_CHAR, 0, TAG, MPI_COMM_WORLD, &stat);
     sprintf(idstr, "Processor %d ", myid);
     strcat(buff, idstr);
     strcat(buff, "reporting for duty\n");
     /* send to rank 0: */
     MPI_Send(buff, BUFSIZE, MPI_CHAR, 0, TAG, MPI_COMM_WORLD);
   }
 
   MPI_Finalize(); /* MPI Programs end with MPI Finalize; this is a weak synchronization point */
   return 0;
 }

Edit: Popravil formatiranje kode.

Mi je uspelo zadevo pripravit do izvajanja. Namestil sem tako LAM kot tudi OpenMPI; mimogrede Debian ima to odlično podprto (za ostale distrote ne vem) -- sedaj pa sledi primerjava in testiranje.

@jmakov
Hvala za odgovor in povezavo glede izvajanja na GPU-jih. Vsekakor zanimivo branje. Po guglanju dalje sem ugotovil, da so nekatere algoritme (gledal sem tiste ki so zanimivi zame) uspeli pohitriti med 15 do 100-krat, kar je čista zmaga.

LPr