Re: [petsc-users] Failure of MUMPS

2018-10-09 Thread Zhang, Junchao
OK, I found -ksp_error_if_not_converged will trigger PETSc to fail in this case. --Junchao Zhang On Tue, Oct 9, 2018 at 3:38 PM Junchao Zhang mailto:jczh...@mcs.anl.gov>> wrote: I met a case where MUMPS returned an out-of-memory code but PETSc continued to run. When PETSc calls MUMPS, it

Re: [petsc-users] Failure of MUMPS

2018-10-09 Thread Zhang, Junchao
I met a case where MUMPS returned an out-of-memory code but PETSc continued to run. When PETSc calls MUMPS, it checks if (A->erroriffailure). I added -mat_error_if_failure, but it did not work since it was overwritten by MatSetErrorIfFailure(pc->pmat,pc->erroriffailure) Does it suggest we

Re: [petsc-users] Problems about Block Jocabi and Matrix-Free Method in SNES

2018-10-09 Thread Matthew Knepley
On Mon, Oct 8, 2018 at 11:33 PM Yingjie Wu wrote: > Dear Petsc developer: > Hi, > > I've been studying Petsc recently about Precontioner and Metrix-Free, and > I have some questions that puzzle me. > > 1. I want to test block Jacobi preconditioner, so I chose > /snes/example/tutorial/ex3.c as an

Re: [petsc-users] Increasing norm with finer mesh

2018-10-09 Thread Mark Adams
To reiterate what Matt is saying, you seem to have the exact solution on a 10x10 grid. That makes no sense unless the solution can be represented exactly by your FE space (eg, u(x,y) = x + y). On Mon, Oct 8, 2018 at 9:33 PM Matthew Knepley wrote: > On Mon, Oct 8, 2018 at 9:28 PM Weizhuo Wang