[petsc-dev] bug in shell.c

[gouarin]

Hi,

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

VecWAXPY(z,1.0,shell->right_add_work,y);

if z and y are the same Vec.

I have attached a patch.

Best,
Loic

-- 
Loic Gouarin
Laboratoire de Mathématiques
Université Paris-Sud
Bâtiment 425
91405 Orsay Cedex
France
Tel: (+33) 1 69 15 60 14
Fax: (+33) 1 69 15 67 18

-------------- next part --------------
A non-text attachment was scrubbed...
Name: 0001-fix-a-bug-in-matmultadd_shell-when-y-and-z-are-the-s.patch
Type: text/x-patch
Size: 981 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20150329/78f6bb71/attachment.bin>