[petsc-users] ERROR: MatPartitioningApply_Parmetis() ... Does not support any processor with 0 entries

[Dominik Szczerba]

The following used to work with Petsc 3.1, but fails in runtime at
MatPartitioningApply with Petsc 3.2 with the above error.
I do not seem to find any pointers in the Changelog. Can you please advise?

Regards,
Dominik

// Connectivity arrays ia, ja allocated and filled on master, zero otherwise
// mlocal is number of cells on master, zero otherwise
Mat adj;
ierr = MatCreateMPIAdj(PETSC_COMM_WORLD, mlocal, mlocal, ia, ja,
PETSC_NULL, &adj); CHKERRQ(ierr);
MatPartitioning part;
ierr = MatPartitioningCreate(PETSC_COMM_WORLD, &part); CHKERRQ(ierr);
ierr = MatPartitioningSetAdjacency(part, adj); CHKERRQ(ierr);
ierr = MatPartitioningSetNParts(part, np); CHKERRQ(ierr);
IS isAssignment;
ierr = MatPartitioningApply(part, &isAssignment); CHKERRQ(ierr);