4{
5 MPI_Comm bcomm, ucomm;
6 MPI_Group bgrp, ugrp;
8
10 for (i=0; i<Np; i++) mpmap[i] = pmap[i];
11
12 ucomm = MPI_Comm_f2c(Ucomm);
13 i=MPI_Comm_group(ucomm, &ugrp);
14 i=MPI_Group_incl(ugrp, Np, mpmap, &bgrp);
15 i=MPI_Comm_create(ucomm, bgrp, &bcomm);
16 i=MPI_Group_free(&ugrp);
17 i=MPI_Group_free(&bgrp);
18
19 free(mpmap);
20
21 return(bcomm);
22}