Hi all,
This is not really a Fortran question, but maybe comebody can help me.
The most common vector update operation in numerical methods is the
axpy:
y<--y+alpha*x
So the BLAS provides this as a level 1 routine.
But what happens if I need (as I do in my project):
y<--beta*y+x
or the most general gaxpy:
y<--beta*y+alpha*x
Why aren't these in the BLAS? How can I emulate them most efficiently
with the level 1 routines?
Thanks,
Aleksandar
--
_____________________________________________
Aleksandar Donev
Physics Department
Michigan State University
East Lansing, MI 48824-1116
E-mail: [log in to unmask]
Work phone: (517) 432-6770
_____________________________________________
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|