Re: [petsc-users] question on PCLSC with matrix blocks of type MATNEST

2024-02-12 Thread Pierre Jolivet
> On 13 Feb 2024, at 12:33 AM, Zhang, Hong wrote: > > Pierre, > I just modified the error message of MatProductSymbolic() and added a testing > segment in src/mat/tests/ex195.c. I have not pushed my change yet. > > Your fix at >

Re: [petsc-users] question on PCLSC with matrix blocks of type MATNEST

2024-02-12 Thread Zhang, Hong via petsc-users
Pierre, I just modified the error message of MatProductSymbolic() and added a testing segment in src/mat/tests/ex195.c. I have not pushed my change yet. Your fix at https://gitlab.com/petsc/petsc/-/commit/9dcea022de3b0309e5c16b8c554ad9c85dea29cf?merge_request_iid=7283 is more general. Has this

Re: [petsc-users] Near nullspace lost in fieldsplit

2024-02-12 Thread Barry Smith
7279 does change the code for MatZeroRowsColumns_MPIAIJ(). But perhaps that does not resolve the problem you are seeing? If that is the case we will need a reproducible example so we can determine exactly what else is happening in your code to cause the difficulties. Here is the diff for

Re: [petsc-users] Near nullspace lost in fieldsplit

2024-02-12 Thread Matthew Knepley
On Mon, Feb 12, 2024 at 7:10 AM Jeremy Theler (External) < jeremy.theler-...@ansys.com> wrote: > Hi Barry > > > The bug fix for 2 is availabel in > https://gitlab.com/petsc/petsc/-/merge_requests/7279 > > Note that our code goes through MatZeroRowsColumns_MPIAIJ() not through >

Re: [petsc-users] Near nullspace lost in fieldsplit

2024-02-12 Thread Jeremy Theler (External)
Hi Barry > The bug fix for 2 is availabel in > https://gitlab.com/petsc/petsc/-/merge_requests/7279 Note that our code goes through MatZeroRowsColumns_MPIAIJ() not through MatZeroRows_MPIAIJ() so this fix does not change anything for the case I mentioned. -- jeremy