MatMult_Scatter woes

[Jed Brown]

On Tue 2009-04-07 21:02, Barry Smith wrote:
>
>    Jed,
>
>     Please remind me of what the simplest fix is? Can we just zero the 
> vector before the scatter for the MatMult and MatMultTranspose?
> What about for MatMultAdd and MatMultTransposeAdd, does anything need to 
> be done?

The simplest fix is to zero the vector before and always use ADD_VALUES.
I pushed this a week ago and cited this thread in the commit message,
see 6d3d60abcf7d.  It's just worth remembering that it can be optimized
if someone ever finds that this is a bottleneck.

Jed
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 197 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20090408/ff384fb0/attachment.sig>