November 2014 - discuss - mailman3.cels.anl.gov

Re: [mpich-discuss] mpich run error
by نازنین 29 Nov '14

29 Nov '14

please give me more detail . you mean to compile mpich-3.0.4 with what flags? ./configure --prefix=/home/nazanin/program_install/mpich-3.0.4 make make install do u mean to add another flag to the command ? thanks nazi On Sat, Nov 29, 2014 at 8:53 AM, Junchao Zhang <jczhang(a)mcs.anl.gov> wrote: > Try to compile your program with mpicc or mpicxx > On Nov 29, 2014 8:15 AM, "نازنین" <nazanin.mirshokraei(a)gmail.com> wrote: > >> thank u for your reply .... can you please tell me what is the reason? it >> is my first time running a mpi .... >> i compiled a program with cpp which tells to use : >> export USE_MPI=on # distributed-memory parallelism >> export USE_MPIF90=on # compile with mpif90 script >> export which_MPI=mpich # compile with MPICH library >> #export which_MPI=mpich2 # compile with MPICH2 library >> # export which_MPI=openmpi # compile with OpenMPI library >> >> #export USE_OpenMP=on # shared-memory parallelism >> >> # export FORT=ifort >> export FORT=gfortran >> #export FORT=pgi >> >> #export USE_DEBUG= # use Fortran debugging flags >> export USE_LARGE=on # activate 64-bit compilation >> #export USE_NETCDF4=on # compile with NetCDF-4 library >> #export USE_PARALLEL_IO=on # Parallel I/O with Netcdf-4/HDF5 >> >> #export USE_MY_LIBS=on # use my library paths below >> >> >> >> and then it gives me a file(coawstM) and the manual tells me to run that >> file in this way : >> mpirun -np 2 ./coawstM swan_only.in >> >> what will be the reason? is it because mpich was not compiled truely or >> the coawstM is not working for a reason? >> thank you >> >> On Sat, Nov 29, 2014 at 6:07 AM, Junchao Zhang <jczhang(a)mcs.anl.gov> >> wrote: >> >>> As indicated by the message, you passed an invalid communicator to MPI_Comm_rank() >>> in your program. >>> >>> --Junchao Zhang >>> >>> On Sat, Nov 29, 2014 at 4:32 AM, نازنین <nazanin.mirshokraei(a)gmail.com> >>> wrote: >>> >>>> hi does any one know about this error : >>>> Fatal error in PMPI_Comm_rank: Invalid communicator, error stack: >>>> PMPI_Comm_rank(108): MPI_Comm_rank(comm=0x0, rank=0x7fff24a405bc) failed >>>> PMPI_Comm_rank(66).: Invalid communicator >>>> Fatal error in PMPI_Comm_rank: Invalid communicator, error stack: >>>> PMPI_Comm_rank(108): MPI_Comm_rank(comm=0x0, rank=0x7fff63e1598c) failed >>>> PMPI_Comm_rank(66).: Invalid communicator >>>> >>>> it is when i run mpich like : >>>> mpirun -np 2 ./coawstM swan_only.in >>>> >>>> i don't know if this is a problem in my executable file (coawstM) or >>>> this is from mpich or any other hardware issues or anything else >>>> thank you for your help >>>> nazi >>>> >>>> _______________________________________________ >>>> discuss mailing list discuss(a)mpich.org >>>> To manage subscription options or unsubscribe: >>>> https://lists.mpich.org/mailman/listinfo/discuss >>>> >>> >>> >>> _______________________________________________ >>> discuss mailing list discuss(a)mpich.org >>> To manage subscription options or unsubscribe: >>> https://lists.mpich.org/mailman/listinfo/discuss >>> >> >> >> _______________________________________________ >> discuss mailing list discuss(a)mpich.org >> To manage subscription options or unsubscribe: >> https://lists.mpich.org/mailman/listinfo/discuss >> > > _______________________________________________ > discuss mailing list discuss(a)mpich.org > To manage subscription options or unsubscribe: > https://lists.mpich.org/mailman/listinfo/discuss > _______________________________________________ discuss mailing list discuss(a)mpich.org To manage subscription options or unsubscribe: https://lists.mpich.org/mailman/listinfo/discuss

1 0

[patch] fix bug in NEWMAD and MXM netmods
by Kuleshov Aleksey 29 Nov '14

29 Nov '14

diff --git a/src/mpid/ch3/channels/nemesis/netmod/mxm/mxm_poll.c b/src/mpid/ch3/channels/nemesis/netmod/mxm/mxm_poll.c index e8bddc3..752a1f9 100644 --- a/src/mpid/ch3/channels/nemesis/netmod/mxm/mxm_poll.c +++ b/src/mpid/ch3/channels/nemesis/netmod/mxm/mxm_poll.c @@ -482,8 +482,8 @@ static int _mxm_process_rdtype(MPID_Request ** rreq_p, MPI_Datatype datatype, *iov_count = n_iov; } else { - int packsize = 0; - MPIR_Pack_size_impl(rreq->dev.user_count, rreq->dev.datatype, (MPI_Aint *) & packsize); + MPI_Aint packsize = 0; + MPIR_Pack_size_impl(rreq->dev.user_count, rreq->dev.datatype, &packsize); rreq->dev.tmpbuf = MPIU_Malloc((size_t) packsize); MPIU_Assert(rreq->dev.tmpbuf); rreq->dev.tmpbuf_sz = packsize; diff --git a/src/mpid/ch3/channels/nemesis/netmod/newmad/newmad_poll.c b/src/mpid/ch3/channels/nemesis/netmod/newmad/newmad_poll.c index 2dba872..5a32515 100644 --- a/src/mpid/ch3/channels/nemesis/netmod/newmad/newmad_poll.c +++ b/src/mpid/ch3/channels/nemesis/netmod/newmad/newmad_poll.c @@ -575,7 +575,7 @@ int MPID_nem_newmad_process_rdtype(MPID_Request **rreq_p, MPID_Datatype * dt_ptr } else { - int packsize = 0; + MPI_Aint packsize = 0; MPIR_Pack_size_impl(rreq->dev.user_count, rreq->dev.datatype, &packsize); rreq->dev.tmpbuf = MPIU_Malloc((size_t) packsize); MPIU_Assert(rreq->dev.tmpbuf); _______________________________________________ discuss mailing list discuss(a)mpich.org To manage subscription options or unsubscribe: https://lists.mpich.org/mailman/listinfo/discuss

1 0

Re: [mpich-discuss] mpich run error
by Junchao Zhang 29 Nov '14

29 Nov '14

Try to compile your program with mpicc or mpicxx On Nov 29, 2014 8:15 AM, "نازنین" <nazanin.mirshokraei(a)gmail.com> wrote: > thank u for your reply .... can you please tell me what is the reason? it > is my first time running a mpi .... > i compiled a program with cpp which tells to use : > export USE_MPI=on # distributed-memory parallelism > export USE_MPIF90=on # compile with mpif90 script > export which_MPI=mpich # compile with MPICH library > #export which_MPI=mpich2 # compile with MPICH2 library > # export which_MPI=openmpi # compile with OpenMPI library > > #export USE_OpenMP=on # shared-memory parallelism > > # export FORT=ifort > export FORT=gfortran > #export FORT=pgi > > #export USE_DEBUG= # use Fortran debugging flags > export USE_LARGE=on # activate 64-bit compilation > #export USE_NETCDF4=on # compile with NetCDF-4 library > #export USE_PARALLEL_IO=on # Parallel I/O with Netcdf-4/HDF5 > > #export USE_MY_LIBS=on # use my library paths below > > > > and then it gives me a file(coawstM) and the manual tells me to run that > file in this way : > mpirun -np 2 ./coawstM swan_only.in > > what will be the reason? is it because mpich was not compiled truely or > the coawstM is not working for a reason? > thank you > > On Sat, Nov 29, 2014 at 6:07 AM, Junchao Zhang <jczhang(a)mcs.anl.gov> > wrote: > >> As indicated by the message, you passed an invalid communicator to MPI_Comm_rank() >> in your program. >> >> --Junchao Zhang >> >> On Sat, Nov 29, 2014 at 4:32 AM, نازنین <nazanin.mirshokraei(a)gmail.com> >> wrote: >> >>> hi does any one know about this error : >>> Fatal error in PMPI_Comm_rank: Invalid communicator, error stack: >>> PMPI_Comm_rank(108): MPI_Comm_rank(comm=0x0, rank=0x7fff24a405bc) failed >>> PMPI_Comm_rank(66).: Invalid communicator >>> Fatal error in PMPI_Comm_rank: Invalid communicator, error stack: >>> PMPI_Comm_rank(108): MPI_Comm_rank(comm=0x0, rank=0x7fff63e1598c) failed >>> PMPI_Comm_rank(66).: Invalid communicator >>> >>> it is when i run mpich like : >>> mpirun -np 2 ./coawstM swan_only.in >>> >>> i don't know if this is a problem in my executable file (coawstM) or >>> this is from mpich or any other hardware issues or anything else >>> thank you for your help >>> nazi >>> >>> _______________________________________________ >>> discuss mailing list discuss(a)mpich.org >>> To manage subscription options or unsubscribe: >>> https://lists.mpich.org/mailman/listinfo/discuss >>> >> >> >> _______________________________________________ >> discuss mailing list discuss(a)mpich.org >> To manage subscription options or unsubscribe: >> https://lists.mpich.org/mailman/listinfo/discuss >> > > > _______________________________________________ > discuss mailing list discuss(a)mpich.org > To manage subscription options or unsubscribe: > https://lists.mpich.org/mailman/listinfo/discuss > _______________________________________________ discuss mailing list discuss(a)mpich.org To manage subscription options or unsubscribe: https://lists.mpich.org/mailman/listinfo/discuss

1 0

Re: [mpich-discuss] mpich run error
by نازنین 29 Nov '14

29 Nov '14

thank u for your reply .... can you please tell me what is the reason? it is my first time running a mpi .... i compiled a program with cpp which tells to use : export USE_MPI=on # distributed-memory parallelism export USE_MPIF90=on # compile with mpif90 script export which_MPI=mpich # compile with MPICH library #export which_MPI=mpich2 # compile with MPICH2 library # export which_MPI=openmpi # compile with OpenMPI library #export USE_OpenMP=on # shared-memory parallelism # export FORT=ifort export FORT=gfortran #export FORT=pgi #export USE_DEBUG= # use Fortran debugging flags export USE_LARGE=on # activate 64-bit compilation #export USE_NETCDF4=on # compile with NetCDF-4 library #export USE_PARALLEL_IO=on # Parallel I/O with Netcdf-4/HDF5 #export USE_MY_LIBS=on # use my library paths below and then it gives me a file(coawstM) and the manual tells me to run that file in this way : mpirun -np 2 ./coawstM swan_only.in what will be the reason? is it because mpich was not compiled truely or the coawstM is not working for a reason? thank you On Sat, Nov 29, 2014 at 6:07 AM, Junchao Zhang <jczhang(a)mcs.anl.gov> wrote: > As indicated by the message, you passed an invalid communicator to MPI_Comm_rank() > in your program. > > --Junchao Zhang > > On Sat, Nov 29, 2014 at 4:32 AM, نازنین <nazanin.mirshokraei(a)gmail.com> > wrote: > >> hi does any one know about this error : >> Fatal error in PMPI_Comm_rank: Invalid communicator, error stack: >> PMPI_Comm_rank(108): MPI_Comm_rank(comm=0x0, rank=0x7fff24a405bc) failed >> PMPI_Comm_rank(66).: Invalid communicator >> Fatal error in PMPI_Comm_rank: Invalid communicator, error stack: >> PMPI_Comm_rank(108): MPI_Comm_rank(comm=0x0, rank=0x7fff63e1598c) failed >> PMPI_Comm_rank(66).: Invalid communicator >> >> it is when i run mpich like : >> mpirun -np 2 ./coawstM swan_only.in >> >> i don't know if this is a problem in my executable file (coawstM) or this >> is from mpich or any other hardware issues or anything else >> thank you for your help >> nazi >> >> _______________________________________________ >> discuss mailing list discuss(a)mpich.org >> To manage subscription options or unsubscribe: >> https://lists.mpich.org/mailman/listinfo/discuss >> > > > _______________________________________________ > discuss mailing list discuss(a)mpich.org > To manage subscription options or unsubscribe: > https://lists.mpich.org/mailman/listinfo/discuss > _______________________________________________ discuss mailing list discuss(a)mpich.org To manage subscription options or unsubscribe: https://lists.mpich.org/mailman/listinfo/discuss

1 0

Re: [mpich-discuss] mpich run error
by Junchao Zhang 29 Nov '14

29 Nov '14

As indicated by the message, you passed an invalid communicator to MPI_Comm_rank() in your program. --Junchao Zhang On Sat, Nov 29, 2014 at 4:32 AM, نازنین <nazanin.mirshokraei(a)gmail.com> wrote: > hi does any one know about this error : > Fatal error in PMPI_Comm_rank: Invalid communicator, error stack: > PMPI_Comm_rank(108): MPI_Comm_rank(comm=0x0, rank=0x7fff24a405bc) failed > PMPI_Comm_rank(66).: Invalid communicator > Fatal error in PMPI_Comm_rank: Invalid communicator, error stack: > PMPI_Comm_rank(108): MPI_Comm_rank(comm=0x0, rank=0x7fff63e1598c) failed > PMPI_Comm_rank(66).: Invalid communicator > > it is when i run mpich like : > mpirun -np 2 ./coawstM swan_only.in > > i don't know if this is a problem in my executable file (coawstM) or this > is from mpich or any other hardware issues or anything else > thank you for your help > nazi > > _______________________________________________ > discuss mailing list discuss(a)mpich.org > To manage subscription options or unsubscribe: > https://lists.mpich.org/mailman/listinfo/discuss > _______________________________________________ discuss mailing list discuss(a)mpich.org To manage subscription options or unsubscribe: https://lists.mpich.org/mailman/listinfo/discuss

1 0

Re: [mpich-discuss] (no subject)
by Junchao Zhang 29 Nov '14

29 Nov '14

Don't know what you meant by "add a piece "B" of code aims to read a .txt file nothing appear". The code you added in B in a normal file read/write. You can easily debug it with gdb. --Junchao Zhang On Fri, Nov 28, 2014 at 6:29 PM, Chafik sanaa <san.chafik(a)gmail.com> wrote: > when I execute the "A" program the result is displayed , but when i trying > to add a piece "B" of code aims to read a .txt file nothing appear: > > program A: > --------------------------------------- > #include <malloc.h> > #include <stdlib.h> > #include <stdio.h> > #include <time.h> > #include "math.h" > #include "mpi.h" > > #define MAX_LIGNE 20 /* au maximum, 10 lignes */ > #define MAX_COL 20 /* au maximum 8 colonnes */ > #define NOM_A_LIRE "fichier1.txt" /*matrice a lire*/ > typedef int Matrice[MAX_LIGNE][MAX_COL]; > int main(int argc, char** argv) > { > int taskid, ntasks; > int ierr, i, itask; > int sendcounts[2048], displs[2048], recvcount; > double **sendbuff, *recvbuff, buffsum, buffsums[2048]; > double inittime, totaltime; > const int nbr_etat = 10; > double tab[nbr_etat]; > > > > for (int i = 0; i < nbr_etat; i++) > tab[i] = i; > > for (int i = 0; i < nbr_etat; i++) > printf("%0.0f ", tab[i]); > printf("\n"); > int nbr_elm[2] = { 4, 6 }; > int dpl[2] = { 0, 4 }; > > MPI_Init(&argc, &argv); > MPI_Comm_rank(MPI_COMM_WORLD, &taskid); > MPI_Comm_size(MPI_COMM_WORLD, &ntasks); > > recvbuff = (double *)malloc(sizeof(double)*nbr_etat); > if (taskid == 0) > { > sendbuff = (double **)malloc(sizeof(double *)*ntasks);// on execute pour > deux proc > sendbuff[0] = (double *)malloc(sizeof(double)*ntasks*nbr_etat); > for (i = 1; i < ntasks; i++) > { > sendbuff[i] = sendbuff[i - 1] + nbr_etat; > } > } > else > { > sendbuff = (double **)malloc(sizeof(double *)* 1); > sendbuff[0] = (double *)malloc(sizeof(double)* 1); > } > > if (taskid == 0){ > > srand((unsigned)time(NULL) + taskid); > for (itask = 0; itask<ntasks; itask++) > { > int k; > displs[itask] = itask*nbr_etat; > int s = dpl[itask]; > sendcounts[itask] = nbr_elm[itask]; > > > for (i = 0; i<sendcounts[itask]; i++) > { > k = i + s; > sendbuff[itask][i] = tab[k]; > printf("+ %0.0f ", sendbuff[itask][i]); > } > printf("\n"); > } > } > > recvcount = nbr_elm[taskid]; > > inittime = MPI_Wtime(); > > ierr = MPI_Scatterv(sendbuff[0], sendcounts, displs, MPI_DOUBLE, > recvbuff, recvcount, MPI_DOUBLE, > 0, MPI_COMM_WORLD); > > totaltime = MPI_Wtime() - inittime; > > printf("\n >>>> \n"); > buffsum = 0.0; > printf("\n > %d < \n", taskid); > for (i = 0; i<recvcount; i++) > { > buffsum = buffsum + recvbuff[i]; > printf("* %0.0f ", recvbuff[i]); > } > printf("\n"); > printf("(%d) %0.3f ", taskid, buffsum); > > ierr = MPI_Gather(&buffsum, 1, MPI_DOUBLE, > buffsums, 1, MPI_DOUBLE, > 0, MPI_COMM_WORLD); > if (taskid == 0) > { > printf("\n"); > printf("##########################################################\n\n"); > printf(" --> APReS COMMUNICATION <-- \n\n"); > for (itask = 0; itask<ntasks; itask++) > { > printf("Processus %d : Somme du vecteur reçu : %0.0f\n", > itask, buffsums[itask]); > } > printf("\n"); > printf("##########################################################\n\n"); > printf(" Temps total de communication : %f secondes\n\n", totaltime); > printf("##########################################################\n\n"); > } > > /* Libération de la mémoire */ > if (taskid == 0) > { > free(sendbuff[0]); > free(sendbuff); > } > else > { > free(sendbuff[0]); > free(sendbuff); > free(recvbuff); > } > > /* Finalisation de MPI */ > MPI_Finalize(); > } > > program B: > ------------------------- > #include <malloc.h> > #include <stdlib.h> > #include <stdio.h> > #include <time.h> > #include "math.h" > #include "mpi.h" > > #define MAX_LIGNE 20 /* au maximum, 10 lignes */ > #define MAX_COL 20 /* au maximum 8 colonnes */ > #define NOM_A_LIRE "fichier1.txt" /*matrice a lire*/ > typedef int Matrice[MAX_LIGNE][MAX_COL]; > int main(int argc, char** argv) > { > int taskid, ntasks; > int ierr, i, itask; > int sendcounts[2048], displs[2048], recvcount; > double **sendbuff, *recvbuff, buffsum, buffsums[2048]; > double inittime, totaltime; > const int nbr_etat = 10; > double tab[nbr_etat]; > > > > for (int i = 0; i < nbr_etat; i++) > tab[i] = i; > > for (int i = 0; i < nbr_etat; i++) > printf("%0.0f ", tab[i]); > printf("\n"); > int nbr_elm[2] = { 4, 6 }; > int dpl[2] = { 0, 4 }; > > MPI_Init(&argc, &argv); > MPI_Comm_rank(MPI_COMM_WORLD, &taskid); > MPI_Comm_size(MPI_COMM_WORLD, &ntasks); > > recvbuff = (double *)malloc(sizeof(double)*nbr_etat); > if (taskid == 0) > { > sendbuff = (double **)malloc(sizeof(double *)*ntasks);// on execute pour > deux proc > sendbuff[0] = (double *)malloc(sizeof(double)*ntasks*nbr_etat); > for (i = 1; i < ntasks; i++) > { > sendbuff[i] = sendbuff[i - 1] + nbr_etat; > } > } > else > { > sendbuff = (double **)malloc(sizeof(double *)* 1); > sendbuff[0] = (double *)malloc(sizeof(double)* 1); > } > > if (taskid == 0){ > > srand((unsigned)time(NULL) + taskid); > for (itask = 0; itask<ntasks; itask++) > { > int k; > displs[itask] = itask*nbr_etat; > int s = dpl[itask]; > sendcounts[itask] = nbr_elm[itask]; > > > for (i = 0; i<sendcounts[itask]; i++) > { > k = i + s; > sendbuff[itask][i] = tab[k]; > printf("+ %0.0f ", sendbuff[itask][i]); > } > printf("\n"); > } > } > > recvcount = nbr_elm[taskid]; > > inittime = MPI_Wtime(); > > ierr = MPI_Scatterv(sendbuff[0], sendcounts, displs, MPI_DOUBLE, > recvbuff, recvcount, MPI_DOUBLE, > 0, MPI_COMM_WORLD); > > totaltime = MPI_Wtime() - inittime; > > printf("\n >>>> \n"); > buffsum = 0.0; > printf("\n > %d < \n", taskid); > for (i = 0; i<recvcount; i++) > { > buffsum = buffsum + recvbuff[i]; > printf("* %0.0f ", recvbuff[i]); > } > printf("\n"); > printf("(%d) %0.3f ", taskid, buffsum); > > ierr = MPI_Gather(&buffsum, 1, MPI_DOUBLE, > buffsums, 1, MPI_DOUBLE, > 0, MPI_COMM_WORLD); > if (taskid == 0) > { > printf("\n"); > printf("##########################################################\n\n"); > printf(" --> APReS COMMUNICATION <-- \n\n"); > for (itask = 0; itask<ntasks; itask++) > { > printf("Processus %d : Somme du vecteur reçu : %0.0f\n", > itask, buffsums[itask]); > } > printf("\n"); > printf("##########################################################\n\n"); > printf(" Temps total de communication : %f secondes\n\n", totaltime); > printf("##########################################################\n\n"); > } > > /*===============================================================*/ > /* Libération de la mémoire */ > if (taskid == 0) > { > free(sendbuff[0]); > free(sendbuff); > } > else > { > free(sendbuff[0]); > free(sendbuff); > free(recvbuff); > } > > /*===============================================================*/ > /* Finalisation de MPI */ > MPI_Finalize(); > Matrice A; > int m, n; > /* ouvrir le fichier Matrices.dta sur le disque réseau pour la lecture */ > FILE * aLire = fopen(NOM_A_LIRE, "r"); > if (aLire == NULL) > { > printf("Le fichier n'existe pas"); > } > else > { > fscanf(aLire, "%d%d\n", &m, &n); > ////void lire(FILE * aLire, Matrice mat, int nbLigne, int nbCol) > for (int i = 0; i < m; i++) { > for (int j = 0; j < n; j++) > fscanf(aLire, "%d", &A[i][j]); > fscanf(aLire, "\n"); > } > ////void afficher(Matrice mat, char nom[], int nbLigne, int nbCol) > > printf("\nContenu de la matrice de %d ligne(s) et %d colonne(s) :\n", m, > n); > for (int i = 0; i < m; i++) { > for (int j = 0; j < n; j++) > printf("%d ", A[i][j]); > printf("\n"); > } > printf("\n"); > > > } > fclose(aLire); > > } > > _______________________________________________ discuss mailing list discuss(a)mpich.org To manage subscription options or unsubscribe: https://lists.mpich.org/mailman/listinfo/discuss

1 0

mpich run error
by نازنین 29 Nov '14

29 Nov '14

hi does any one know about this error : Fatal error in PMPI_Comm_rank: Invalid communicator, error stack: PMPI_Comm_rank(108): MPI_Comm_rank(comm=0x0, rank=0x7fff24a405bc) failed PMPI_Comm_rank(66).: Invalid communicator Fatal error in PMPI_Comm_rank: Invalid communicator, error stack: PMPI_Comm_rank(108): MPI_Comm_rank(comm=0x0, rank=0x7fff63e1598c) failed PMPI_Comm_rank(66).: Invalid communicator it is when i run mpich like : mpirun -np 2 ./coawstM swan_only.in i don't know if this is a problem in my executable file (coawstM) or this is from mpich or any other hardware issues or anything else thank you for your help nazi _______________________________________________ discuss mailing list discuss(a)mpich.org To manage subscription options or unsubscribe: https://lists.mpich.org/mailman/listinfo/discuss

1 0

(no subject)
by Chafik sanaa 28 Nov '14

28 Nov '14

when I execute the "A" program the result is displayed , but when i trying to add a piece "B" of code aims to read a .txt file nothing appear: program A: --------------------------------------- #include <malloc.h> #include <stdlib.h> #include <stdio.h> #include <time.h> #include "math.h" #include "mpi.h" #define MAX_LIGNE 20 /* au maximum, 10 lignes */ #define MAX_COL 20 /* au maximum 8 colonnes */ #define NOM_A_LIRE "fichier1.txt" /*matrice a lire*/ typedef int Matrice[MAX_LIGNE][MAX_COL]; int main(int argc, char** argv) { int taskid, ntasks; int ierr, i, itask; int sendcounts[2048], displs[2048], recvcount; double **sendbuff, *recvbuff, buffsum, buffsums[2048]; double inittime, totaltime; const int nbr_etat = 10; double tab[nbr_etat]; for (int i = 0; i < nbr_etat; i++) tab[i] = i; for (int i = 0; i < nbr_etat; i++) printf("%0.0f ", tab[i]); printf("\n"); int nbr_elm[2] = { 4, 6 }; int dpl[2] = { 0, 4 }; MPI_Init(&argc, &argv); MPI_Comm_rank(MPI_COMM_WORLD, &taskid); MPI_Comm_size(MPI_COMM_WORLD, &ntasks); recvbuff = (double *)malloc(sizeof(double)*nbr_etat); if (taskid == 0) { sendbuff = (double **)malloc(sizeof(double *)*ntasks);// on execute pour deux proc sendbuff[0] = (double *)malloc(sizeof(double)*ntasks*nbr_etat); for (i = 1; i < ntasks; i++) { sendbuff[i] = sendbuff[i - 1] + nbr_etat; } } else { sendbuff = (double **)malloc(sizeof(double *)* 1); sendbuff[0] = (double *)malloc(sizeof(double)* 1); } if (taskid == 0){ srand((unsigned)time(NULL) + taskid); for (itask = 0; itask<ntasks; itask++) { int k; displs[itask] = itask*nbr_etat; int s = dpl[itask]; sendcounts[itask] = nbr_elm[itask]; for (i = 0; i<sendcounts[itask]; i++) { k = i + s; sendbuff[itask][i] = tab[k]; printf("+ %0.0f ", sendbuff[itask][i]); } printf("\n"); } } recvcount = nbr_elm[taskid]; inittime = MPI_Wtime(); ierr = MPI_Scatterv(sendbuff[0], sendcounts, displs, MPI_DOUBLE, recvbuff, recvcount, MPI_DOUBLE, 0, MPI_COMM_WORLD); totaltime = MPI_Wtime() - inittime; printf("\n >>>> \n"); buffsum = 0.0; printf("\n > %d < \n", taskid); for (i = 0; i<recvcount; i++) { buffsum = buffsum + recvbuff[i]; printf("* %0.0f ", recvbuff[i]); } printf("\n"); printf("(%d) %0.3f ", taskid, buffsum); ierr = MPI_Gather(&buffsum, 1, MPI_DOUBLE, buffsums, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); if (taskid == 0) { printf("\n"); printf("##########################################################\n\n"); printf(" --> APReS COMMUNICATION <-- \n\n"); for (itask = 0; itask<ntasks; itask++) { printf("Processus %d : Somme du vecteur reçu : %0.0f\n", itask, buffsums[itask]); } printf("\n"); printf("##########################################################\n\n"); printf(" Temps total de communication : %f secondes\n\n", totaltime); printf("##########################################################\n\n"); } /* Libération de la mémoire */ if (taskid == 0) { free(sendbuff[0]); free(sendbuff); } else { free(sendbuff[0]); free(sendbuff); free(recvbuff); } /* Finalisation de MPI */ MPI_Finalize(); } program B: ------------------------- #include <malloc.h> #include <stdlib.h> #include <stdio.h> #include <time.h> #include "math.h" #include "mpi.h" #define MAX_LIGNE 20 /* au maximum, 10 lignes */ #define MAX_COL 20 /* au maximum 8 colonnes */ #define NOM_A_LIRE "fichier1.txt" /*matrice a lire*/ typedef int Matrice[MAX_LIGNE][MAX_COL]; int main(int argc, char** argv) { int taskid, ntasks; int ierr, i, itask; int sendcounts[2048], displs[2048], recvcount; double **sendbuff, *recvbuff, buffsum, buffsums[2048]; double inittime, totaltime; const int nbr_etat = 10; double tab[nbr_etat]; for (int i = 0; i < nbr_etat; i++) tab[i] = i; for (int i = 0; i < nbr_etat; i++) printf("%0.0f ", tab[i]); printf("\n"); int nbr_elm[2] = { 4, 6 }; int dpl[2] = { 0, 4 }; MPI_Init(&argc, &argv); MPI_Comm_rank(MPI_COMM_WORLD, &taskid); MPI_Comm_size(MPI_COMM_WORLD, &ntasks); recvbuff = (double *)malloc(sizeof(double)*nbr_etat); if (taskid == 0) { sendbuff = (double **)malloc(sizeof(double *)*ntasks);// on execute pour deux proc sendbuff[0] = (double *)malloc(sizeof(double)*ntasks*nbr_etat); for (i = 1; i < ntasks; i++) { sendbuff[i] = sendbuff[i - 1] + nbr_etat; } } else { sendbuff = (double **)malloc(sizeof(double *)* 1); sendbuff[0] = (double *)malloc(sizeof(double)* 1); } if (taskid == 0){ srand((unsigned)time(NULL) + taskid); for (itask = 0; itask<ntasks; itask++) { int k; displs[itask] = itask*nbr_etat; int s = dpl[itask]; sendcounts[itask] = nbr_elm[itask]; for (i = 0; i<sendcounts[itask]; i++) { k = i + s; sendbuff[itask][i] = tab[k]; printf("+ %0.0f ", sendbuff[itask][i]); } printf("\n"); } } recvcount = nbr_elm[taskid]; inittime = MPI_Wtime(); ierr = MPI_Scatterv(sendbuff[0], sendcounts, displs, MPI_DOUBLE, recvbuff, recvcount, MPI_DOUBLE, 0, MPI_COMM_WORLD); totaltime = MPI_Wtime() - inittime; printf("\n >>>> \n"); buffsum = 0.0; printf("\n > %d < \n", taskid); for (i = 0; i<recvcount; i++) { buffsum = buffsum + recvbuff[i]; printf("* %0.0f ", recvbuff[i]); } printf("\n"); printf("(%d) %0.3f ", taskid, buffsum); ierr = MPI_Gather(&buffsum, 1, MPI_DOUBLE, buffsums, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); if (taskid == 0) { printf("\n"); printf("##########################################################\n\n"); printf(" --> APReS COMMUNICATION <-- \n\n"); for (itask = 0; itask<ntasks; itask++) { printf("Processus %d : Somme du vecteur reçu : %0.0f\n", itask, buffsums[itask]); } printf("\n"); printf("##########################################################\n\n"); printf(" Temps total de communication : %f secondes\n\n", totaltime); printf("##########################################################\n\n"); } /*===============================================================*/ /* Libération de la mémoire */ if (taskid == 0) { free(sendbuff[0]); free(sendbuff); } else { free(sendbuff[0]); free(sendbuff); free(recvbuff); } /*===============================================================*/ /* Finalisation de MPI */ MPI_Finalize(); Matrice A; int m, n; /* ouvrir le fichier Matrices.dta sur le disque réseau pour la lecture */ FILE * aLire = fopen(NOM_A_LIRE, "r"); if (aLire == NULL) { printf("Le fichier n'existe pas"); } else { fscanf(aLire, "%d%d\n", &m, &n); ////void lire(FILE * aLire, Matrice mat, int nbLigne, int nbCol) for (int i = 0; i < m; i++) { for (int j = 0; j < n; j++) fscanf(aLire, "%d", &A[i][j]); fscanf(aLire, "\n"); } ////void afficher(Matrice mat, char nom[], int nbLigne, int nbCol) printf("\nContenu de la matrice de %d ligne(s) et %d colonne(s) :\n", m, n); for (int i = 0; i < m; i++) { for (int j = 0; j < n; j++) printf("%d ", A[i][j]); printf("\n"); } printf("\n"); } fclose(aLire); } _______________________________________________ discuss mailing list discuss(a)mpich.org To manage subscription options or unsubscribe: https://lists.mpich.org/mailman/listinfo/discuss

1 0

Re: [mpich-discuss] the SCATTERV function
by Chafik sanaa 28 Nov '14

28 Nov '14

I was able to correct, thanks #include <malloc.h> #include <stdlib.h> #include <stdio.h> #include <time.h> #include "math.h" #include "mpi.h" int main(int argc, char** argv) { int taskid, ntasks; int ierr, i, itask; int sendcounts[2048], displs[2048], recvcount; double **sendbuff, *recvbuff, buffsum, buffsums[2048]; double inittime, totaltime; const int nbr_etat = 10; double tab[nbr_etat]; for (int i = 0; i < nbr_etat; i++) tab[i] = i; for (int i = 0; i < nbr_etat; i++) printf("%0.0f ", tab[i]); printf("\n"); int nbr_elm[2] = { 4, 6 }; int dpl[2] = { 0, 4 }; MPI_Init(&argc, &argv); MPI_Comm_rank(MPI_COMM_WORLD, &taskid); MPI_Comm_size(MPI_COMM_WORLD, &ntasks); recvbuff = (double *)malloc(sizeof(double)*nbr_etat); if (taskid == 0) { sendbuff = (double **)malloc(sizeof(double *)*ntasks);// on execute pour deux proc sendbuff[0] = (double *)malloc(sizeof(double)*ntasks*nbr_etat); for (i = 1; i < ntasks; i++) { sendbuff[i] = sendbuff[i - 1] + nbr_etat; } } else { sendbuff = (double **)malloc(sizeof(double *)* 1); sendbuff[0] = (double *)malloc(sizeof(double)* 1); } if (taskid == 0){ srand((unsigned)time(NULL) + taskid); for (itask = 0; itask<ntasks; itask++) { int k; displs[itask] = itask*nbr_etat; int s = dpl[itask]; sendcounts[itask] = nbr_elm[itask]; for (i = 0; i<sendcounts[itask]; i++) { k = i + s; sendbuff[itask][i] = tab[k]; printf("+ %0.0f ", sendbuff[itask][i]); } printf("\n"); } } recvcount = nbr_elm[taskid]; inittime = MPI_Wtime(); ierr = MPI_Scatterv(sendbuff[0], sendcounts, displs, MPI_DOUBLE, recvbuff, recvcount, MPI_DOUBLE, 0, MPI_COMM_WORLD); totaltime = MPI_Wtime() - inittime; printf("\n >>>> \n"); buffsum = 0.0; printf("\n > %d < \n", taskid); for (i = 0; i<recvcount; i++) { printf("* %0.0f ", recvbuff[i]); } printf("\n"); if (taskid == 0){ free(sendbuff[0]); free(sendbuff); } else{ free(sendbuff[0]); free(sendbuff); free(recvbuff); } /*===============================================================*/ /* Finalisation de MPI */ MPI_Finalize(); } 2014-11-28 20:52 GMT+01:00 Chafik sanaa <san.chafik(a)gmail.com>: > Hi > when i execute this program (I use two processes) i have a error in the > display part : > * RESULT OF EXECUTION: > 0 1 2 3 4 5 6 7 8 9 > >> (1,6) << > > >>>> > > > 1 < > * -6277438562204192500000000000000000000000000000000000000000000000000 * > -6277438562204192500000000000000000000000000000000000000000000000000 * > -6277438562204192500000000000000000000000000000000000000000000000000 * > -6277438562204192500000000000000000000000000000000000000000000000000 * > -6277438562204192500000000000000000000000000000000000000000000000000 * > -6277438562204192500000000000000000000000000000000000000000000000000 > 0 1 2 3 4 5 6 7 8 9 > 0 1 2 3 > 4 5 6 7 8 9 > >> (0,4) << > > >>>> > > > 0 < > * 0 * 1 * 2 * 3 > * PROGRAM > #include <malloc.h> > #include <stdlib.h> > #include <stdio.h> > #include <time.h> > #include "math.h" > #include "mpi.h" > > int main(int argc, char** argv) > { > int taskid, ntasks; > int ierr, i, itask; > int sendcounts[2048], displs[2048], recvcount; > double **sendbuff, *recvbuff, buffsum, buffsums[2048]; > double inittime, totaltime; > const int nbr_etat = 10; > double tab[nbr_etat]; > > for (int i = 0; i < nbr_etat; i++) > tab[i] = i; > > for (int i = 0; i < nbr_etat; i++) > printf("%0.0f ", tab[i]); > printf("\n"); > int nbr_elm[2] = { 4, 6 }; > int dpl[2] = { 0, 4 }; > > MPI_Init(&argc, &argv); > MPI_Comm_rank(MPI_COMM_WORLD, &taskid); > MPI_Comm_size(MPI_COMM_WORLD, &ntasks); > > > recvbuff = (double *)malloc(sizeof(double)*nbr_etat); > if (taskid == 0) > { > sendbuff = (double **)malloc(sizeof(double *)*ntasks);// on execute pour > deux proc > sendbuff[0] = (double *)malloc(sizeof(double)*ntasks*nbr_etat); > for (i = 1; i < ntasks; i++) > { > sendbuff[i] = sendbuff[i - 1] + nbr_etat; > } > } > else > { > sendbuff = (double **)malloc(sizeof(double *)* 1); > sendbuff[0] = (double *)malloc(sizeof(double)* 1); > } > > if (taskid == 0){ > > > > > srand((unsigned)time(NULL) + taskid); > for (itask = 0; itask<ntasks; itask++) > { > int k; > displs[itask] = dpl[itask]; > int s = displs[itask]; > sendcounts[itask] = nbr_elm[itask]; > > for (i = 0; i<sendcounts[itask]; i++) > { > k = i + s; > sendbuff[itask][i] = tab[k]; > printf("%0.0f ", sendbuff[itask][i]); > } > printf("\n"); > } > > > } > > > recvcount = nbr_elm[taskid]; > > inittime = MPI_Wtime(); > > ierr = MPI_Scatterv(sendbuff[0], sendcounts, displs, MPI_DOUBLE, > recvbuff, recvcount, MPI_DOUBLE, > 0, MPI_COMM_WORLD); > > totaltime = MPI_Wtime() - inittime; > > printf("\n >>>> \n"); > buffsum = 0.0; > printf("\n > %d < \n",taskid); > for (i = 0; i<recvcount; i++) > { > printf("* %0.0f ", recvbuff[i]); > } > > printf("\n"); > if (taskid == 0){ > free(sendbuff[0]); > free(sendbuff); > } > else{ > free(sendbuff[0]); > free(sendbuff); > free(recvbuff); > } > > /*===============================================================*/ > /* Finalisation de MPI */ > MPI_Finalize(); > > } > > > 2014-11-28 20:53 GMT+01:00 Chafik sanaa <san.chafik(a)gmail.com>: > Hi > when i execute this program (I use two processes) i have a error in the > display part : > * RESULT OF EXECUTION: > 0 1 2 3 4 5 6 7 8 9 > >> (1,6) << > > >>>> > > > 1 < > * -6277438562204192500000000000000000000000000000000000000000000000000 * > -6277438562204192500000000000000000000000000000000000000000000000000 * > -6277438562204192500000000000000000000000000000000000000000000000000 * > -6277438562204192500000000000000000000000000000000000000000000000000 * > -6277438562204192500000000000000000000000000000000000000000000000000 * > -6277438562204192500000000000000000000000000000000000000000000000000 > 0 1 2 3 4 5 6 7 8 9 > 0 1 2 3 > 4 5 6 7 8 9 > >> (0,4) << > > >>>> > > > 0 < > * 0 * 1 * 2 * 3 > * PROGRAM > #include <malloc.h> > #include <stdlib.h> > #include <stdio.h> > #include <time.h> > #include "math.h" > #include "mpi.h" > > int main(int argc, char** argv) > { > int taskid, ntasks; > int ierr, i, itask; > int sendcounts[2048], displs[2048], recvcount; > double **sendbuff, *recvbuff, buffsum, buffsums[2048]; > double inittime, totaltime; > const int nbr_etat = 10; > double tab[nbr_etat]; > > for (int i = 0; i < nbr_etat; i++) > tab[i] = i; > > for (int i = 0; i < nbr_etat; i++) > printf("%0.0f ", tab[i]); > printf("\n"); > int nbr_elm[2] = { 4, 6 }; > int dpl[2] = { 0, 4 }; > > MPI_Init(&argc, &argv); > MPI_Comm_rank(MPI_COMM_WORLD, &taskid); > MPI_Comm_size(MPI_COMM_WORLD, &ntasks); > > > recvbuff = (double *)malloc(sizeof(double)*nbr_etat); > if (taskid == 0) > { > sendbuff = (double **)malloc(sizeof(double *)*ntasks);// on execute pour > deux proc > sendbuff[0] = (double *)malloc(sizeof(double)*ntasks*nbr_etat); > for (i = 1; i < ntasks; i++) > { > sendbuff[i] = sendbuff[i - 1] + nbr_etat; > } > } > else > { > sendbuff = (double **)malloc(sizeof(double *)* 1); > sendbuff[0] = (double *)malloc(sizeof(double)* 1); > } > > if (taskid == 0){ > > > > > srand((unsigned)time(NULL) + taskid); > for (itask = 0; itask<ntasks; itask++) > { > int k; > displs[itask] = dpl[itask]; > int s = displs[itask]; > sendcounts[itask] = nbr_elm[itask]; > > for (i = 0; i<sendcounts[itask]; i++) > { > k = i + s; > sendbuff[itask][i] = tab[k]; > printf("%0.0f ", sendbuff[itask][i]); > } > printf("\n"); > } > > > } > > > recvcount = nbr_elm[taskid]; > > inittime = MPI_Wtime(); > > ierr = MPI_Scatterv(sendbuff[0], sendcounts, displs, MPI_DOUBLE, > recvbuff, recvcount, MPI_DOUBLE, > 0, MPI_COMM_WORLD); > > totaltime = MPI_Wtime() - inittime; > > printf("\n >>>> \n"); > buffsum = 0.0; > printf("\n > %d < \n",taskid); > for (i = 0; i<recvcount; i++) > { > printf("* %0.0f ", recvbuff[i]); > } > > printf("\n"); > if (taskid == 0){ > free(sendbuff[0]); > free(sendbuff); > } > else{ > free(sendbuff[0]); > free(sendbuff); > free(recvbuff); > } > > /*===============================================================*/ > /* Finalisation de MPI */ > MPI_Finalize(); > > } > > _______________________________________________ discuss mailing list discuss(a)mpich.org To manage subscription options or unsubscribe: https://lists.mpich.org/mailman/listinfo/discuss

1 0

the SCATTERV function
by Chafik sanaa 28 Nov '14

28 Nov '14

Hi when i execute this program (I use two processes) i have a error in the display part : * RESULT OF EXECUTION: 0 1 2 3 4 5 6 7 8 9 >> (1,6) << >>>> > 1 < * -6277438562204192500000000000000000000000000000000000000000000000000 * -6277438562204192500000000000000000000000000000000000000000000000000 * -6277438562204192500000000000000000000000000000000000000000000000000 * -6277438562204192500000000000000000000000000000000000000000000000000 * -6277438562204192500000000000000000000000000000000000000000000000000 * -6277438562204192500000000000000000000000000000000000000000000000000 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 >> (0,4) << >>>> > 0 < * 0 * 1 * 2 * 3 * PROGRAM #include <malloc.h> #include <stdlib.h> #include <stdio.h> #include <time.h> #include "math.h" #include "mpi.h" int main(int argc, char** argv) { int taskid, ntasks; int ierr, i, itask; int sendcounts[2048], displs[2048], recvcount; double **sendbuff, *recvbuff, buffsum, buffsums[2048]; double inittime, totaltime; const int nbr_etat = 10; double tab[nbr_etat]; for (int i = 0; i < nbr_etat; i++) tab[i] = i; for (int i = 0; i < nbr_etat; i++) printf("%0.0f ", tab[i]); printf("\n"); int nbr_elm[2] = { 4, 6 }; int dpl[2] = { 0, 4 }; MPI_Init(&argc, &argv); MPI_Comm_rank(MPI_COMM_WORLD, &taskid); MPI_Comm_size(MPI_COMM_WORLD, &ntasks); recvbuff = (double *)malloc(sizeof(double)*nbr_etat); if (taskid == 0) { sendbuff = (double **)malloc(sizeof(double *)*ntasks);// on execute pour deux proc sendbuff[0] = (double *)malloc(sizeof(double)*ntasks*nbr_etat); for (i = 1; i < ntasks; i++) { sendbuff[i] = sendbuff[i - 1] + nbr_etat; } } else { sendbuff = (double **)malloc(sizeof(double *)* 1); sendbuff[0] = (double *)malloc(sizeof(double)* 1); } if (taskid == 0){ srand((unsigned)time(NULL) + taskid); for (itask = 0; itask<ntasks; itask++) { int k; displs[itask] = dpl[itask]; int s = displs[itask]; sendcounts[itask] = nbr_elm[itask]; for (i = 0; i<sendcounts[itask]; i++) { k = i + s; sendbuff[itask][i] = tab[k]; printf("%0.0f ", sendbuff[itask][i]); } printf("\n"); } } recvcount = nbr_elm[taskid]; inittime = MPI_Wtime(); ierr = MPI_Scatterv(sendbuff[0], sendcounts, displs, MPI_DOUBLE, recvbuff, recvcount, MPI_DOUBLE, 0, MPI_COMM_WORLD); totaltime = MPI_Wtime() - inittime; printf("\n >>>> \n"); buffsum = 0.0; printf("\n > %d < \n",taskid); for (i = 0; i<recvcount; i++) { printf("* %0.0f ", recvbuff[i]); } printf("\n"); if (taskid == 0){ free(sendbuff[0]); free(sendbuff); } else{ free(sendbuff[0]); free(sendbuff); free(recvbuff); } /*===============================================================*/ /* Finalisation de MPI */ MPI_Finalize(); } _______________________________________________ discuss mailing list discuss(a)mpich.org To manage subscription options or unsubscribe: https://lists.mpich.org/mailman/listinfo/discuss

1 0