/// <summary>
/// Performs a symmetric rank-2k update. C = alpha * A * B.T + alpha * B * A.T + beta * C.
/// </summary>
/// <param name="alpha"></param>
/// <param name="A"></param>
/// <param name="beta"></param>
/// <param name="C"></param>
public static void Symmetric_Rank_2k_Update(double alpha, vector A, vector B, double beta, matrix C)
{
Blas.syr2k(Layout.ColMajor, UpLo.Lower, Trans.No, A.Length, 1, alpha, A.Array, A.Length, B.Array, B.Length, beta, C.Array, C.Rows);
}
/// <summary>
/// Performs a symmetric rank-2k update. C = alpha * A * B.T + alpha * B * A.T + beta * C.
/// </summary>
/// <param name="alpha"></param>
/// <param name="A"></param>
/// <param name="beta"></param>
/// <param name="C"></param>
public static void Symmetric_Rank_2k_Update(double alpha, matrix A, matrix B, double beta, matrix C)
{
Blas.syr2k(Layout.ColMajor, UpLo.Lower, Trans.No, A.Rows, A.Cols, alpha, A.Array, A.Rows, B.Array, B.Rows, beta, C.Array, C.Rows);
}