Alvaro Fernandez wrote: > <snipped> > As for the gaxpy (again, my syntax may be transposed) > > call scal (beta,y) > call axpy(alpha,x) > where alpha is now some number not necessarily = 1.0. > > Alvaro Fernandez My goof: axpy should have both alpha, x, _and_ y as inputs, I think. (Like I said - it's been a while.) But the point is to use scal first and then axpy, either with alpha = 1.0 in the first case or alpha = <whatever> in the second. Alvaro Fernandez %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%