[petsc-dev] LSC question

Barry Smith bsmith at mcs.anl.gov
Thu Mar 8 15:23:41 CST 2012



  Jed,

     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?

   Barry




More information about the petsc-dev mailing list