Function rgsl::blas::level3::dsyr2k

pub fn dsyr2k(
    uplo: Uplo,
    trans: Transpose,
    alpha: f64,
    A: &MatrixF64,
    B: &MatrixF64,
    beta: f64,
    C: &mut MatrixF64
) -> Value

This function computes a rank-2k update of the symmetric matrix C, C = \alpha A B^T + \alpha B A^T + \beta C when Trans is NoTrans and C = \alpha A^T B + \alpha B^T A + \beta C when Trans is Trans. Since the matrix C is symmetric only its upper half or lower half need to be stored. When Uplo is Upper then the upper triangle and diagonal of C are used, and when Uplo is Lower then the lower triangle and diagonal of C are used.