<div class="gmail_extra">On Sun, Apr 22, 2012 at 3:55 PM, TAY wee-beng <span dir="ltr"><<a href="mailto:zonexo@gmail.com" target="_blank">zonexo@gmail.com</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 bgcolor="#FFFFFF" text="#000000">
Hi, <br>
<br>
I changed the vs2008 and linux source code to
VecView(x,PETSC_VIEWER_STDOUT_WORLD,ierr) but the results for both
outputs are 0 and 0.0416667 (except for one value which is
0.0833333) respectively.</div></blockquote><div><br></div><div>It looks to me like there is a bug in KSPSetUp(). The DMKSPGetContext() is returning a null struct,</div><div>so the RHS is never formed, and the Mat is formed by the default call to DMComputeJacobian. What</div>
<div>is intended here with KSPDM?</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 bgcolor="#FFFFFF" text="#000000">
Yours sincerely,<pre cols="72">TAY wee-beng</pre>
<br>
On 22/4/2012 9:35 PM, Jed Brown wrote:
<blockquote type="cite">
<div class="gmail_quote">On Sun, Apr 22, 2012 at 14:30, TAY
wee-beng <span dir="ltr"><<a href="mailto:zonexo@gmail.com" target="_blank">zonexo@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 bgcolor="#FFFFFF" text="#000000">Sorry for the mistake. I
tried again, changing to "call VecView(x,0,ierr)"<br>
<br>
but still got the same error:</div>
</blockquote>
</div>
<br>
<div>0 has the wrong type. Use
VecView(x,PETSC_VIEWER_STDOUT_WORLD,ierr)</div>
</blockquote>
</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<br>
</div>