On Thu, Mar 8, 2012 at 15:23, Barry Smith <bsmith at mcs.anl.gov> wrote:
> ierr = PCLSCAllocate_Private(pc);CHKERRQ(ierr); > ierr = > PetscObjectQuery((PetscObject)pc->pmat,"LSC_L",(PetscObject*)&L);CHKERRQ(ierr); > ierr = > PetscObjectQuery((PetscObject)pc->pmat,"LSC_Lp",(PetscObject*)&Lp);CHKERRQ(ierr); > > > It uses pc->pmat in both places. Is this really what you want? Wouldn't > you sometimes want mat? > This is just attached data for an auxiliary operator that should have some properties that relate somehow to the original operator. The whole interface is dirty, but I don't know a better approach (although this one could have a proper API). You could put these auxiliary operators anywhere, but I think attaching both to the preconditioning operator is natural. I don't mind if you want it to look other places for these things. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120308/c44cca4e/attachment.html>
