[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