Barry,<br><br>I did as you said and now I receive these errors once I try to run the program:<br><br>[3]PETSC ERROR: MatHYPRE_IJMatrixCreate() line 76 in src/dm/da/utils/mhyp.c<br>[3]PETSC ERROR: PCSetUp_HYPRE() line 112 in src/ksp/pc/impls/hypre/hypre.c<br>
[3]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c<br>[3]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c<br>[3]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c<br>[3]PETSC ERROR: UStep() line 511 in PFMAT-FD.cpp<br>
[4]PETSC ERROR: MatHYPRE_IJMatrixCreate() line 76 in src/dm/da/utils/mhyp.c<br>[4]PETSC ERROR: PCSetUp_HYPRE() line 112 in src/ksp/pc/impls/hypre/hypre.c<br>[4]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c<br>
[4]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c<br>[3]PETSC ERROR: main() line 95 in PFMAT-main.cpp<br>[4]PETSC ERROR: KSPSolve() line 353 in src/ksp/ksp/interface/itfunc.c<br>[4]PETSC ERROR: UStep() line 511 in PFMAT-FD.cpp<br>
[4]PETSC ERROR: main() line 95 in PFMAT-main.cpp<br><br>Does right preconditiong work with <b>HYPRE</b>?<br>Best,<br>Ata<br><br><div class="gmail_quote">On Tue, Feb 1, 2011 at 2:21 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><br>
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.<br>
<br>
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<br>
<font color="#888888"><br>
Barry<br>
</font><div><div></div><div class="h5"><br>
<br>
<br>
On Feb 1, 2011, at 2:14 PM, Ataollah Mesgarnejad wrote:<br>
<br>
> Dear all,<br>
><br>
> 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?<br>
><br>
> Best,<br>
> A. Mesgarnejad<br>
<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>A. Mesgarnejad<br>PhD Student, Research Assistant<br>Mechanical Engineering Department<br>Louisiana State University<br>2203 Patrick F. Taylor Hall<br>Baton Rouge, La 70803<br>