Function dspr

matrix-vector operation defined as A := alpha * x * x' + A where A is a symmetric matrix, supplied in packed form