Create group communicator in mpi (mpi4py)












0















I am trying to create an MPI group using mpi4py. My current attempt looks like this:



from mpi4py import MPI

comm = MPI.COMM_WORLD
newGroup = comm.group.Excl([0, 1, 2])
print newGroup.size

newComm = comm.Create_group(newGroup)
print newComm
print newComm.Get_size()


The newGroup.size call does return 2 (I started the application with 5 processes) and the newComm variable indicates to be a communicator. But as soon as I try to call Get_size an exception is thrown:



mpi4py.MPI.Exception: MPI_ERR_COMM: invalid communicator



How can I create a new communicator based on a predefined sequence of process IDs?










share|improve this question



























    0















    I am trying to create an MPI group using mpi4py. My current attempt looks like this:



    from mpi4py import MPI

    comm = MPI.COMM_WORLD
    newGroup = comm.group.Excl([0, 1, 2])
    print newGroup.size

    newComm = comm.Create_group(newGroup)
    print newComm
    print newComm.Get_size()


    The newGroup.size call does return 2 (I started the application with 5 processes) and the newComm variable indicates to be a communicator. But as soon as I try to call Get_size an exception is thrown:



    mpi4py.MPI.Exception: MPI_ERR_COMM: invalid communicator



    How can I create a new communicator based on a predefined sequence of process IDs?










    share|improve this question

























      0












      0








      0








      I am trying to create an MPI group using mpi4py. My current attempt looks like this:



      from mpi4py import MPI

      comm = MPI.COMM_WORLD
      newGroup = comm.group.Excl([0, 1, 2])
      print newGroup.size

      newComm = comm.Create_group(newGroup)
      print newComm
      print newComm.Get_size()


      The newGroup.size call does return 2 (I started the application with 5 processes) and the newComm variable indicates to be a communicator. But as soon as I try to call Get_size an exception is thrown:



      mpi4py.MPI.Exception: MPI_ERR_COMM: invalid communicator



      How can I create a new communicator based on a predefined sequence of process IDs?










      share|improve this question














      I am trying to create an MPI group using mpi4py. My current attempt looks like this:



      from mpi4py import MPI

      comm = MPI.COMM_WORLD
      newGroup = comm.group.Excl([0, 1, 2])
      print newGroup.size

      newComm = comm.Create_group(newGroup)
      print newComm
      print newComm.Get_size()


      The newGroup.size call does return 2 (I started the application with 5 processes) and the newComm variable indicates to be a communicator. But as soon as I try to call Get_size an exception is thrown:



      mpi4py.MPI.Exception: MPI_ERR_COMM: invalid communicator



      How can I create a new communicator based on a predefined sequence of process IDs?







      python mpi mpi4py






      share|improve this question













      share|improve this question











      share|improve this question




      share|improve this question










      asked Nov 20 '18 at 22:05









      WoltanWoltan

      8,776126285




      8,776126285
























          1 Answer
          1






          active

          oldest

          votes


















          1














          newComm is a legit communicator on tasks [3-4], but is MPI_COMM_NULL on tasks [0-2]. The standard does not allow you to invoke MPI_Comm_size(MPI_COMM_NULL, ...), and hence your error.



          The solution is to call newComm.Get_size() explicitly on tasks [3-4], or on any tasks where newComm is not MPI_COMM_NULL.






          share|improve this answer























            Your Answer






            StackExchange.ifUsing("editor", function () {
            StackExchange.using("externalEditor", function () {
            StackExchange.using("snippets", function () {
            StackExchange.snippets.init();
            });
            });
            }, "code-snippets");

            StackExchange.ready(function() {
            var channelOptions = {
            tags: "".split(" "),
            id: "1"
            };
            initTagRenderer("".split(" "), "".split(" "), channelOptions);

            StackExchange.using("externalEditor", function() {
            // Have to fire editor after snippets, if snippets enabled
            if (StackExchange.settings.snippets.snippetsEnabled) {
            StackExchange.using("snippets", function() {
            createEditor();
            });
            }
            else {
            createEditor();
            }
            });

            function createEditor() {
            StackExchange.prepareEditor({
            heartbeatType: 'answer',
            autoActivateHeartbeat: false,
            convertImagesToLinks: true,
            noModals: true,
            showLowRepImageUploadWarning: true,
            reputationToPostImages: 10,
            bindNavPrevention: true,
            postfix: "",
            imageUploader: {
            brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
            contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
            allowUrls: true
            },
            onDemand: true,
            discardSelector: ".discard-answer"
            ,immediatelyShowMarkdownHelp:true
            });


            }
            });














            draft saved

            draft discarded


















            StackExchange.ready(
            function () {
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53402294%2fcreate-group-communicator-in-mpi-mpi4py%23new-answer', 'question_page');
            }
            );

            Post as a guest















            Required, but never shown

























            1 Answer
            1






            active

            oldest

            votes








            1 Answer
            1






            active

            oldest

            votes









            active

            oldest

            votes






            active

            oldest

            votes









            1














            newComm is a legit communicator on tasks [3-4], but is MPI_COMM_NULL on tasks [0-2]. The standard does not allow you to invoke MPI_Comm_size(MPI_COMM_NULL, ...), and hence your error.



            The solution is to call newComm.Get_size() explicitly on tasks [3-4], or on any tasks where newComm is not MPI_COMM_NULL.






            share|improve this answer




























              1














              newComm is a legit communicator on tasks [3-4], but is MPI_COMM_NULL on tasks [0-2]. The standard does not allow you to invoke MPI_Comm_size(MPI_COMM_NULL, ...), and hence your error.



              The solution is to call newComm.Get_size() explicitly on tasks [3-4], or on any tasks where newComm is not MPI_COMM_NULL.






              share|improve this answer


























                1












                1








                1







                newComm is a legit communicator on tasks [3-4], but is MPI_COMM_NULL on tasks [0-2]. The standard does not allow you to invoke MPI_Comm_size(MPI_COMM_NULL, ...), and hence your error.



                The solution is to call newComm.Get_size() explicitly on tasks [3-4], or on any tasks where newComm is not MPI_COMM_NULL.






                share|improve this answer













                newComm is a legit communicator on tasks [3-4], but is MPI_COMM_NULL on tasks [0-2]. The standard does not allow you to invoke MPI_Comm_size(MPI_COMM_NULL, ...), and hence your error.



                The solution is to call newComm.Get_size() explicitly on tasks [3-4], or on any tasks where newComm is not MPI_COMM_NULL.







                share|improve this answer












                share|improve this answer



                share|improve this answer










                answered Nov 21 '18 at 0:04









                Gilles GouaillardetGilles Gouaillardet

                4,52951323




                4,52951323
































                    draft saved

                    draft discarded




















































                    Thanks for contributing an answer to Stack Overflow!


                    • Please be sure to answer the question. Provide details and share your research!

                    But avoid



                    • Asking for help, clarification, or responding to other answers.

                    • Making statements based on opinion; back them up with references or personal experience.


                    To learn more, see our tips on writing great answers.




                    draft saved


                    draft discarded














                    StackExchange.ready(
                    function () {
                    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53402294%2fcreate-group-communicator-in-mpi-mpi4py%23new-answer', 'question_page');
                    }
                    );

                    Post as a guest















                    Required, but never shown





















































                    Required, but never shown














                    Required, but never shown












                    Required, but never shown







                    Required, but never shown

































                    Required, but never shown














                    Required, but never shown












                    Required, but never shown







                    Required, but never shown







                    Popular posts from this blog

                    Biblatex bibliography style without URLs when DOI exists (in Overleaf with Zotero bibliography)

                    ComboBox Display Member on multiple fields

                    Is it possible to collect Nectar points via Trainline?