<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Apr 30, 2018 at 1:38 AM, Stefano Zampini <span dir="ltr"><<a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><br>
> On Apr 30, 2018, at 3:14 AM, Smith, Barry F. <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br>
> <br>
> <br>
>   I don't understand<br>
> <br>
>    The code has<br>
> <br>
>  ierr = MatGetDiagonal(A,W2);CHKERRQ(<wbr>ierr);<br>
>  ierr = VecView(W2,viewer);CHKERRQ(<wbr>ierr);<br>
> <br>
>   the output from this VecView seems to match the diagonal of the matrix printed earlier.<br>
> <br>
>> On Apr 29, 2018, at 10:02 AM, Stefano Zampini <<a href="mailto:stefano.zampini@gmail.com">stefano.zampini@gmail.com</a>> wrote:<br>
>> <br>
>> By looking at src/mat/examples/tests/output/<wbr>ex88_1.out, it seems that the output of MatGetDiagonal is wrong for MATSHELL. <br>
>> <br>
>> it should be <br>
>> 2.96782e+08<br>
>> 1.41731e+09<br>
>> and not (as computed) <br>
>> 1.10918e+08<br>
>> 2.06459e+08<br>
> <br>
>   What do you mean "as computed"? The output file does not contain<br>
>> 1.10918e+08<br>
>> 2.06459e+08<br>
> anywhere.<br>
> <br>
>   Are you modifying the example and running it again to get the incorrect value?<br>
<br>
</span>Yes, my bad. I did not explain it correctly<br>
>From maint, if you just run ex88 you will get the incorrect output, which is not caught in the nightlyies because floating point values are completely ignored.<br></blockquote><div><br></div><div>Now that we ignore all values, tests need to be changed to check the distance from the correct vector.</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">
<br>
[szampini@localhost tests]$ ./ex88<br>
Matrix of type: shell<br>
Mat Object: 1 MPI processes<br>
  type: seqdense<br>
2.9678190300000000e+08 3.9288008100000000e+08 <br>
6.1374683500000000e+08 1.4173141580000000e+09 <br>
Vec Object: 1 MPI processes<br>
  type: seq<br>
3.92973e+19<br>
8.54197e+19<br>
Vec Object: 1 MPI processes<br>
  type: seq<br>
1.10918e+08<br>
2.06459e+08<br>
MatDiagonalSet not tested on MATSHELL<br>
….<br>
<br>
> <br>
>  Barry<br>
> <br>
<span class="HOEnZb"><font color="#888888">>> <br>
>> -- <br>
>> Stefano<br>
> <br>
<br>
</font></span></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>