On Tue, Oct 30, 2012 at 4:33 PM, Mark F. Adams <span dir="ltr"><<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word"><br><div><div>On Oct 25, 2012, at 2:03 PM, Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>> wrote:</div><br><blockquote type="cite">On Thu, Oct 25, 2012 at 12:52 PM, Mark F. Adams <span dir="ltr"><<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>></span> wrote:<br>
<div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
How should I go about getting -ksp_monitor to report the infinity norm?<br></blockquote><div><br></div><div>Preconditioned or unpreconditioned residual?</div><div><br></div><div>You would use KSPMonitorSet(). This may be handy enough that PETSc should have -ksp_monitor_max. </div>
</div></blockquote><div><br></div><div>OK, I'm going to start putting in -ksp_monitor_max, speak now or forever hold your peace.</div><br><blockquote type="cite"><div class="gmail_quote"><div>You can see the setup code in KSPSetFromOptions(). You could copy KSPMonitorTrueResidualNorm() or refactor so that both call a common helper that takes an additional argument (the norm type). Or put the norm type in the context data structure, I guess that's what it's there for, but then you need a destructor.</div>
</div></blockquote></div></div></blockquote></div>We are getting enough of these that they might make more sense as flags. Right now, true_residual prints<div>the true and preconditioned residual. I think it would be nicer to print all the different norms I want on one line.</div>
<div><br></div><div>   Matt<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<br>
</div>