[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