[petsc-dev] VecAXPY() in loops in TS

Jed Brown jed at 59A2.org
Fri Feb 18 17:34:09 CST 2011

On Sat, Feb 19, 2011 at 00:24, Barry Smith <bsmith at mcs.anl.gov> wrote:

>  Even adjusting the array passed in will be worth it.

Yes, the work to futz around with the small array is tiny compared to the
cost of streaming the vector through extra times.

>  Should MAXPY have an optional constant scale factor, does it come up
> enough?

It looks like about half of the uses in PETSc are preceded by a loop to
apply some constant factor, so I think such a parameter would be worthwhile.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20110219/4eb751f4/attachment.html>

More information about the petsc-dev mailing list