Re: [petsc-users] snes matrix-free jacobian fails with preconditioner shape mismatch

2024-03-10 Thread Barry Smith
> On Mar 10, 2024, at 10:16 AM, Yi Hu wrote: > > This Message Is From an External Sender > This message came from outside your organization. > Dear Mark, > > Thanks for your reply. I see this mismatch. In fact my global DoF is 324. It > seems like I always get the local size = global Dof /

Re: [petsc-users] snes matrix-free jacobian fails with preconditioner shape mismatch

2024-03-10 Thread Yi Hu
Dear Mark, Thanks for your reply. I see this mismatch. In fact my global DoF is 324. It seems like I always get the local size = global Dof / np^2, np is my processor number. By the way, I used DMDASNESsetFunctionLocal() to set my form function. Is it eligible to mix

Re: [petsc-users] snes matrix-free jacobian fails with preconditioner shape mismatch

2024-03-10 Thread Mark Adams
It looks like your input vector is the global vector, size 162, and the local matrix size is 81. Mark [1]PETSC ERROR: Preconditioner number of local rows 81 does not equal input vector size 162 On Sun, Mar 10, 2024 at 7:21 AM Yi Hu wrote: > Dear petsc team, I implemented a matrix-free

[petsc-users] snes matrix-free jacobian fails with preconditioner shape mismatch

2024-03-10 Thread Yi Hu
Dear petsc team, I implemented a matrix-free jacobian, and it can run sequentially. But running parallel I got the pc error like this (running with mpirun -np 2, only error from rank1 is presented here) [1]PETSC ERROR: - ZjQcmQRYFpfptBannerStart