<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Sep 24, 2014 at 2:31 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
  Brian,<br>
<br>
    Your convergence test routine could call KSPBuildResidual() and then compute the 1-norm of the residual and make any decision it likes based on that norm.  See KSPSetConvergenceTest()</blockquote><div><br></div><div>Note that what Barry is saying is that the convergence theory for CG guarantees monotonicity in the energy (A) norm, but says nothing</div><div>about L1 so you might get crap.</div><div><br></div><div>  Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="HOEnZb"><font color="#888888"><br>
   Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
On Sep 24, 2014, at 2:24 PM, Brian Yang <<a href="mailto:brianyang1106@gmail.com">brianyang1106@gmail.com</a>> wrote:<br>
<br>
> Thanks Mat and Barry,<br>
><br>
> Yes I want to use the L1 norm in the convergence test. For the ||Ax-b||, I always want to calculate L1 norm instead of L2. Is Mat's way gonna work on this?<br>
><br>
> On Wed, Sep 24, 2014 at 2:17 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br>
><br>
>   Brian,<br>
><br>
>    Do you wish to monitor the convergence of the L1 norm, or do you wish to use the L1 norm in the convergence test?   Regardless you cannot change the LSQR or CG to be a different algorithm (with different convergence properties) with the L1 norm.<br>
><br>
><br>
>   Barry<br>
><br>
> On Sep 24, 2014, at 2:08 PM, Brian Yang <<a href="mailto:brianyang1106@gmail.com">brianyang1106@gmail.com</a>> wrote:<br>
><br>
> > Hi all,<br>
> ><br>
> > For example, I am using LSQR or CG to solve Ax=b.<br>
> ><br>
> > When we turn on the monitor option of solving the linear system, we could see iteration number and residual for each iteration. I believe it's L2 norm of the objective function?<br>
> ><br>
> > If yes, is there a way that I could set it to L1 norm for solving the system?<br>
> ><br>
> > Thanks.<br>
> ><br>
> > --<br>
> > Brian Yang<br>
> > U of Houston<br>
> ><br>
> ><br>
> ><br>
><br>
><br>
><br>
><br>
> --<br>
> Brian Yang<br>
> U of Houston<br>
><br>
><br>
><br>
<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener
</div></div>