<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Message: 2<br>
Date: Tue, 24 Aug 2010 10:56:53 -0500<br>
From: Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>><br>
Subject: Re: [petsc-dev] Problem with conjugate gradient solver<br>
To: For users of the development version of PETSc<br>
        <<a href="mailto:petsc-dev@mcs.anl.gov">petsc-dev@mcs.anl.gov</a>>,        Victor Eijkhout <<a href="mailto:eijkhout@tacc.utexas.edu">eijkhout@tacc.utexas.edu</a>><br>
Message-ID: <<a href="mailto:A4C6A1ED-CDE7-4653-8046-C2797AD5AF31@mcs.anl.gov">A4C6A1ED-CDE7-4653-8046-C2797AD5AF31@mcs.anl.gov</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
<br>
   Well you are doing nothing wrong. When I run the code I get the same results as you. I dump the matrix into Matlab and it has a negative eigenvalue so the example is wrong. The code is handling it correctly.<br>
<br>
>> eig(full(Mat_0))<br>
<br>
ans =<br>
<br>
   -0.8042<br>
    0.1716<br>
    0.6489<br>
    1.0000<br>
    1.0000<br>
    1.0000<br>
    1.0000<br></blockquote><div><br>So are you saying that trying to solve the system of equations in the 
example will always result in divergence? I tried the same program with a
 different system that I know has a solution, and even that diverges. I get a proper solution if I use the default KSP solver.<br>
<br>
Anush<br> <br></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
   I am cc:iing Victor who supposedly wrote the original example.  <br></blockquote><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
   Barry<br>
<br>
Attached is the matrix saved with -mat_view_matlab<br></blockquote></div>