Function rgsl::blas::level2::dsyr2

pub fn dsyr2(
    uplo: Uplo,
    alpha: f64,
    x: &VectorF64,
    y: &VectorF64,
    A: &mut MatrixF64
) -> Value

These functions compute the symmetric rank-2 update A = \alpha x y^T + \alpha y x^T + A of the symmetric matrix A. Since the matrix A is symmetric only its upper half or lower half need to be stored. When Uplo is CblasUpper then the upper triangle and diagonal of A are used, and when Uplo is CblasLower then the lower triangle and diagonal of A are used.