<div dir="ltr">Hello again!<div><br></div><div>After some tweaking the code is giving right answers now. However it start to disagree with MATLAB results ('traditional' way using matrix inverse) when the grid is larger than 100*100. My PhD advisor and I suspects that the default dimension of the Krylov subspace is 100 in the test case we are running. If so, is there a way to increase the size of the subspace?</div><div><br></div><div><div><img src="cid:ii_jnccwbst0" alt="Disagrees.png" width="547" height="289"><br></div></div><div><br></div><div>Thanks!</div><div><br></div><div>Wang Weizhuo</div></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Oct 9, 2018 at 2:50 AM Mark Adams <<a href="mailto:mfadams@lbl.gov">mfadams@lbl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">To reiterate what Matt is saying, you seem to have the exact solution on a 10x10 grid. That makes no sense unless the solution can be represented exactly by your FE space (eg, u(x,y) = x + y).</div><br><div class="gmail_quote"><div dir="ltr">On Mon, Oct 8, 2018 at 9:33 PM Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Mon, Oct 8, 2018 at 9:28 PM Weizhuo Wang <<a href="mailto:weizhuo2@illinois.edu" target="_blank">weizhuo2@illinois.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>The code is attached in case anyone wants to take a look, I will try the high frequency scenario later.</div></div></blockquote><div><br></div><div>That is not the error. It is superconvergence at the vertices. The real solution is trigonometric, so your</div><div>linear interpolants or whatever you use is not going to get the right value in between mesh points. You</div><div>need to do a real integral over the whole interval to get the L_2 error.</div><div><br></div><div>  Thanks,</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"><div class="gmail_quote"><div dir="ltr">On Mon, Oct 8, 2018 at 7:58 PM Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br><br><div class="gmail_quote"><div dir="ltr">On Mon, Oct 8, 2018 at 6:58 PM Weizhuo Wang <<a href="mailto:weizhuo2@illinois.edu" target="_blank">weizhuo2@illinois.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr">The first plot is the norm with the flag -pc_type lu with respect to number of grids in one axis (n), and the second plot is the norm without the flag -pc_type lu. </div></div></div></blockquote><div><br></div><div>So you are using the default PC w/o LU. The default is ILU. This will reduce high frequency effectively but is not effective on the low frequency error. Don't expect your algebraic error reduction to be at the same scale as the residual reduction (what KSP measures). </div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
</blockquote></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="m_6685943040659866736m_3945712423757607757m_-3838818614766066975gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Wang Weizhuo</div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="m_6685943040659866736m_3945712423757607757gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.cse.buffalo.edu_-7Eknepley_&d=DwMFaQ&c=OCIEmEwdEq_aNlsP4fF3gFqSN-E3mlr2t9JcDdfOZag&r=hsLktHsuxNfF6zyuWGCN8x-6ghPYxhx4cV62Hya47oo&m=EFM29ATgv4U8PjXEtfgMkuxKr5DGscMlH-j769W5W_4&s=grgSL2LaDCthvYvvFITmeOOWPCwgmNfYRPs94N8kmOs&e=" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>
</blockquote></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Wang Weizhuo</div></div>