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^Ty, the elementary vector operation, y :=αx+y, and the euclidean vector norm, η = (x^Tx)^1/2.