[petsc-dev] Lies in PCFieldSplit
Matthew Knepley
knepley at gmail.com
Sun Dec 2 23:00:45 CST 2012
On Sun, Dec 2, 2012 at 10:57 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
>
> Jed, Matt, or I are responsible for what are likely lies in the field split naming terminology
>
> We have
>
> const char *const PCFieldSplitSchurPreTypes[] = {"SELF","DIAG","USER","PCFieldSplitSchurPreType","PC_FIELDSPLIT_SCHUR_PRE_",0};
>
> Diag here actually means the A11 matrix, correct?
>
> case PC_FIELDSPLIT_SCHUR_PRE_DIAG:
> ierr = PetscViewerASCIIPrintf(viewer," Preconditioner for the Schur complement formed from the block diagonal part of A11\n");CHKERRQ(ierr);break;
> case PC_FIELDSPLIT_SCHUR_PRE_USER:
> if (jac->schur_user) {
>
> block diagonal part of A11 actually means A11, correct? Otherwise what blocks????
Processor diagonal. This is not ASM where what we select is
automatically process diagonal.
Matt
> Once I get confirmation that these are phrased wrong I will update the source to be less confusing. Otherwise please example what diagonal means and how it is used here.
>
> Thanks
>
> Barry
>
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which
their experiments lead.
-- Norbert Wiener
More information about the petsc-dev
mailing list