[petsc-dev] PCFieldSplit broken, what was intended?

Barry Smith bsmith at mcs.anl.gov
Fri May 27 12:56:09 CDT 2011


On May 27, 2011, at 11:44 AM, Jed Brown wrote:

> On Fri, May 27, 2011 at 18:35, Barry Smith <bsmith at mcs.anl.gov> wrote:
> I quick and dirty fix would be to add a reset flag into PC_Fieldsplit and then do this block of code only if that reset flag is true then turn off the reset flag.
> 
> What about putting that block of code in a different function and calling it directly from PCReset_FieldSplit()?

   I cannot be called from the PCReset() it needs to be called during the next PCSetUp(). The information is not available yet at the reset time.

   Barry

>  
> 
>     BTW: once we have PetscLayout containing field information we should totally rewrite the PCFIELDSPLIT to simplify and rationalize it,
>  
> Yes, absolutely.




More information about the petsc-dev mailing list