<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sun, Mar 29, 2015 at 11:32 AM, gouarin <span dir="ltr"><<a href="mailto:loic.gouarin@math.u-psud.fr" target="_blank">loic.gouarin@math.u-psud.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hi,<br>
<br>
I think that there is a bug in the MatMultAdd_Shell function when entries y and z are the same Vec. It is not possible to call<br>
<br>
VecWAXPY(z,1.0,shell->right_<u></u>add_work,y);<br>
<br>
if z and y are the same Vec.<br>
<br>
I have attached a patch.<br></blockquote><div><br></div><div>Good catch, but this has been fixed in master:</div><div><br></div><div>  <a href="https://bitbucket.org/petsc/petsc/src/23524bc85f88d7c2c867e276c52fc9d9816c7ba0/src/mat/impls/shell/shell.c?at=master#cl-259">https://bitbucket.org/petsc/petsc/src/23524bc85f88d7c2c867e276c52fc9d9816c7ba0/src/mat/impls/shell/shell.c?at=master#cl-259</a></div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
Best,<br>
Loic<span class=""><font color="#888888"><br>
<br>
-- <br>
Loic Gouarin<br>
Laboratoire de Mathématiques<br>
Université Paris-Sud<br>
Bâtiment 425<br>
91405 Orsay Cedex<br>
France<br>
Tel: <a href="tel:%28%2B33%29%201%2069%2015%2060%2014" value="+33169156014" target="_blank">(+33) 1 69 15 60 14</a><br>
Fax: <a href="tel:%28%2B33%29%201%2069%2015%2067%2018" value="+33169156718" target="_blank">(+33) 1 69 15 67 18</a><br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>