[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