def DSYR_blas(A, x, alpha=<NUM_LIT:1.>):
At = blas.dsyr(lower=<NUM_LIT:0>, x=x, a=A, alpha=alpha, overwrite_a=False) <EOL>A[:] = At<EOL>symmetrify(A, upper=True)<EOL>
Performs a symmetric rank-1 update operation:
A <- A + alpha * np.dot(x,x.T)

:param A: Symmetric NxN np.array
:param x: Nx1 np.array
:param alpha: scalar