Thank You very much! I don't know Why I was avoiding the ownprecond variant. May be just was scared to try. I think I'll use it. I also could just set the maximum number of iterations and watch the -ksp_monitor_true_residual, but it's not so convinient. Thank You!<br>
<br><div class="gmail_quote">2011/7/17 Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"><br>
On Jul 16, 2011, at 4:55 AM, áÌÅËÓÅÊ òÑÚÁÎÏ× wrote:<br>
<br>
> Hi, Im solving a linear system with different methods for comparison goals.<br>
> And I want to make Krylov iterative methods work until they reach the predefined value of the TRUE residual norm.<br>
> So Im setting up tolerances (in fact only atol) with KSPSetTolerances(KSP, 1.e-50, 1.e-10, PETSC_DEFAULT, PETSC_DEFAULT) and watching for -ksp_monitor (or -ksp_monitor_true_residual).<br>
> By default most of methods monitor their preconditioned residual norm and stop working then the preconditioned residual norm reaches the predefined value.<br>
> That of course doesn't mean true residual norm also reached the value.<br>
><br>
><br>
> So I need to change this behavior and use KSPSetNormType(KSP, KSP_NORM_UNPRECONDITIONED) to use unpreconditioned residual or KSPSetPreconditionerSide(KSP, PC_RIGHT) to make preconditioned residual be the same as unpreconditioned and achieve my goal anyway.<br>
><br>
> So what I have:<br>
><br>
><br>
> 1)<br>
> KSPSetNormType(KSP, KSP_NORM_UNPRECONDITIONED) is only supported by CG, RICHARDSON AND CHEBYCHEV (a said in user manual).<br>
> So as it is expected, CG, RICHARDSON AND CHEBYCHEV print error message: "No right preconditioning for KSPCG, KSPRICHARDSON, KSPCHEBYCHEV!" and work šgreat with KSPSetNormType(KSP, KSP_NORM_UNPRECONDITIONED)<br>
> Great!<br>
><br>
> 2)<br>
> Practically the same behavior does KSPLSQR (except the great work, actually it doesn't work at all with my linear system, but I don't care about lsqr)<br>
> So it's ok.<br>
><br>
><br>
> 3)<br>
> KSPBCGS and KSPBICG don't support right preconditioning, so their true residual norm values cant be set as the finish point for iterating.<br>
> Am I right?<br>
<br>
</div> š šActually in petsc-dev <a href="http://www.mcs.anl.gov/petsc/petsc-as/developers/index.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/developers/index.html</a> KSPBCGS does support both left and right preconditioning. KSPBICG is coded only for left preconditioning<br>
<div class="im">><br>
><br>
> 4)<br>
> KSPGMRES print error message: "Use right preconditioning -ksp_right_pc if want unpreconditioned norm", when I use it with KSPSetNormType(KSP, KSP_NORM_UNPRECONDITIONED).<br>
> And works great with šKSPSetPreconditionerSide(KSP, PC_RIGHT).<br>
> As expected.<br>
><br>
> 5)<br>
> KSPCGS prints no error, when I use it with KSPSetNormType(KSP, KSP_NORM_UNPRECONDITIONED) and works with preconditioned norm.<br>
> And works great with šKSPSetPreconditionerSide(KSP, PC_RIGHT).<br>
> So it's ok.<br>
><br>
> 6)<br>
> KSPTFQMR and KSPTCQMR print no errors with any settings and works with preconditioned norm.<br>
> Is there any methods to make KSPTFQMR use unpreconditioned norm or right preconditioning?<br>
<br>
</div> š These two algorithms actually do work with both left and right preconditioner but they are different than all the other methods in that:<br>
<br>
š šThe "residual norm" computed in this algorithm is actually just an upper bound on the actual residual norm.<br>
š š š š šThat is for left preconditioning it is a bound on the preconditioned residual and for right preconditioning<br>
š š š š šit is a bound on the true residual.<br>
<br>
š šI have added this information to their manual pages in petsc-dev for clarity.<br>
<br>
š š Note that our GMRES is implemented only with left preconditioning and preconditioned residual norm but you can use the FGMRES which is implemented only for right preconditioning and the true residual norm for testing.<br>
<br>
<br>
š šIf you not particularly interested in speed and just want the iterations to stop when the true residual norm reaches some tolerance you can easily provide your own KSPConvergedTest() by calling KSPSetConvergenceTest() and writing a routine that calls KSPBuildResidual() (which always builds the true residual) and compute the norm of the result.<br>
<font color="#888888"><br>
š šBarry<br>
</font><div><div></div><div class="h5"><br>
<br>
><br>
><br>
> Please help me to understand situations 3) and 6)<br>
><br>
> Thank in advance!<br>
><br>
> --<br>
> Best regards,<br>
> Alexey Ryazanov<br>
> ______________________________________<br>
> Nuclear Safety Institute of Russian Academy of Sciences<br>
><br>
><br>
<br>
</div></div></blockquote></div><br>