<div dir="ltr">Thanks, those functions are what I was looking for.<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sat, Mar 15, 2014 at 6:29 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div class="">On Sat, Mar 15, 2014 at 6:24 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br>

</div><div class=""><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><br>
<br>
   You can write a custom monitor and set it with TSMonitorSet()<br>
<br>
    This routine would call TSGetSNES() then SNESGetSolution() then call SNESComputeFunction() then call VecView() on the result.<br>
<br>
    But note that just because the residual is big somewhere doesn’t mean the error need be.<br>
<br>
    You could also run with -snes_monitor_residual to see how the residual is being reduced inside the nonlinear solve (that is, what parts of the residual are most stubborn).</blockquote><div><br></div></div><div>If you really want to play with the residual, inside your monitor you can use:</div>

<div><br></div><div>  TSGetSNES()</div><div>  <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/SNESGetFunction.html#SNESGetFunction" target="_blank">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/SNESGetFunction.html#SNESGetFunction</a></div>

<div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div class=""><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">

<span><font color="#888888"><br>
    Barry<br>
</font></span><div><div><br>
<br>
On Mar 15, 2014, at 5:47 PM, Mani Chandra <<a href="mailto:mc0710@gmail.com" target="_blank">mc0710@gmail.com</a>> wrote:<br>
<br>
> Hi,<br>
><br>
> Is there anyway I can VecView the residual after TS has completed an implicit time step? I'd like to see where in my domain most of the errors are coming from. I looked at TSMonitor but that doesn't seem to give me access to the residual at the end of the current time step.<br>


><br>
> Thanks,<br>
> Mani<br>
<br>
</div></div></blockquote></div></div><span class="HOEnZb"><font color="#888888"><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
</font></span></div></div>
</blockquote></div><br></div>