[petsc-dev] Note about CPR preconditioner in PCFIELDSPLIT man page?

2019-02-20 Thread Mills, Richard Tran via petsc-dev
Folks, There is a note in the PCFIELDSPLIT manual page about implementing a so-called "Constrained Pressure Residual" (CPR) preconditioner: "The Constrained Pressure Preconditioner (CPR) can be implemented using

[petsc-dev] Errors in MatSetValuesBlockedLocal() when using AIJ, but not BAIJ!

2019-02-20 Thread Mills, Richard Tran via petsc-dev
Folks, I'm working with some PFLOTRAN examples, and I've been scratching my head for a while at errors that are being generated in MatSetValuesBlockedLocal(), related to new nonzeros requiring a malloc. E.g., [0]PETSC ERROR: Argument out of range [0]PETSC ERROR: New nonzero at (2,0) caused a

Re: [petsc-dev] Regression in Vec assembly

2019-02-20 Thread Jed Brown via petsc-dev
https://bitbucket.org/petsc/petsc/pull-requests/1389/vecstashsortcompress_private-fix/diff Stefano Zampini writes: > Jed, > > thanks. It fixed my issue too. > > > Il giorno mer 20 feb 2019 alle ore 09:33 Lisandro Dalcin > ha scritto: > >> I tested your fix (on top of master) with my original

Re: [petsc-dev] Regression in Vec assembly

2019-02-20 Thread Stefano Zampini via petsc-dev
Jed, thanks. It fixed my issue too. Il giorno mer 20 feb 2019 alle ore 09:33 Lisandro Dalcin ha scritto: > I tested your fix (on top of master) with my original PetIGA reproducer, > and the issue seems to be fixed. > > On Wed, 20 Feb 2019 at 08:32, Jed Brown wrote: > >> Yuck, bad bug. I