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
