Improving the efficiency of portable software for linear algebra

In algorithms for linear algebraic computations there are a fairly small number of basic operations which are generally responsible for a significant percentage of the total execution time. We mention particularly the dot product operation, ω := x T y, the elementary vector operation, y :=αx+y, and the euclidean vector norm, η = (x T x) 1/2 .