<div dir="ltr"><div><div>Hi Barry,<br><br></div>So for the matrix C that I am currently testing (size 162x162), the condition number is roughly 10^4.<br> <br>For reference, I'm porting MATLAB code into PETSc, and for one 
processor, the PETSc b vector is roughly equivalent to the MATLAB b vector. So I
 know that for one processor, my program is performing as expected.<br><br></div><div>I've included examples below of values for b (also of size 162), ranging from indices 131 to 141.<br><br></div><div>#processors=1: <br>
                         0<br>     1.315217173959314e-20<br>     1.315217173959314e-20<br>     4.843201487740107e-17<br>     4.843201487740107e-17<br>     8.166104700666665e-14<br>     8.166104700666665e-14<br>     6.303834267553249e-11<br>
     6.303834267553249e-11<br>     2.227932688485483e-08<br>     2.227932688485483e-08<br></div><div><br># processors=2:<br>     5.480410831461926e-22<br>     2.892553944350444e-22<br>     2.892553944350444e-22<br>     7.524038923310717e-24<br>
     7.524038923214420e-24<br>    -3.340766769043093e-26<br>    -7.558372155761972e-27<br>     5.551561288838557e-25<br>     5.550551546879874e-25<br>    -1.579397982093437e-22<br>     2.655766754178065e-22<br></div><div>
<br># processors = 4:<br>     5.480410831461926e-22<br>     2.892553944351728e-22<br>     2.892553944351728e-22<br>     7.524092205125593e-24<br>     7.524092205125593e-24<br>    -2.584939414228212e-26<br>    -2.584939414228212e-26<br>
                         0<br>                         0<br>    -1.245940797657998e-23<br>    -1.245940797657998e-23<br></div><div><br># processors = 8:<br>     5.480410831461926e-22<br>     2.892553944023035e-22<br>     2.892553944023035e-22<br>
     7.524065744581494e-24<br>     7.524065744581494e-24<br>    -2.250265175188197e-26<br>    -2.250265175188197e-26<br>    -6.543127892265160e-26<br>    1.544288143499193e-317<br>     8.788794008375919e-25<br>     8.788794008375919e-25<br>
<br><br></div><div>Thanks,<br>Arthur<br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jun 25, 2014 at 4:06 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
   How different are the values in b? Can you send back a few examples of the different b’s? Any idea of the condition number of C?<br>
<span class="HOEnZb"><font color="#888888"><br>
   Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
On Jun 25, 2014, at 3:10 PM, Arthur Kurlej <<a href="mailto:akurlej@gmail.com">akurlej@gmail.com</a>> wrote:<br>
<br>
> Hi all,<br>
><br>
> While running my code, I have found that MatMult() returns different values depending on the number of processors I use (and there is quite the variance in the values).<br>
><br>
> The setup of my code is as follows (I can go into more depth/background if needed):<br>
> -Generate parallel AIJ matrix of size NxN, denoted as A<br>
> -Retrieve parallel AIJ submatrix from the last N-1 rows&columns from A, denoted as C<br>
> -Generate vector of length N-1, denoted as x<br>
> -Find C*x=b<br>
><br>
> I have already checked that A, C, and x are all equivalent when ran for any number of processors, it is only the values of vector b that varies.<br>
><br>
> Does anyone have an idea about what's going on?<br>
><br>
><br>
> Thanks,<br>
> Arthur<br>
><br>
<br>
</div></div></blockquote></div><br></div>