In PCSetUp_BJacobi() it has at the end

     /* ------
     Setup code depends on the number of blocks
  */
  if (jac->n_local == 0) {
    ierr = PCSetUp_BJacobi_Singleblock(pc,mat,pmat);CHKERRQ(ierr);
  } else {
    ierr = PCSetUp_BJacobi_Multiblock(pc,mat,pmat);CHKERRQ(ierr);
  }

does this make sense? No. In the normal case it is not using 
PCSetUp_BJacobi_Singleblock() at all! Something has gone wrong in the logic of 
the code when multiblock was added.

    Barry


Reply via email to