<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">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>
<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>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">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> </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 class=""><font color="#888888"><br>
Barry<br>
</font></span><div class=""><div class="h5"><br>
<br>
On Mar 15, 2014, at 5:47 PM, Mani Chandra <<a href="mailto:mc0710@gmail.com">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><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>