1! 2! Program to test PetscSubcomm. 3! 4#include <petsc/finclude/petscsys.h> 5program main 6 use petscsys 7 implicit none 8 9 PetscErrorCode ierr 10 PetscSubcomm r 11 PetscMPIInt rank, size 12 MPI_Comm scomm 13 14 PetscCallA(PetscInitialize(ierr)) 15 PetscCallMPIA(MPI_Comm_size(PETSC_COMM_WORLD, size, ierr)) 16 PetscCallMPIA(MPI_Comm_rank(PETSC_COMM_WORLD, rank, ierr)) 17 PetscCallA(PetscSubcommCreate(PETSC_COMM_WORLD, r, ierr)) 18 PetscCallA(PetscSubcommSetFromOptions(r, ierr)) 19 PetscCallA(PetscSubcommSetTypeGeneral(r, rank, rank, ierr)) 20 21 PetscCallA(PetscSubcommGetChild(r, scomm, ierr)) 22 PetscCallA(PetscSubcommView(r, PETSC_VIEWER_STDOUT_WORLD, ierr)) 23 PetscCallA(PetscSubcommDestroy(r, ierr)) 24 PetscCallA(PetscFinalize(ierr)) 25end 26 27! 28!/*TEST 29! 30! test: 31! nsize: 2 32! 33!TEST*/ 34