[petsc-users] KspTrueResidualNorm

[Barry Smith]

  Right preconditioner shouldn't matter to the preconditioner at all. What is the complete error message. Does it run on one process? 

   Barry

On Feb 1, 2011, at 2:32 PM, Ataollah Mesgarnejad wrote:

> Barry,
> 
> I did as you said and now I receive these errors once I try to run the program:
> 
> [3]PETSC ERROR: MatHYPRE_IJMatrixCreate() line 76 in src/dm/da/utils/mhyp.c
> [3]PETSC ERROR: PCSetUp_HYPRE() line 112 in src/ksp/pc/impls/hypre/hypre.c
> [3]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
> [3]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
> [3]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c
> [3]PETSC ERROR: UStep() line 511 in PFMAT-FD.cpp
> [4]PETSC ERROR: MatHYPRE_IJMatrixCreate() line 76 in src/dm/da/utils/mhyp.c
> [4]PETSC ERROR: PCSetUp_HYPRE() line 112 in src/ksp/pc/impls/hypre/hypre.c
> [4]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
> [4]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
> [3]PETSC ERROR: main() line 95 in PFMAT-main.cpp
> [4]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c
> [4]PETSC ERROR: UStep() line 511 in PFMAT-FD.cpp
> [4]PETSC ERROR: main() line 95 in PFMAT-main.cpp
> 
> Does right preconditiong work with HYPRE?
> Best,
> Ata
> 
> On Tue, Feb 1, 2011 at 2:21 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
> 
>   The simplest thing is to simply run with -ksp_monitor_true_residual and see how the convergence is going in the true residual norm also.
> 
>   You can also switch to right preconditioning with gmres and then the residual used by gmres is the true residual norm.  Use -ksp_preconditioner_side right -ksp_norm_type unpreconditioned
> 
>   Barry
> 
> 
> 
> On Feb 1, 2011, at 2:14 PM, Ataollah Mesgarnejad wrote:
> 
> > Dear all,
> >
> > I'm using Ksp gmres with boomeramg preconditioning and I suspect that even though it converges in preconditioned norm it doesn't converge in true norm, but as I understand KSPSetNormType gmres does not support true residual norm? Is that correct and If it is, is there any other way to monitor the true norm? Do I need to Introduce my own convergence test?
> >
> > Best,
> > A. Mesgarnejad
> 
> 
> 
> 
> -- 
> A. Mesgarnejad
> PhD Student, Research Assistant
> Mechanical Engineering Department
> Louisiana State University
> 2203 Patrick F. Taylor Hall
> Baton Rouge, La 70803