Re: [petsc-users] Issue updating MUMPS ictnl after failed solve

2016-09-20 Thread Hong
David: I'll patch petsc-maint (v3.7), then merge it to petsc-master. It might take 1-2 days for regression tests. Hong On Tue, Sep 20, 2016 at 2:01 PM, Hong wrote: > >> David : >> This is a bug in PETSc. A change >> $ git diff ../../../pc/impls/factor/cholesky/cholesky.c >>

Re: [petsc-users] Issue updating MUMPS ictnl after failed solve

2016-09-20 Thread David Knezevic
On Tue, Sep 20, 2016 at 2:01 PM, Hong wrote: > David : > This is a bug in PETSc. A change > $ git diff ../../../pc/impls/factor/cholesky/cholesky.c > diff --git a/src/ksp/pc/impls/factor/cholesky/cholesky.c > b/src/ksp/pc/impls/factor/cholesky/cholesky.c > index

Re: [petsc-users] Issue updating MUMPS ictnl after failed solve

2016-09-20 Thread Hong
David : This is a bug in PETSc. A change $ git diff ../../../pc/impls/factor/cholesky/cholesky.c diff --git a/src/ksp/pc/impls/factor/cholesky/cholesky.c b/src/ksp/pc/impls/factor/cholesky/cholesky.c index 953d551..cc28369 100644 --- a/src/ksp/pc/impls/factor/cholesky/cholesky.c +++

Re: [petsc-users] Issue updating MUMPS ictnl after failed solve

2016-09-20 Thread David Knezevic
On Mon, Sep 19, 2016 at 11:04 PM, Hong wrote: > David : > I did following: > > PC pc; > Mat F; > ierr = KSPGetPC(ksp,);CHKERRQ(ierr); > ierr = PCReset(pc);CHKERRQ(ierr); > ierr = KSPSetOperators(ksp,A,A);CHKERRQ(ierr); > ierr =

Re: [petsc-users] Computing condition number

2016-09-20 Thread Matthew Knepley
On Tue, Sep 20, 2016 at 6:44 AM, Florian Lindner wrote: > Hello, > > to compute / approximate the condition number of a MATSBAIJ, I put > > -ksp_view > > # Conditon number estimate > -pc_type none > -ksp_type gmres > -ksp_monitor_singular_value > -ksp_gmres_restart 1000 > >

[petsc-users] Computing condition number

2016-09-20 Thread Florian Lindner
Hello, to compute / approximate the condition number of a MATSBAIJ, I put -ksp_view # Conditon number estimate -pc_type none -ksp_type gmres -ksp_monitor_singular_value -ksp_gmres_restart 1000 in my .petscrc Output is like: [...] 566 KSP Residual norm 1.241765807317e-07 % max