[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