[petsc-dev] PETSc communicators with MPI in Fortran

Satish Balay balay at mcs.anl.gov
Tue Oct 31 20:59:30 CDT 2017


You should be able to use any MPI communicator with PETSc objects.

PETSC_COMM_WORLD/PETSC_COMM_SELF are just a couple of convinence
communicators PETSc provides.

Satish

On Tue, 31 Oct 2017, Neelam Patel wrote:

> Hello PETSc users,
> Working in Fortran, I created 2 disjoint communicators with MPI_Group operations using PETSC_COMM_WORLD as the "base" comm. I created parallel vectors on each communicator, and set values in them equal to their ranks on PETSC_COMM_WORLD. Everything seems to be working fine. My question is: Is this a valid thing to do? Or will it cause some undefined behaviour that I'm not seeing right now?
> I somehow thought that communicators created this way were not 'visible' to PETSc. Or is it that using PETSC_COMM_WORLD above eliminates that problem?
> Thank you in advance,Neelam
> 



More information about the petsc-dev mailing list