[petsc-dev] bjacobi bugs ?

Barry Smith bsmith at mcs.anl.gov
Wed Jan 25 11:34:32 CST 2012


   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




More information about the petsc-dev mailing list