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

[Dominik Szczerba]

On Fri, Oct 21, 2011 at 11:21 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
> On Fri, Oct 21, 2011 at 16:14, Dominik Szczerba <dominik at itis.ethz.ch>
> wrote:
>>
>> Hmm, I just want to partition my mesh that I read in serial. So
>> adjacency vectors are only on root. How can I provide distributed
>> adjacency information, before partitioning? This is an endless circle.
>> Or do you mean I should just broadcast my adjacency vectors ia,ja to
>> all the other processes?
>
> The usual thing is to distribute the mesh naively to begin with, then
> partition and move the mesh to the correct place with respect to the new
> partition. It is possible that the case you describe works reliably, in
> which case the guard was to blunt. The crashes I'm familiar with occur when
> there are very few nodes such that some processors don't get any. You can
> open up pmetis.c and remove the guard entirely or write a better guard.

I can not find anything resembling any guards in pmetis.c as pet Petsc
3.2. Could you please refer me to the specific section?

Thanks a lot,
Dominik