Function rgsl::blas::level3::dsyrk

pub fn dsyrk(
    uplo: CblasUplo, 
    trans: CblasTranspose, 
    alpha: f64, 
    A: &MatrixF64, 
    beta: f64, 
    C: &mut MatrixF64
) -> Value

This function computes a rank-k update of the symmetric matrix C, C = \alpha A A^T + \beta C when Trans is NoTrans and C = \alpha A^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.