<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Oct 2, 2017 at 6:21 AM, Florian Lindner <span dir="ltr"><<a href="mailto:mailinglists@xgm.de" target="_blank">mailinglists@xgm.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
<br>
I have a matrix and vector that live on 4 ranks, but only rank 2 and 3 have values:<br>
<br>
e.g.<br>
<br>
Vec Object: 4 MPI processes<br>
  type: mpi<br>
Process [0]<br>
Process [1]<br>
1.1<br>
2.5<br>
3.<br>
4.<br>
Process [2]<br>
5.<br>
6.<br>
7.<br>
8.<br>
Process [3]<br>
<br>
<br>
Doing a simple LSQR solve does not converge. However, when the values are distributed equally, it converges within 3<br>
iterations.<br>
<br>
What can I do about that?<br>
<br>
I have attached a simple program and creates the matrix and vector or loads them from a file.<br></blockquote><div><br></div><div>There are a few problems with this program. I am attaching a cleaned up version. However, convergence still differs starting</div><div>at iteration 2. It appears that LSQR has a problem with this system, or we have a bug that I cannot see.</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">
Thanks,<br>
Florian<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><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="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>