<div class="gmail_quote">On Mon, Jul 9, 2012 at 11:23 AM, Colin McAuliffe <span dir="ltr"><<a href="mailto:cjm2176@columbia.edu" target="_blank">cjm2176@columbia.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Ok thanks, I can see how introducing units to the residual is necessary, but when would it be necessary to introduce units to the state? I would imagine this would only be needed if a numerically differentiated jacobian was being used.<br>
</blockquote><div><br></div><div>It depends if you care about symmetry and if you want the step tolerance to be meaningful. It can also affect scaling of boundary conditions.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<br>
Thanks<span class="HOEnZb"><font color="#888888"><br>
Colin</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
Quoting Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>>:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
On Tue, Jun 26, 2012 at 2:31 PM, Colin McAuliffe <<a href="mailto:cjm2176@columbia.edu" target="_blank">cjm2176@columbia.edu</a>>wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello,<br>
<br>
I would like to use GMRES for a monolithic solution of a multiphysics<br>
problem. There are drastically different units in the problem and so the<br>
residual norm is dominated by the terms with the largest units. Is there a<br>
way to tell petsc to scale the residual before taking the norm to avoid<br>
this?<br>
<br>
</blockquote>
<br>
You can "fix" it with preconditioning, but then you can only use the<br>
preconditioned residual. You should really fix the units when formulating<br>
the problem. (This does not require formal nondimensionalization, just<br>
introduced units so that state and residuals are of order 1.)<br>
<br>
</blockquote>
<br>
<br>
<br></div></div><div class="HOEnZb"><div class="h5">
-- <br>
Colin McAuliffe<br>
PhD Candidate<br>
Columbia University<br>
Department of Civil Engineering and Engineering Mechanics<br>
</div></div></blockquote></div><br>