Dev Faq
8
votes

Problèmes OpenMPI MPI_Barrier

J'ai des problèmes de synchronisation à l'aide de la mise en œuvre OpenMPI de MPI_Barrier: 

int rank;
int nprocs;

int rc = MPI_Init(&argc, &argv);

if(rc != MPI_SUCCESS) {
    fprintf(stderr, "Unable to set up MPI");
    MPI_Abort(MPI_COMM_WORLD, rc);
}

MPI_Comm_size(MPI_COMM_WORLD, &nprocs);
MPI_Comm_rank(MPI_COMM_WORLD, &rank);


printf("P%d\n", rank);
fflush(stdout);

MPI_Barrier(MPI_COMM_WORLD);

printf("P%d again\n", rank);

MPI_Finalize();

0 commentaires

4 Réponses :

14
votes

L'ordre dans lequel vos lignes d'impression apparaissent sur votre terminal n'est pas nécessairement l'ordre dans lequel les choses sont imprimées. Vous utilisez une ressource partagée ( stdout ) pour cela afin qu'il doit toujours être un problème de commande. (Et fflush ne vous aide pas ici, stdout est tamponné de la ligne tamponnée de toute façon.)

Vous pouvez essayer de préfixer votre sortie avec un horodatage et de sauvegarder tout cela à différents fichiers, un par processus MPI.

puis pour inspecter votre journal, vous pouvez fusionner les deux fichiers ensemble et trier en fonction de l'horodatage.

Votre problème devrait disparaître, alors.

4 commentaires

S'appuyant sur des horodatés peut ne pas être idéal si les processus MPI fonctionnent sur différents nœuds, sauf si vous pouvez garantir que les horloges sont synchronisées.

@Shawn: Il y a MPI_WTime () pour cela.

@suszterpatt: MPI_WTime () n'est généralement pas une horloge globale / synchronisée! (C'est seulement si mpi_wtime_is_global est défini et vrai)

AFAIK, MPI_WTIME IN OPENMPI n'est pas synchronisé.

12
votes

Il n'y a rien de mal avec MPI_Barrier ().

comme Jens mentionné , la raison pour laquelle vous ne voyez pas la sortie que vous attendez, c'est que STDOUT est tamponnée sur chaque processus. Il n'y a aucune garantie que les impressions à partir de plusieurs processus seront affichées sur le processus d'appel dans l'ordre. (Si STDOUT de chaque processus est transféré au processus principal d'impression en temps réel, cela conduira à de nombreuses communications inutiles!)

Si vous voulez vous convaincre que la barrière fonctionne, vous pouvez essayer écrire dans un fichier à la place. Avoir de multiples processus écrivant dans un seul fichier peut entraîner des complications supplémentaires, vous pouvez donc avoir chaque processus d'écriture sur un fichier, puis après la barrière, échangez les fichiers qu'ils écrivent à. Par exemple: xxx

Exemple de mise en oeuvre: xxx

Après avoir exécuté le code, vous devez obtenir les résultats suivants: xxx

pour tous les fichiers, les relevés "After Barrière" apparaîtront toujours plus tard.

0 commentaires
3
votes

La commande de sortie n'est pas garantie dans les programmes MPI.

Ceci n'est pas lié à MPI_Barrier du tout.

En outre, je ne passerais pas trop de temps à vous inquiéter de la commande de sortie avec les programmes MPI.

Le moyen le plus élégant d'y parvenir, si vous voulez vraiment, est de laisser les processus envoyer leurs messages à un rang, par exemple, classez 0 et laissez le rang 0 imprimer la sortie dans l'ordre qu'il a reçu ou commandé par rangs.

Encore une fois, ne passez pas trop de temps à essayer de commander la production des programmes MPI. Ce n'est pas pratique et a peu d'utilité.

0 commentaires
0
votes

Ajout aux réponses précédentes ici, votre MPI_Barrier fonctionne bien.

Cependant, si vous avez l'intention de le voir fonctionner, vous pouvez forcer une pause l'exécution ( Sleep (1) ) pour un instant pour que la sortie soit rattraper la sortie.

0 commentaires

Articles qui pourrait vous intéresser :

 Que puis-je utiliser pour la conversion d'entrée au lieu de scanf?
Comment puis-je calculer a * b / c lorsque a et b sont tous deux plus petits que c, mais que a * b déborde?
Comment vérifier si un pointeur void (void *) est l'un des deux types de données?
L'écran LCD affiche des caractères aléatoires (MPLAB)